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FORWORD 


This report was prepared for NASA as part of contract NASl-10139 with 
Langley Research Center under the Advanced Applications Flight Experiments 
(AAFE) Program. The objective of this contract is the development of the Carbon 
Monoxide Pollution Experiment (COPE) . This experiment is designed primarily to 
obtain data for the investigation of mechanisms by which CO is removed from the 
earth's atmosphere. The approach uses an orbiting platform to remotely map 
global CO concentrations and determine vertical CO profiles using a correlation 
interferometer measurement technique developed by Barringer Research Ltd. The 
instrument is to be capable of measuring CO over the range of expected atmo- 
spheric burdens and of measuring other trace atmospheric constituents. 

Since various parts of the program have been reported in detail in other 
reports they are only briefly described to herein. These parts include the CO 
problem (ref. 1), the Breadboard (ref. 2), and the Engineering Model (ref. 3). 
The major topic covered in this report is that of flight measurements of carbon 
monoxide and methane . 
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EXPERIMENT DEFINITION 


1.1 CO in the Atmosphere 

Carbon monoxide (CO) is a toxic pollutant which is continually introduced 
into the Earth's atmosphere in significant quantities. It is the most widely 
distributed and most commonly occurring air pollutant. Although in recent years 
the effort to measure CO in the atmosphere has increased significantly, there is 
still a need for a much expanded program of CO measurements both in urban and 
rural areas and on a more global scale. Such data are needed because of the 
need to have sufficient data to evaluate pollution hazards as such, to quantita- 
tively evaluate sources of CO, to determine transport effects and to investigate 
sink mechanisms. That is, data are needed both for the investigation of local- 
ized short-term problems and for large-scale long-term effects. 

One of the latter questions is that of the CO "sink". There are apparently 
some mechanisms operating which destroy most of the CO in the atmosphere, i.e., 
a CO "sink". These mechanisms have not as yet been established in a quantita- 
tive sense. If there were no mechanism of destruction, then the CO concentra- 
tions would be expected to double in two or three years. This apparently is not 
occurring at present, but it is impossible to extrapolate CO concentrations 
without a knowledge of the sink mechanism. Thus, the potential of CO as an eco- 
logical problem in the future can not be determined. There is, therefore, sub- 
stantial interest in ascertaining the global balance of CO. Such a considera- 
tion involves a knowledge of the sources, the present global concentrations and 
the sinks of CO. The CO balance in the atmosphere has been discussed in detail 
elsewhere (ref. 1) . 

The problem of carbon monoxide in the atmosphere is still replete with un- 
settled questions. Although over the past several years significant new infor- 
mation has been obtained relating to sinks, to sources, and to concentrations, 
none of these questions have been definitely and quantitatively answered. The 
present status of information concerning the sink question indicates that the 
most likely sinks are interaction with soil bacteria, reaction with OH near the 
ground, and transport to and reaction with OH in the stratosphere. With these 
singly, or in combination, lifetimes of a few tenths of a year could be ac- 
counted for although for the definite establishment of a sink mechanism both 
better data on various parameters involved and better CO concentration data are 
needed. 

Sinks are most likely rather extensive covering significant global areas 
Since they are extensive and could occur anywhere on the globe, a satellite- 
based remote measurement technique is ideally suited for the measurement. Of 
course, if the sink covers the entire globe, the measurement of concentrations 
will not show the sink region. If it is localized to some extent to regions 
which are rather extensive, as would be effected by most sink mechanisms, the 
measurement of CO densities with spatial resolution of the order of 25 to 50 
miles should provide the information needed to locate the sink. 

The carbon monoxide produced by the various sources has time to circulate 
around the globe and at least mix in each hemisphere, although perhaps not 



enough time for appreciable inter-hemispheric transfer. There may or may not 
be enough time for transfer of significant amounts to the stratosphere . Just 
how much mixing is possible depends on the sink and the resultant lifetime which 
is probably somewhere between 0.1 and 1 year, most likely around 0.4 year. 

There are logical indications that there may be some unidentified signifi- 
cant natural source of CO. If this is the case, it suggests that the sink 
strength is much larger than previously supposed. Data are needed to determine 
this. 

' All parts of the CO question require further investigation. Data are 
needed on both localized and global CO measurements, at ground level, in the 
troposphere and in the stratosphere. Data best found by remote and in some 
cases by remote global measurements are needed on other sources of CO, on both 
urban and rural variations, on hemispheric differences, and on CO transport. 

Data are needed related to certain sink processes. Also important are various 
laboratory data such as those required for the analysis of atmospheric chemis- 
try and for those required to quantify biological CO removal processes. The 
sink mechanism problem can be solved only by obtaining sufficient atmospheric 
CO measurement data, supporting atmospheric measurements of related species, 
supporting laboratory data, and ground measurements of the characteristics of 
any sink region found, thus permitting the quantitative determination of the 
sink strength and allowing extrapolation to determine the trend of world-wide 
CO concentration in the future. 


1.2 Measurements Needed 

^I-n-order-to-obtai-n-suffi-ci-ent-data-to-understand" and~eventuaTly“b'e— ab^l'^fd 

solve both short-term and long-term CO problems, measurements of CO are needed 
on a global scale. That is, measurements are needed on a scale which permits 
areas of about the size of cities and of the size of significantly large, rath- 
er uniform areas which may in some way affect the CO balance to be isolated. 
These measurements must also cover the entire globe so as not to overlook the 
possibility of any major sources or sinks. Thus an ideal platform is a space- 
craft with an orbit such that essentially all the globe can be covered with a 
spatial resolution of the order of 25 miles and with a reasonable frequency. 

It is also advisable to be able to investigate certain areas, especially those 
areas found from spacecraft-based measurements to be interesting, by more fre- 
quent measurements with higher spatial resolution, such as those which could be 
obtained from an aircraft. 

The study described herein is on an instrument which seems most ideally 
suited for such measurements - both spacecraft-based and aircraft-based. 
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2 . 


FEASIBILITY OF THE MEASUREfCNT 


2.1 The Feasibility Problem 

An important question in the application of the correlation interferometric 
technique to the measurement of atmospheric burdens of CO concerns the effect of 
variations in atmospheric conditions and profiles. Specifically, variations in 
atmospheric temperature and its profile including inversion layers, atmospheric 
pressure, surface temperature, surface emissivity, surface reflectivity, the at- 
mospheric profile of the gas to be measured, the atmospheric profile of gases 
which interfere with the measurement, cloud cover, the velocity of the platform, 
atmospheric scintillations, and certain other properties were studied to deter- 
mine their effect on the radiation incident on an instrument on a spacecraft. 
This was done by a simulation which consisted of defining a representative mod- 
el atmosphere and subsequently solving the radiative transfer equation to obtain 
the spectrum falling on the instrument. Effects such as reflected sunlight, 
earthshine, atmospheric absorption and atmospheric emission were considered in 
order to yield physically meaningful solutions to the equation. For carbon mon- 
oxide transmission was computed for the spectral regions of the fundamental and 
first overtone bands, i.e., the 4.6 and 2.3 pm regions. The calculations to 
that point were applicable to any instrument. The action of the instrument on 
the incident radiation was then calculated to determine the signal produced. 

For the case in which the instrument is an interferometer, the Fourier transform 
of the spectrum is then calculated to produce an interferogram. Treatment of 
the interferogram produces the measurement of the pollutant. 

The model has been used to investigate the feasibility of the use of the 
correlation interferometer to measure CO, to consider the usefulness of the mea- 
surement in looking for a CO sink (ref. 4), and to define the exact spectral 
band which should be used for most reliable and useful results. The following 
section summarizes the work done on the feasibility study. It is covered in de- 
tail elsewhere (ref. 5) . 


2.2 Atmospheric Models 

2.2.1 Composition Models .- The analytical work performed employed several 
variations of the carbon monoxide, water, and temperature profiles and one car- 
bon dioxide and one methane profile. 

Specifically, ten CO profiles were used. These are shown graphically in 
Figure 2.1. Profile 1 is that for a constant mixing ratio of 0.1 ppm; profile 
2 is that for 10 ppm up to 2 km and 0.1 ppm above that; profile 3 is that for a 
constant mixing ratio of 0.01 ppm; profile 4 is that representing a sink in the 
20-45 km range with the mixing ratio dropping from 0.1 to 0.01 over this alti- 
tude range; profile 5 is that representing a low altitude sink with an effect 
up to 9 km, having a mixing ratio of 0.01 ppm at 0 km and 0.1 ppm at 9 km; pro- 
file 6 represents a low altitude sink with an effect up to 1 km, having a mix- 
ing ratio of 0.01 ppm at 0 km increasing to 0.1 ppm at 1 km; profile 7 also 
represents a low altitude sink with an effect up to 3 km, having a mixing ratio 
of 0.005 ppm at 0 km and 0.1 ppm at 3 km; profile 8 represents a low altitude 
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sink with an effect up to 1 km as in profile 6, but having a mixing ratio of 
0.01 ppm up to 1 km and 0.1 ppm above. Models 9 and 10 representing 0.025 and 
0.25 ppm CO were used for other purposes discussed later. It is important for 
the mapping measurements to show the effects of these differences except for 
profile 4 , the effect of which should be shown in the limb experiment. The to- 
tal number densities from which the CO number densities were obtained with the 
above mixing ratios were taken from Bortner and Kummler (ref. 6) in which they 
were derived mainly from CIRA (ref. 7) representing mean conditions throughout 
the year for latitude near 30°. The CO optical thickness (total number of CO 
molecules per cm^ column of sight) is given for the mapping experiment in Fig- 
ure 2.2 and for the limb experiment . in Figure 2.3. 

The CO2 profile used was based on a constant mixing ratio of 320 ppm. The 
CH4 profile used was based on a constant mixing ratio of 2 ppm (ref. 8) . The 
three water profiles, corresponding to a dry model (1), normal (2), and a wet 
(3) atmosphere, were taken from Gutnick (ref. 9), Oppel (ref. 10), and Linquist 
(ref. 11), respectively, as summarized by Anding (ref. 12). These are given in 
Figure 2.4. These are equivalent to approximately 0.2, 1.5, and 3 percipitable 
cm H2O . 


2.2.2 Temperature and Pressure ftodels .- Four specific temperature pro- 
files were used. These are shown in Figure 2.5. Profile 1 corresponds to a 
cold atmosphere (ref. 13); profile 2 corresponds to an average atmosphere (ref. 
13); profile 3 corresponds to a hot atmosphere (ref. 13); profile 4 involves a 
low altitude inversion layer (up to 2 km), based on average values for Vanden- 
berg AFB, June 1970. It is important that these variations should have no large 
effect on the CO measurements. Any such effects require detailed accurate sup- 
plemental measurements and make the interpretation of the data much more complex 
aiidTTor practical purposes, means t'Kat many fewer ‘da'ta can be^iiit^fpref^C 
Analysis of such effects is covered in subsequent sections of this report. 


2.3 Vehicle Effects 


The planned satellite experiment may be best carried out in an approximate- 
ly polar orbit. Without noting reasons, advantages, and disadvantages, this is 
excellent for the mapping experiment, but restricts the limb measurements to the 
polar regions. This is not a problem, however, since above the tropopause, it 
is to be expected that there is no appreciable latitude effect on the concentra- 
tions. Keneshea (ref. 14) has calculated the OH concentration at different lat- 
itudes and found little variation. Since the chief CO removal mechanism at 
these altitudes is, in all probability, CO + OH reaction, the rate of CO removal 
should be about constant with latitude. The small temperature effect on the 
rate constant should not cause significant variation of CO with latitude. 

One effect which must be considered in remote sensing from a fast moving 
vehicle is the possibility of a Doppler shift. If there is a large velocity 
component in the direction of the observed radiation, there is a Doppler shift 

V 

c 

which is determined by Av = where v^ is the velocity component in the line- 
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Figure 2.2 Optical Thickness as a Function 
of Angle from Vertical 



ALTITUDE (KM) 


Figure 2.3 Optical Thickness for Limb Transmission Experiment 

as a Function of Grazing Altitude 
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Figure 2.4 H 2 O Concentrations 
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o£-sight, or Av = cos 0 where 0 is the angle between the line-of-sight and 

the vehicle movement. In the limb experiment, the satellite may be traveling 
very nearly direct toward or away from the sun so that cos 0^1, and v^ is of 

the order of 7 x 10^ cm sec"^. Thus — = — is of the order of 2.3 x 10^. Thus 

c V 

at 2.3 pm, Av = 0.1 cm“} and at 4.6 pm, Av = 0.05 cm~^. These shifts are of the 
order of a Doppler line-width, which is the width of the line at the upper alti- 
tudes, i.e., those of interest in the limb experiment. If the technique is 
spectral and employs a reference cell of any sort on board this is extremely im- 
portant and probably means that the technique is not suitable for the limb mea- 
surement. For the correlation interferometric technique a Doppler shift essen- 
tially appears as a phase shift which can be readily taken into account. 

As mentioned above, a desirable spatial resolution is of the order of 25 to 
50 miles. From spacecraft altitudes, such as those of the Nimbus vehicle, a 
2° field-of-view gives approximately a 25 mile spatial resolution, a distance 
which is traversed in about four seconds of satellite flight. For the instru- 
ment described herein, a 2° telescope was therefore used in some measurements to 
narrow the normal field-of-view of the instrument from 7°. 


2.4 Spectral Consideration 


2.4.1 CO Spectra .- The problem of measuring a variety of CO burdens or 
of locating and explaining a CO sink depends on having a reliable technique ap- 
plicable to remote sensing. Early in this work, correlation interferometry was 
compared with various other remote techniques considering specificity to CO, 
sensitivity, signal/noise, dynamic range, observation time, spatial resolution, 
stability and data interpretation. Correlation interferometry was found to be 
the most suitable. Such techniques are mainly spectral and involve the measure- 
ment of the amount of absorption or of emission of spectral radiation. Consid- 
eration of the energy levels of CO and their populations shows that the poten- 
tially useable spectral features lie in the ultraviolet, infrared, and micro- 
wave regions. However, in the region of the ultraviolet at which there is sig- 
nificant CO absorption or emission, the atmosphere does not transmit appreciably 
due to the presence of various dissociation and ionization continua. Microwave 
techniques are also not practical for this measurement. The infrared has no 
such continua but has many lines and bands. There are thus a number of windows 
in the infrared, particularly in high resolution spectra of which the most nota- 
ble are the 10-12, 8-9.5, 3.5-4, 2-2.5 and 4.5-5 ym windows (ref. 15). CO bands 
include the following: 


r 

upper 

^lower 

Band 

1 

0 


1 

1 


1 

0 

Fundamental 

2 

1 


2 

0 

1st Overtone 

3 

0 

2nd Overtone 
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The pure rotational lines (V^ = Vj) of CO occur in the microwave region. The 
other bands (V^ > V^) are in the infrared (between about 1 and 5 microns) . 

The fundamental and the first overtone are the bands of most interest here. 
The energies and wavelengths for the lines of these systems are given in Tables 
2.1 and 2.2, respectively. The fundamental band is the strongest. The first 
overtone is about two orders of magnitude weaker and the second overtone is ap- 
proximately that much weaker than the first overtone. The transitions involving 
only excited vibrational states (e.g., 2-1, 3-2, 3-1) are weak because of the 
low populations of the lower state, that of the first excited state being some 
four orders of magnitude lower than the ground vibrational state and the second 
excited state being another four or more orders of magnitude lower. 

The R branch consists of transitions in which the rotational state of the 
upper vibrational level is one greater than that of the lower vibrational level. 
Similarly, the P branch consists of those transitions in which the rotational 
state of the upper vibrational level is one less than that of the lower vibra- 
tional level. That is for v' > v" , the R branch has 3' = 3" + 1, the P branch 

has 3' = 3" - 1. The lines, R or P, are given a subscript denoting the rota- 

tional level of the lower vibrational state. Thus the Rq line of the fundamen- 
tal is between the v = 1, J = 1 and v = 0, J = 0 states. Thus the lines in the 

R branch have a higher energy (lower wavelength) than those in the P branch of 

the same band. 

The strengths of the individual lines in both the fundamental and first 
overtone bands are given in Table 2.3. The integrated band intensities are: 

Fundamental (ref. 16): 260 cm"^ atm"^ 

Overtone (refs. 17,18): 1.64 cm“^ atm"^ 

The advantages and disadvantages of the two bands are discussed later. 

It should be mentioned that the energies and wavelengths given are for the 

normal, most prevalent, isotope of CO, i.e., C^^O^^. Other isotopes have simi- 
lar levels and lines with the lines shifted slightly but significantly. For 
example, the center of the fundamental and first overtone bands of several iso- 
topes of CO are given in Table 2.4, together with the relative abundances (ref. 

19). They play essentially no role in this work. 

2.4.2 Intensity of Reflected Solar Radiation and Earthshine . - 

2.4. 2.1 Solar reflection: The solar radiation reaching the earth and 

being reflected can be calculated for any wavelength as follows. From black- 
body radiation tables (ref. 20) the radiant intensity, J, for a 5900 K black- 
body is calculated to be 9.3 x 10^ watt cm"^ sr“^ ym”^ at 2.3 pm and 8.0 x 10° 
at 4.6 pm. The average irradiance, H, incident on the earth's atmosphere (and 
incident on the earth's surface in the absence of any atmospheric attenuation) 
is 6.8 X 10“^ watts cm"^ at 2.3 pm and 5.5 x 10“** watts cm“^ pm“^ at 4.6 

pm (ref. 21) . 

The reflected solar radiation, N, from the earth's surface is 


12 



TABLE 2.1 ENERGY AND WAVELENGTH OF VIBRATION- 
ROTATION LINES Ol- CO (1-0) 


R BRANCH 


J -J 

cm' 

X ( pm) 

I - 3 

21 47.283 

4. 65748 

2- 1 

2 1 51 . a 53 

4. 64931 

3- 2 

21 54. 598 

4. 641 24 

A- 3 

21 58.331 

4. 63327 

5- A 

21 61 .970 

4. 62541 

6- 5 

21 65. 632 

4. 61 765 

7- 6 

21 69. 1 98 

4. 61003 

8- 7 

21 72. 7 59 

4. 632 4 4 

9- 8 

21 76.282 

4. 59499 

13- 9 

2119. m 

4. 53764 

11-13 

21 83.223 

4. 58339 

12-11 

21 36. 634 

4. 57324 

13-12 

2193.011 

4. 5661 9 

M-13 

21 93. 351 

4. 55923 

15-14 

21 96. 653 

4. 55233 

16-15 

21 99.918 

4. 54562 

17-16 

2233. 1 45 

4. 53397 

18-17 

2236. 334 

4. 53243 

19-18 

2239. 486 

4. 52 59 4 

23-19 

2212. 599 

4. 51957 

21 -23 

221 5. 674 

4. 51333 

22-21 

2218. 713 

4. 53712 

23-22 

2221 . 738 

4. 53134 

24-23 

2224.667 

4. 49 53 5 

25-24 

2227. 537 

4. 4891 6 

26-25 

2233.469 

4. 48336 

.27-26 

2233.310 

4 . 47766 

28-27 

2236. 1 1 3 

4. 472-3 5 

29-28 

2238.875 

4 . 46653 

3'3-29 

2241 . 593 

4. 461 10 

31-33 

2244.281 

4. 45577 

32-31 

2246. 924 

4. 450 53 

33-32 

2249. 527 

4. 44538 

34-33 

2252.389 

4. 44332 

35-34 

2254. 61 1 

4.43536 

36-35 

2257.092 

4. 43048 

37-36 

22 59. 532 

4. 42573 

38-37 

2261 .931 

4. 42100 

39-38 

2264.289 

4. 41 640 

43-39 

2266. 635 

4. 41 1 88 


P BRANCH 


r r t 

J -J 

cm' ‘ 

X(uin) 

0- 1 

2139. 428 

4. 6741 5 

1 - 2 

2135. 548 

4. 68264 

2- 3 

2131. 634 

4. 69 1 24 

3- 4 

2127. 684 

4. 69994 

4- 5 

2123-701 

4. 70876 

5- 6 

2119. 683 

4. 71 769. 

6- 7 

21 I 5. 631 

4. 72672 

7- 8 

2111. 545 

4. 73587 

8- 9 

2107. 426 

4. 7451 3 

9- 1 0 

2103.273 

4. 75449 

13-11 

2099.086 

4. 76398 

11-12 

2094. 867 

4. 77357 

12-13 

2090. 61 4 

4. 78328 

13-14 

2036. 328 

4. 7931 1 

14-15 

2082.010 

4. 80305 

15-16 

2077. 659 

4.81311 

16-17 

2373-275 

4.82329 

17-18 

20 68- 8 60 

4-83358 

18-19 

2364. 412 

4. 84399 

19-2-3 

2359.932 

4. 8 54 53 

20-21 

23 55. 421 

4.8651 8 

21-22 

2353. 873 

4.87596 

22-23 

2346. 303 

4.38686 

23-24 

20 41 . 69 7 

4. 89789 

24-25 

2337.060 

4.93903 

25-26 

2032. 393 

4.92031 

26-27 

2027. 694 

4.931 71 

27-28 

2022.965 

4.94324 

28-29 

2-01 8.20 5 

4.95490 

29-30 

231 3. 41 5 

4.96669 

30-31 

2338. 595 

4.97860 

31-32 

2303. 745 

4.99065 

32-33 

1 998.865 

5.03284 

33-34 

1993.956 

5.01 51 6 

34-35 

1 989.-31 7 

5.02 761 

35-36 

1 984.049 

5.04320 

3 6-37 

1 979.351 

5.35293 

37-38 

1 974.025 

5.06579 

38-39 

1 968.970 

5.07880 

39-40 

1 963-886 

5.09 1 94 



TABLE 2.2 ENERGY AND WAVELENGTH OF VIBRATION- 
ROTATION LINES OF CO (2-0) 


R BRANCH ^ P BRANCH 


J -J 

CM 

- 1 

X(ym) 

f 

J -J 


X(ym) 

1- 0 

42 63. 

831 

2. 34531 

0- 1 

4256.21 1 

2.34951 

2- 1 

4267. 

536 

2.34327 

1- 2 

42 52. 29 6 

2.351 67 

3- 2 

4271 . 

1 73 

2.34128 

2- 3 

4248.31 1 

2.35388 

4- 3 

4274. 

734 

2.33933 

3- 4 

4244.257 

2.35612 

5- A 

4278. 

227 

2.33742 

4- 5 

4240. 133 

2.35842 

6- 5 

4281 . 

6 49 

2.33555 

5- 6 

4235. 949 

2.36075 

7- 6 

4285. 

000 

2. 33372 

6- 7 

4231 . 678 

2.36313 

a- 7 

4288. 

280 

2.33194 

7- 8 

4227.347 

2.36555 

9- 6 

4291 . 

488 

2.33019 

8- 9 

4222.947 

2.36801 

10- 9 

4294. 

625 

2.32849 

9-10 

4218.479 

2.370 52 

11-10 

4297. 

690 

2.32683 

10-11 

421 3. 942 

2.37307 

12-11 

4300. 

684 

2.32 521 

11-12 

4209.337 

2.37567 

13-13 

4303. 

60 5 

2. 32363 

12-13 

4204* 664 

2.37831 

14-13 

430 6. 

4 54 

2. 32210 

13-14 

4199.923 

2.38100 

15-14 

4309. 

231 

2.32063 

14-15 

4195. 1 1 4 

2.38373 

16-15 

4311. 

935 

2.3191 4 

15-16 

4190.237 

2.38650 

17-16 

4314. 

567 

2.31 773 

16-17 

41 85.293 

2.38932 

18-17 

431 7. 

125 

2. 31 636 

17-18 

4180.282 

2.39218 

19-18 

4319. 

611 

2.31 502 

18-19 

41 75.203 

2.39509 

23-19 

4322. 

324 

2.31373 

19-20 

41 70.058 

2.39805 

21-23 

4324. 

3 63 

2.31248 

20-21 

41 64. 845 

2 . 40 1 0 5 

22-21 

4326. 

628 

2. 31 127 

21-22 

41 59. 567 

2.40410 

23-22 

4328. 

821 

2.31019 

22-23 

41 54.221 

2. 4071 9 

24-33 

4333. 

9 39 

2. 30697 

23-24 

41 48* 810 

2.41033 

25-24 

4332. 

983 

2. 30788 

24-25 

4143. 332 

2. 41 352 

26-25 

4334. 

9 53 

2. 30 683 

2 5-26 

41 37. 788 

2.41675 

27-26 

4336. 

8 49 

2.33 582 

26-27 

41 32. 1 79 

2.42303 

28-27 

4338. 

670 

2. 39485 

27-28 

4126. 504 

2.42336 

29-28 

4340. 

417 

2.30393 

28-29 

4120. 763 

2.42674 

33-39 

4342. 

0 89 

2.30 334 

29-30 

41 1 4.957 

2.4301 6 

31-39 

4343. 

686 

2.33219 

33-31 

4109.086 

2.43363 

32-31 

4345. 

208 

2. 33139 

31-32 

4103. 1 50 

2.4371 5 

33-33 

4346. 

654 

2.30362 

32- 33 

4097. 1 49 

2. 440 72 

34-33 

4348. 

026 

2.29989 

33- 34 

4091 .083 

2. 44434 

3 5-34 

4349. 

321 

2.29921 

34-35 

4384. 953 

2. 4480 1 

36-35 

43 50. 

541 

2.29856 

3 5-36 

4073. 759 

2. 451 73 

37-36 

4351 . 

665 

2.29796 

3 6-37 

4072. 501 

2. 45549 

38-37 

4352. 

7 53 

2.29740 

37-38 

4066. 1 78 

2.45931 

39-38 

43 53. 

745 

2.29687 

38-39 

4059. 792 

2.46318 

40-39 

4354. 

669 

2.29639 

39-40 

4053. 342 

2.46713 



TABLE 2.3 CO LINE STRENGTHS (cm-^ atm"^) 


FUNDAMENTAL OVERTONE 



( ref . 

19 ) 

( ref . 

22 ) 


R 

P 

R 

P 

0 

1.970 


. 0163 


1 

3 . 878 

1 . 927 

. 0321 

.0158 

2 

5 . 643 

3 . 707 

. 0467 

. 0303 

3 

7.100 

5.251 

. 0592 

.0428 

4 

8.261 

6.491 

. 0691 

. 0527 

5 

9.059 

7 . 386 

.0760 

. 0598 

6 

9.482 

7.920 

. 0797 

. 0638 

7 

9 . 544 

8 . 105 

.0805 

.0651 

8 

9:286 

7.980 

. 0785 , 

.0638 

9 

8 . 760 

7.591 

.0742 

.0604 

10 

8 . 032 

7 . 002 

. 0682 

.0555 

1 1 

7 . 170 

6 . 278 

.0610 

.0495 

12 

6 . 240 

5 . 480 

. 0532 

. 0430 

13 

5 . 301 

4 . 665 

.0453 

. 0365 

14 

4 . 400 

3 . 876 

. 0376 

,0301 

15 

3 . 568 

3 . 145 

.0306 

.0243 

16 

2.830 

2.494 

.0243 

.0192 

17 

2 . 196 

1. 934 

. 0189 

.0148 

18 

1. 668 

1.486 

.0144 

.0112 

19 

1. 241 

1. 091 

.0107 

.00828 

20 

0 . 904 

0.793 

.00779 

. 00597 

21 

0 . 645 

0 . 566 

. 00557 

. 00423 

22 

0.451 

0 . 395 

. 00390 

. 00294 

23 

0 . 309 

0 . 290 

. 00268 

. 00200 

24 

0 . 208 

0 . 181 

. 00180 

. 00133 

25 

0 . 137 

0 . 119 

. 00119 

. 000872 

26 

0.0884 

0.0766 

. 000765 

.000556 

27 

0. 0589 

0. 0483 

. 000486 

. 000350 

28 

0. 0349 

0. 0300 

. 000302 

.000215 

29 

0. 0212 

0. 0182 

. 000184 

. 000130 

30 

0. 0127 

0. 0108 

. 000110 

.0000770 

31 

0 . 00743 

0 . 00634 

. 0000642 

. 0000446 

32 

0. 00427 

0 . 00364 

. 0000371 

, 0000254 

33 

0 . 00241 

0 . 00205 

. 0000209 

. 0000142 

34 

0 . 00134 

0 . 00113 

, 0000116 

. 0000078 



TABLE 2.4 CO BAND CENTERS FOR ISOTOPES 


V V 

o 0 


Isotope 

(Fundamental) 

(1st Overtone) 

Abundance 

C12Q16 

2143.273 

4260.056 

.98654 

C12Q17 

2116.111 

4206.675 

.00037 

C12Q18 

2092.124 

4161.021 

.00202 

Cl3ol6 

2116.075 

4206.825 

.01106 

Cl30l7 

2068.460 

4112.264 

.0000041 

Cl3ol8 

2043.701 

4062.639 

.0000226 



N(watts cm~^ sr-^ Mm 


-1) = pH 


The values of N are thus 2.16 x 10“^ p watt cm~^ sr-^ pm-^ at 2.3 pm and 
1.68 X 10-^* p watt cm-^ sr-^ pm-^ at 4.6 pm. The intensities (in watt cm-^ 
sr-^ pm-^) are shown in Figure 2.6 for various reflectivities. 

2. 4. 2. 2 Earthshine: From blackbody radiation tables (ref. 20) the radiant 

intensity is given for a 300 K blackbody by the ratio 



1.439 X 10*t 
e 

X X 300 ^ 

^300 

1.439 X lO** 
e 


the intensities (watts cm”^ sr"^ pm“^) 
culated giving the following, as shown 

for 273.2 and 288 K blackbodies were cal- 
in Figure 2.6 

T 

_g 

2.3 pm 

4.6 pm 

273.2 

2.0 X 10-® 

6.15 X 10-5 

288 

7.54 X 10-® 

1.11 X 10-‘t 

300 

1.60 X 10-'=' 

1.71 X lO-** 

2. 4. 2. 3 Ratio of reflected solar radiation and earthshine intensities: As 

seen in Figure 2.6 the radiation seen in this experiment in the overtone band 
region would consist almost exclusively of reflected sunlight while that radia- 
tion for the fundamental band region would be mainly earthshine with some effect 


of reflected sunlight and serious atmospheric emission effects. These factors 
will be discussed later. 

From the above numbers the ratios of reflected sunlight to earthshine at 
2.3 pm are found to be 


p_ 

T : 

273.2 K 

288 K 

300K 

0.02 

g 

2 X 1Q3 

5.4 X 10^ 

2.5 X 102 

0.1 


1 X lO** 

2.7 X 103 

1.2 X 103 

1.0 


1 X 10^ 

2.7 X lO** 

1.2 X 10** 

earthshine 

to reflected 

sunlight at 46 

pm are : 

P_ 

T : 

272. 2K 

288 K 

300 K 

0.02 

g 

1.78 X 10^ 

3.20 X lOl 

4.95 X lOl 

0.1 


3.56 X 10° 

6.40 X 10° 

9.90 X 10° 

1.0 


3.56 X 10-1 

6.40 X 10-1 

9.90 X 10- 
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2.5 Radiative Theory 


To study the feasibility of the CO density measurements and the effect of 
various atmospheric parameters on such measurements, calculations have been made 
using the radiative transfer theory described in this section. The model used 
includes absorption and emission of photons along a ray from the source to the 
detector. The two geometries considered are direct solar observation - limb 
experiment, and absorption of the earth's albedo - mapping experiment. Provi- 
sions are made for adding earthshine to reflection of sunlight from the earth's 
surface. However, phenomena such as scattering resonance fluorescence and radi- 
ation trapping have been neglected. Boltzmann populations, as described by the 
local atmospheric temperature, have been assumed for the upper and lower states 
of all transitions and the line profiles are given by the Voigt function at all 
altitudes. Doppler shifts as caused by atmospheric winds and by relative motion 
between the atmosphere and the satellite have been neglected. 

The monochromatic radiative transfer equation (ref. 23) is given by: 





dl 

V 

dS 

e - I a 

V V V 

Where : 

I 

V 

= spectral 

intensity 

(watts/cm^ -sr-cm" ^ ) 


S 

= distance 

along ray 

(cm) 


e 

V 

= spectral emission coefficient (watts/cm^-sr-cm"^) , 
as a function of S 


a 

V 

= spectral 
tion of S 

absorption 

coefficient (cm"^), as a func- 


Formally, 

= C exp(-/a^dS") + exp(-/a^dS^) J exp(/a^dS^'') dS^ (2.2) 
For the limb boundary condition: 

z „ ^ 

I (z) = I exp(-/^a dS"") + exp(-/^a dS"") f e exp(/ a dS'") dS" (2.3) 

V ^ o,v ^ o^v ^ ov J^v^ ov 

or 

z ^ z 

I (z) = I exp(-/ a dS"") + f e exp(-/^,a dS"") dS" (2.4) 

o 

The geometry for this situation is shown in Figure 2.7. 

For the mapping experiment, two conditions are required. The solar energy 
is attenuated as it penetrates to the earth's surface. At this point, it is 
reflected and the total intensity increased by earthshine. The resulting flux 
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Figure 2.7 Geometry for Limb Conditions 





is then further attenuated by the atmosphere until the satellite is reached. 

It is pointed out that atmospheric emission which is reflected by the earth in- 
to the field-of-view of the detector is negligible and therefore not included 
in the model. The geometry of the mapping experiment is given in Figure 2.8 
and is described mathematically by Equation 2.5. 

Ix = Pv f^o,v expC-/^ladS'0 + expC/gla^dS dS^J (2.5) 

+ Ig expC-/^a^dx'") + exp(-/.a^dx'0 dx" 

Where 

= earthshine (watts/cm^-sr-cm-^) 
p = earth reflection coefficient 

V 


Equations 2.4 and 2.5 have been numerically integrated along the rays indicated 
in Figures 2.7 and 2.8 using Simpson's rule. The properties e and a are func- 
tions of altitude which must be evaluated prior to integration^ This^is con- 
sidered below. 


On a microscopic basis, for a single rotational line (ref. 24), 
n hcv 

e = £(p = — -j A 

V V 4tt uJl V 


( 2 . 6 ) 


01 


V 


f“A ■ “S> 





n 

u 




hcv , 
4tt 


(2.7) 


Where 


u’ Z 
hcv 


uZ 


Zu 



Voigt profile (1/cm"^) 

upper and lower state number densities (particles/cm^) 

energy per transition (watt sec/transtion) 

absorption coefficient (cm“Vcm) 

stimulation emission coefficient (cm“Vcm) 

Einstein coefficient for spontaneous emission 

/ Transitions 
I Sec. Particle 

Einstein coefficient for absorption 

( Transitions Cm-Ster 
Particle ’ Watt Sec. 

Einstein coefficient for stimulated emission 

/ Transitions Cm-Ster 
^Particle ‘ Watt Sec. 



21 








It is noted that the definition of the Einstein coefficients is arbitrary but 
that they are related by Kirchoff's Law for equilibrium at temperature, T. 


Thus , 


or 


e 

V 


B CT) a 

V V 


B (T) 

V 


2 hc^ 

hcv/kT 

e 


(watits/cm^-sr-cm 


n A = (n„ - n B J B (T) 

u ui, Jl Jlu U ui V 


A ./B „ 

u!L u£ 

(n. B„ /n B ) 
^ I 2.U u u£, 


2 hc^ 

exp(hcv/kT) -1 


But, at equilibrium the Boltzmann distribution is: 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


— = — exp(hcv/kT) 

u ®u 


( 2 . 11 ) 


and by comparison 




g B „ 
*u u£ 


( 2 . 12 ) 


A „ = B • 2hc2 v3 (2.13) 

u£ uX, ^ ^ 

The Voigt profile, (|)y, is defined in terms of the Doppler and Lorentz half 
widths, i.e., half the line width at half the peak intensity (ref. 25). 




1/2 


K(x,y) 


K(x,y) 


= L f exp(-t^) 

+ (x - t)2 


(2.14) 

(2.15) 


a 

y = — (ln2)l/2 C2.16) 

“d 

(v - V ) 

X = 2_(in2)l/2 (2.17) 

“d 
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(2.18) 


a„ = Doppler half width Ccm"^) 


\ | 2R(ln2) 
c I M 

Lorentz half width (cm"^) 


= 3.58 X 10"’^ V 


Ti 

Mi 


(2.19) 


= r/27TC 


1 .4 p N 


' 87T U/2 


12 0 \M R T 


c 


Where 

r = gas collision frequency (sec”^) 

p = pressure (dynes/cm^) 

= optical collision cross-section fcm^) 

12 

M = reduced mass of the colliding species (gm/gm mol) 

R = gas constant (erg/gm mol K) 

T = temperature (K) 

Nq = Avagadro's number (particles/gm mol) 

c = velocity of light (cm/sec) 

The Doppler and Lorentz widths are functions of temperature and pressure and 
thus are easily obtained at any point along the ray from the atmospheric prop- 
erties. The integral, K(x,y), is then calculated using a subroutine developed 
by Armstrong (ref. 25). It should be noted that 


The emission coefficient number density product necessary for Equations 2.6 and 
2.7 is most easily obtained by scaling experimentally derived line strengths 
which are defined as (ref. 26) . 


Where 


5.U 


1 ® 
^ n — 

8TTCV2 P g 




exp 


hcv 

kf 



( 2 . 20 ) 
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V 


= line strength (cm“^-atm“^) 
= wavenumber (cm"^) 



P = pressure (atm) 

= lower state population for conditions where S 
is measured (particles/cm^) 

g ,g. = statistical weights (2J + 1) 

U 

A „ = Einstein coefficient 

u£ 

Transitions 
Sec. Particle 

Using equation 2.11 and rearranging 


n A „ 
u u£ 


Sitcv^P 

[exp (hcv/kf)-l] 


Jlu 


( 2 . 21 ) 


Where 


But 


n = upper state population for conditions where 
is measured (particle/cm^) . 


\ \ Q exp (-hcv/kT) 

^u ”t ^ ®^P C-hcv/kf) 


Where 


n^ = total number density 
Q = partition function 
E = upper state energy 


Thus, combining Equations 2.19 and 2.20 


n A . 
u u£ 


Sircv^PS . 

£u 


exp^- 

hcv ri 
k [t 

-i] 


expj^ 

■ j hcv \ 

'Ikf 1 


Qi>t 


(2.23) 


If several lines contribute to the absorption at a specific wavenumber the 
total absorption coefficient is obtained by summation. The computational algo- 
rithm then uses as a basis, n A „ from Equation 2.21 gives (n. B. - n 
£u u uil ^ ^ il ilu u 

B^j^) from Equation 2.9. This, in conjunction with (j>^, leads to ot^ via Equation 
2.7. Equation 2.5 can then be integrated over altitude. 
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2.6 Feasibility Analysis 


2.6.1 Programs . - Calculations have been made with several programs. 
Initial calculations were made with a single-line model. For each of several 
lines of the overtone and fundamental bands, calculations were made of absorp- 
tion across the line, the line shape as a function of altitude, net absorption 
as a function of altitude for various CO density profiles (including sinks and 
various temperature and emissivity, and of other factors which influence the 
net intensity). With a second program, calculations were made, for the limb 
mode, of the effect of instrument error on the inversion of measured total CO 
densities in the path to obtain a CO profile. With a multi-line program which 
computes the spectrum incident on the instrument, generates the corresponding 
interfere gram, and, with another program, inverts this to CO densities in the 
path, calculations were made to determine the sensitivity of the technique and 
the effects of various atmospheric parameters. 

Parameters which can be varied by means of input data include temperature 
and species concentration profiles of the atmosphere, ground temperature and 
emissivity, reflectivity, viewing angles, incident radiation angle, and grazing 
height. 

The net change in intensity for the R7 line of the fundamental is given in 
Figure 2.9. This shows by the solid curve the results obtained by the single- 
line model and by the x's the results obtained by the multi-line model. These 
agree within 1% and it can be concluded that the comparison of the fundamental 
and overtone bands can be carried out using the single-line model. 

For all calculations discussed in the following pages, the following condi- 
tions were used, as standard and variations made therefrom: 

CO Profile 1 
Temperature Model 2 

Ground Temperature = 0 km Atmospheric Temperature 
Emissivity =1.0 
Reflectivity = 0.1 
Bandpass = 20 cm-^ 

Lorentz Half-Width (R7) = 0.06 cm"^ 

Calculations have been made for the CO profile models 1 through 6 described in 
Section 2.2.1 (Figure 2.1). These models were used to evaluate the emission 
and absorption coefficients and compute the CO absorption spectrum which would 
be observed at an altitude of 80 km. 

2.6.2 Single-Line Model . - 

2. 6. 2.1 Line shapes: The absorption and emission coefficients are func- 

tions of pressure, temperature, and CO number density. This is illustrated by 
Figure 2.10 which shows the spectral absorption coefficient of the 4.6 ym P8 
line as a function of altitude for the standard temperature and CO models. At 
low altitudes the lines are Lorentz broadened and then as altitude increases, 
Doppler broadening dominates. The variation in the magnitude of the absorption 
coefficient with altitude is controlled by the variation of the CO number den- 
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sity. IVhen these absorption coefficients are integrated over the atmosphere, 
the line profile variation with altitude is manifested by an emission peak at 
the center of the absorption line. This behavior is illustrated in Figure 2.11. 
Physically, earthshine is absorbed up to an altitude of about 30 km. At this 
point, the atmospheric temperature rises and the light flux is increased by emis- 
sion. However, because the line profiles are very sharp at this altitude, the 
emission is observed only near the line center. 

A second example of the importance of atmospheric emission is illustrated 
in Figure 2.12. Here the net change in absorbed intensity with altitude as inte- 
grated over the P8 line of the fundamental is shown as a function of altitude for 
the four temperature models and the standard CO model. Two temperature profiles 
have inversion layers which cause a net increase in the light flux (negative ab- 
sorption) . When the temperature drops there is less photon emission and the net 
absorption increases. For the case of the P8 line of the overtone this behavior 
is not observed, as shown in Figure 2.13. The reason for this is source temper- 
ature. In the case of the overtone, reflected sunshine dominates the intensity 
at the earth's surface, whereas at 4.6 ym earthshine is the dominant contributor 
to the ground intensity. The atmospheric temperature can be higher than the 
ground temperature and hence atmospheric emission becomes important. 

2.6. 2. 2 CO profile effects: The effect of CO number density profile on 

the intensity has been examined using the standard temperature profile and the 
six CO models. The results for the change of intensity with altitude for the P8 
line of the fundamental are shown as a function of altitude in Figure 2.14. In 

all cases, except the urban atmosphere, the change in intensity peaks between 6 

to 10 km whereas it should peak at ground level except for the low-altitude sink 
models. This peak is caused by a balance between the changing CO concentration 
and the emission. In the case of the urban atmosphere, the CO concentration pro- 
file dominates. This balance should be emphasized in that the change in absorp- 
tion at 4.6 ym is not directly proportional to the CO concentration. On the oth- 

er hand, similar calculations for the P8 line at 2.3 ym are a direct measure of 
the CO number density. This is shown in Figure 2.15. The difference between 
the two wavelengths is again caused by the effective source temperature. 

The data shows that absorption spectra as observed from 80 km depend on the 
atmospheric profiles of the absorbing species and on temperature. The tempera- 
ture dependence on the spectra is controlled by the temperature of the absorbed 
source. For the case where earthshine is important the absorption profiles are 
strongly effected by the atmospheric temperature profile. Also in this case, the 
change in absorption with altitude is not a direct measure of the absorbing spe- 
cies. On the other hand, when the absorbed source is the sun, the absorbed sig- 
nal is almost independent of the temperature profile and the change in absorption 
with respect to altitude is directly proportional to the concentration of the 
absorbing species. 

The results are shown in Table 2.5 for the six models and in Table 2.6 for 
models in which concentrations of the CO were multiples (.1, .2, .4, .8, 1., 1.6, 
3.2, and 6.4) of the standard atmosphere at all altitudes. 

The effect of low altitude sinks on measured CO densities is seen in Tables 
2.7 and 2.8 using temperature models 2 and 4, respectively. The numbers pre- 
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ALTITUDE (km) 


Figure 2,12 Variation of Rate Change in Absorbed 

Intensity Fundamental (4.6 urn) (R7 Line) 
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Figure 2.13 Variation of Rate of Change in Absorbed 
Intensity Overtone (2.3 ym) (R7 Line) 
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Figure 2.14 Variation of Rate Change in Absorbed 
Intensity Fundamental (4.6 pm) 
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Figure 2.15 Variation of Rate of Change in Absorbed 

Intensity Overtone C2.3 um) 


34 


TABLE 2.5 EFFECT OF CO PROFILES ON ABSORPTION 


FUNDAMENTAL (4.6 ym) 


CO MODEL 


R7 

PI 

P8 

P26 

STANDARD 

2. 15 (18) 

. 2002 

. 0810 

. 1700 

. 00104 

LOW ALTITUDE SOURCE 

5.21 (19) 

.6620 

. 2974 

. 5695 

. 0167 

LOW CO 

2. 15 (17) 

. 0439 

. 0119 

. 0350 

. 000107 

HIGH ALTITUDE SINK 

2. 09 (18) 

. 2006 

. 0807 

. 1703 

. 00102 

LOW ALTITUDE SINK 

1. 37 (18) 

. 1580 

. 0642 

. 1335 

. 000614 

(9 kjTi) 

LOW ALTITUDE SINK 

1.93 (18) 

. 1963 . 

.0795 

. 1 666 

. 000965 

(3 km) 


OVERTONE (2. 

3 ym) 



CO MODEL 

nCo(cm‘^ 

R7 

PI 

P8 

P26 

STANDARD 

2. 15 (18) 

. 00647 

. 00144 

. 00462 

. 0000098 

LOW ALTITUDE SOURCE 

5.21 (19) 

. 1226 

.0299 

.0953 

. 000431 

LOW CO 

2. 15 (17) 

.000662 .000146 

.000471 

. 000001 1 

HIGH ALTITUDE SINK 

2.09 (18) 

. 00631 

.00140 

. 00451 

. 0000097 

LOW ALTITUDE SINK 

1. 37 (18) 

.00418 

. 00095 

. 00292 

. 0000037 

(9 km) 

LOW ALTITUDE SINK ' 

1. 93 (18) 

. 00582 

. 00130 

.00 413 

. 0000068 

(3 km) 



TABLE 2.6 EFFECT OF CO COLUMN DENSITY ON FRACTIONAL NET ABSORPTION 



FUNDAMENTAL i 

(4.6 ym) 



R7 

PI 

P8 

P26 

Z. 15 (17) 

. 0439 

. 01 19 

. 0350 

. 000107 

4.3 (17) 

. 0741 

. 0223 

. 0604 

. 00020 

8.6 (17) 

. 1174 

. 0405 

.0978 

. 000410 

1.72 (18) 

. 1769 

. 0692 

. 1496 

. 000820 

2. 15 (18) 

. 2002 

.0810 

. 1710 

.00104 

3. 44 (18) 

. 2585 

. 1108 

.2205 

. 00163 

6. 88 (18) 

. 3714 

. 1683 

. 3183 

. 00322 

1. 37 (19) 

. 5292 

. 2470 

. 4545 

. 00629 


OVERTONE (2 

. 3 ym) 


nco(cm'^ 

R7 

PI 

P8 

P26 

2. 15 (17) 

. 000662 

. 000146 

. 000471 

. 0000011 

4.3 (17) 

.00131 

. 000286 

. 000926 

. 0000015 

8.6 (17) 

. 00260 

.000569 

.00185 

. 0000032 

1. 72 (18) 

.00518 

. 00115 

. 00369 

. 0000071 

2.15 (18) 

. 00647 

. 00144 

. 00462 

. 0000099 

3.44 (18) 

.0102 

. 0C228 

. 00729 

. 000014 

6. 88 (18) 

.0197 

. 00455 

. 0142 

. 000030 

1. 37 (19) 

. 0370 

. 00895 

. 0272 

. 000059 



TABLE 2.7 EFFECT OF CO PROFILES ON DERIVED CO 
CONCENTRATIONS (% DEVIATION) 


FUNDAMENTAL (4.6 pm) 


CO MODEL 


-R7 

PI 

P8 

P26 

HIGH ALTITUDE SINK 

2.09 (18) 

0.2 

0. 0 

8. 3 

2. 5 

LOW ALTITUDE SINK (9 km) 

1.37 (18) 

14.6 

10. 5 

18.0 

4.9 

LOW ALTITUDE SINK (3 km) 

1.93 (18) 

8. 1 

2.7 

12. 7 

1.0 


OVERTONE 

(2.3 pm) 





, -2 





CO MODEL 

"co^^"' > 

R7 

PI 

P8 


HIGH ALTITUDE SINK 

2.09 (18) 

3.5 

4. 5 

3. 5 


LOW ALTITUDE SINK (9 km) 

1. 37 (18) 

6. 1 

6. 1 

6. 1 


LOW ALTITUDE SINK (3 km) 

1.93 (18) 

2.7 

2. 7 

2.7 



TABLE 2.8 EFFECT OF PROFILES ON DERIVED CO CONCENTRATIONS - 
TEMPERATURE INVERSION (% DEVIATION) 

FUNDAMENTAL (4.6 pm) 


CO MODEL 

nco<cm ) 

R7 

PI 

P8 

P26 

LOW ALTITUDE SINK (9 km) 

1.37 (18) 

13. 3 

10.4 

17.0 

0. 04 

LOW ALTITUDE SINK ( 3 km) 

1.93 (18) 

0.9 

29.5 

32.6 

35.2 


OVERTONE 

(2.3 pm) 




CO MODEL 


R7 

PI 

P8 


LOW ALTITUDE SINK (9 km) 

1. 37 (18) 

7.4 

7.4 

7.4 


LOW ALTITUDE SINK (3 km) 

1.93 (18) 

3. 6 

4. 7 

2.6 




sented are those of the precent difference of the actual CO density in the model 
and that calculated using the computed net absorption for that model and con- 
verting that by means of the data of Table 2.6. It can be seen that at 4.6 ym 
the errors are larger than the differences in the amounts of CO in the standard 
and the sink models. To improve these results, an atmospheric CO profile would 
have to be known for measured CO densities; this cannot be assumed to be obtain- 
able. The errors for the temperature model 2 are less drastic than those of 
model 4 which has an inversion layer. The latter errors are extremely large for 
the 4.6 ym band. 

2. 6. 2. 3 Temperature profile effects: Calculations have been made for tem- 

perature models 1, 2, 3, and 4 (Table 2.9, Figure 2.12) using the 0.1 ppm CO 
model. The calculations shown in Figures 2.12 and 2.13 were made for the R7, PI, 
P8, and P26 lines of the fundamental band and the R7, PI, P8 and P26 lines of the 
overtone band, using CO model 1. The results for this are shown in Table 2.9. 

The tabulated results show a great dependence on the atmospheric temperature mod- 
el for the 4.6 ym band and little dependency on the 2.3 ym band. The graphical 
data of Figures 2.12 and 2.13 show the contribution of various altitudes to the 
net absorption. The 2.3 ym band shows absorption as a function of altitude to 

be as would be expected whereas the 4.6 ym band shows much less absorption below 
about 5 km (and, in fact, a significant net emission in some regions) which is 
an important region for CO sink studies. 

2.6. 2.4 Ground temperature effects: The calculations described above in- 

clude a ground temperature of each model at 0 km. Calculations have also been 
made for a different ground temperature. The following cases were carried out. 


T 

Ground 

Atmospheric 

Temperature 

Model 

^Model , 
at 0 ! 

288.2 

2 

288.2 

288.2 

1 

257.3 

288.2 

3 

302.6 

288.2 

4 

283.2 

273.2 

2 

288.2 

284.2 

2 

288.2 

287.2 

2 

288.2 

289.2 

2 

288.2 

291.2 

2 

288.2 

308.2 

2 

288.2 


The results of these calculations are shown in Table 2.10. Drastic effects are 
seen for the fundamental band (4.6 ym) and negligible effects are seen for 2.3 
ym. 


2. 6. 2. 5 Emissivity effects: In calculations described above, an emissiv- 

ity of 1 was used. Since certain ground areas have appreciably different emis- 
sivities, calculations have been made using emissivity values of 0.95, 0.9, 0.8, 
and 0.7 for the 4.6 ym band and a value of 0.7 for the 2.3 ym band, with each of 
the four temperature models. Drastic effects in the fractional net apparent ab- 
sorption are seen by the data in Table 2.11 at 4.6 ym while those at 2.3 ym are 
negligible. 
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TABLE 2.9 EFFECT OF TEMPERATURE PROFILE ON ABSORPTION 


4.6 ym 

R7 

PI 

P8 

P26 

LOW TEMPERATURE 

. 1613 

.0677 

. 1328 

. 000489 

NORMAL TEMPERATURE 

. 2002 

.0810 

. 1700 

. 00104 

HIGH TEMPERATURE 

. 1942 

. 0772 

. 1662 

.00134 

INVERSION LAYER 

. 1270 

. 0521 

. 1063 

. 000433 


2.3 ym 


LOW TEMPERATURE 

. 00660 

.00152 

.00462 

.0000062 

NORMAL TEMPERATURE 

.00647 

. 00144 

. 00462 

. 0000098 

high temperature 

.00637 

. 00140 

. 00461 

. 000014 

INVERSION LAYER 

.00639 

. 00140 

. 00461 

. 000013 



TABLE 2.10 EFFECT OF GROUND TEMPERATURE ON ABSORPTION 


I 


FUNDAMENTAL (4.6 ym) 




R7 

PI 

P8 


257.3 

1613 

0677 

. 1328 

LOW TEMPERATURE 

288.2 

2500 

1030 

. 2117 


302.6 . 

1942 

0772 

. 1662 

HIGH TEMPERATURE 

1 

288.2 

1342 

0550 

. 1128 


283.2 

1270 

0521 

. 1063 

INVERSION LAYER 

288.2 

1540 

0622 

. 1302 


OVERTONE 

(2.3 ym) 





R7 

PI 

P8 


257.3 

00660 

00152 

. 00462 

LOW TEMPERATURE 

288.2 

00657 

00150 

.00459 


302.6 

00637 

00140 

. 00461 

HIGH TEMPERATURE 

288.2 

00634 

00138 

. 00457 


283.2 

00639 

00140 

. 00461 

INVERSION l^YER 

288.2 

00635 

00139 

. 00458 


P26 

.000489 
. 00104 

. 00134 
. 00415 

. 000433 
. 000793 

P26 

. O0OOO61 
. 0000053 

. 0000137 
. 0000122 

. 0000128 


. 0000144 



TABLE 2,11 EFFECTS OF EMISSIVITY ON FRACTIONAL NET ABSORPTION 

FOR R7 LINES 


TEMPERATURE MODEL 


LOW 

TEMPERATURE 


NORMAL 

TEMPERATURE 


NORMAL 
TEMPERATURE 
Tg = 288.2 


HIGH 

TEMPERATURE 


INVERSION 

LAYER 


€ 

4.6 pill 

2.3 pm 

1 

. 161 

. 00660 

0.95 

. 155 


0.9 

. 149 


0.8 

. 133 


0.7 

. 113 

. 00657 


1 

.200 

.00647 

0.95 

. 196 


0.9 

. I9I 


0.8 

. 180 


0.7 

. 165 

.00644 


1 

.200 

0.9 

. 191 

0.7 

. 165 


1 

, 194 

.00637 

0.95 

. 190 


,0.9 

. 184 


0. 8 

. 172 


0.7 

. 156 

. 00634 


1 

. 127 

.00639 

0.95 

.119 


0.9 

. no 


0.6 

.0886 


0,7 

.0710 

, 00635 



It can be seen that atmosphere and ground temperatures and ground emissivi- 
ties would have to be very accurately measured to be able to interpret any CO 
data obtained by remote measurements at 4.6 pm. It appears that the best expect- 
ed temperature data will have an accuracy of ± 2 K (refs. 27, 28, 29). These 
data would not be expected to be this accurate if an inversion layer exists. 

Even ± 2 K would present substantial effects at 4.6 pm, since the emissivity may 
vary with wavelength. 

2. 6. 2. 6 Effects of other parameters; Calculations made to test the effect 
of reflectivity show this effect to be small for the 4.6 pm band and negligible 
for the 2.3 p band. The effect is due to a change in source intensity with no 
accompanying change in atmospheric emission. 

Calculations were made to show the effect of Lorentz half-width on absorp- 
tion. Significant effects were found for some conditions for the 4.6 pm band 
while the effects on the 2.3 pm band were small. Results are shown in Tables 
2.12 and 2.13. This shows that the value of Equation 2.21 must be known 

to accuracies of the order of 10% for the fundamental, while for the overtone 
an accuracy of 50% is sufficient. 

Calculations to test the effect of bandpass (the distance from the center 
of the line over which the absorption is integrated) were made. As expected 
this has an appreciable effect if too small a bandpass is taken. Data for four 
lines of the 4.6 and 2.3 pm bands are given in Tables 2.14 and 2.15 for band- 
passes (v = Vq) of 20, 2, and 0.1 cm"^. 

2. 6. 2. 7 Comparison of overtone and fundamental bands: Using the data giv- 

en in the preceding tables, a summary may be made of the two general wavelength 
regions, the 2.3 pm band of the first overtone of CO and the 4.6 pm band of the 
fundamental of CO, 

The absorption of the 2.3 pm radiation has been found to be sufficient to 
produce the required sensitivities and ranges using the correlation interferom- 
eter technique for optical thicknesses appropriate to both the mapping and the 
limb experiments, considering nominal ambient atmospheric concentrations and 
reasonable fractions thereof. 

The absorption of the 4.6 pm radiation is, of course, sufficient for the 
experiments. One difficulty arises at the higher optical thicknesses appropri- 
ate to the lower altitude limb measurements. Under such conditions the absorp- 
tion is so strong that it is well off the linear portion of the curve -of- growth 
with ambient atmospheric CO concentrations resulting in a loss of sensitivity 
in the CO measurements. 

The data show that there is considerably more error in measurement at 4.6 
pm of the CO column density in the important case of a low altitude sink than 
for the same measurement at 2.3 pm. Since a major objective of the experiment 
is to search for a potential low altitude sink, this indicates a definite advan- 
tage for 2.3 pm. In the case of a temperature inversion drastic errors are 
found in the case of a low altitude sink. Since temperature inversions are com- 
mon the 2.3 pm has another distinct advantage. In general it can be said that 
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TABLE 2.12 EFFECT OF LORENTZ HALF-WIDTH ON APPARENT 
ABSORPTION OF FUNDAMENTAL (4.6 ym) 




R7 

PI 

P8 

P26 

STANDARD CO 

.06 

.2002 

.0810 

. 1700 

.0010 


.09 

.2387 

.0905 

.2012 

.0010 

LOW ALTITUDE 

.06 

.6620 

.2974 

. 5695 

.0167 

SOURCE 

.09 

.8052 

.3555 

. 6921 

.0178 

LOW CO 

.06 

.0439 

.0119 

.0350 

.0001 


.09 

.0470 

.0121 

.0370 

.0001 

HIGH ALTITUDE 

.06 

.2006 

.0807 

. 1703 

.0010 

SINK 

.09 

.2390 

.0898 

.2014 

.0010 


TABLE 2.13 EFFECT OF LORENTZ HALF-WIDTH ON 
ABSORPTION OVERTONE (2.3 ym) 




R7 

PI 

P8 

P26 

STANDARD CO 

.06 

.00647 

.00144 

.00462 

.0000098 


.09 

.00647 

.00142 

.00461 

.0000085 

LOW ALTITUDE 

.06 

. 1226 

.0299 

.0953 

.000431 

SOURCE 







.09 

.1312 

.0303 

.1005 

.000420 

LOW CO 

.06 

.000662 

.000146 

.000471 

.0000011 


.09 

.000649 

.000140 

.000457 

.0000006 

HIGH ALTITUDE 

.06 

.00631 

.00140 

.00451 

.0000097 

SINK 

.09 

.00631 

.00139 

.00450 

.0000083 



TABLE 2.14 EFFECT OF BANDPASS ON FRACTIONAL NET APPARENT. 
ABSORPTION - FUNDAMENTAL (4.6 ym) 


FUNDAMENTAL (4.6 ym) 




R7 

PI 

P8 

P26 


20 

. 1613 

.0677 

. 1328 

. 00049 

LOW 

TEMPERATURE 

2 

. 1593 

. 0670 

.1314 

. 00047 


0. 1 

. 1086 

. 0142 

, 0951 

. 00042 


20 

. 2002 

. 0810 

. 1700 

.00104 

NORMAL 

TEMPERATURE 

2 

. 1979 

. 0804 

. 1682 

. 00102 


0. 1 

. 1286 

.0635 

. 1155 

. 00093 


20 

. 1942 

. 0772 

. 1662 

. 00134 

HIGH 

TEMPERATURE 

2 

. 1921 

.0766 

. 1646 

. 00133 


0. 1 

. 1264 

. 0610 

. 1138 

. 00120 


20 

. 1270 

. 0521 

. 1063 

. 00043 

INVERSION 

LATER 

2 

. 1260 

.0518 

. 1056 

. 00041 


0. 1 

. 0977 

. 0450 

. 0852 

. 00040 



TABLE 2.15 EFFECT OF BANDPASS ON FRACTIONAL NET APPARENT 
ABSORPTION - OVERTONE (2.3 ym) 


OVERTONE (2.3 ym) 




R7 

PI 

P8 

P26 


20 

. 00660 

. 00152. 

. 00462 

. 0000062 

LOW 

TEMPERATURE 

2 

. 00656 

. 00150 

. 00459 

. 0000053 


0. 1 

. 00544 

. 00125 

. 00383 

.0000053 


20 

. 00647 

.00144 

. 00462 

. 0000098 

NORMAL 

TEMPERATURE 

2 

. 00643 

.00143 

. 00459 

.0000089 


0. 1 

.00531 

.00118 

. 00381 

.0000087 


20 

. 00637 

. 00140 

.00461 

.000014 

HIGH 

TEMPERATURE 

2 

.00633 

.00138 

. 00457 

.000012 


0. 1 

.00521 

.00114 

. 00378 

. 000012 


20 

. 00639 

.00140 

. 00461 

. 000013 

INVERSION 

LAYER 

2 

. 00635 

. 00139 

. 00458 

. 000012 


0. 1 

. 00523 

.00115 

. 00378 

.000011 
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the 4.6 ym is greatly affected by the temperature profile. This is as expected 
due to contribution of atmospheric emission. The 2.3 ym radiation is essential- 
ly unaffected., as is desired. Similarly the ground temperature and the ground 
emissivity, have larger effects (not desired) on the CO absorption signal at 4.6 
ym but not at 2.3 ym. 

The calculations show that, assuming measurement of CO absorption can be 
made on the 2.3 ym band, that band is much to be preferred over the 4.6 ym band 
for data interpretation. The advantages are: 

1. Measurements of temperature of the atmosphere, ground 
temperature, and emissivity are not needed for 2.3 ym 
but accurate values are required for 4.6 ym. 

2. The signal at 2.3 ym is affected significantly, as de- 
sired by the CO concentration in the lowest few kilo- 
meters in the atmosphere, but not significantly at 4.6 
ym and in a manner very difficult to interpret. 

3. The measurements at 2.3 ym are readily interpreted di- 
rectly in terms of CO densities whereas those at 4.6 
ym must be interpreted by use of an atmospheric model 
calculation at each point, assuming the required atmo- 
spheric data are available for the latter. 

Calculations using CO models 9 and 10 (Figure 2.1) were made of upwelling 
radiation in the 4.6 ym CO band for conditions as nearly as possible duplicating 
those used elsewhere (ref. 30) . Results of these calculations (ref. 5) show 
very good agreement. It can be concluded that the calculations reported here 
and in more detail (ref. 5) can be used for application to any instrument em- 
ploying radiation in these bands. 

2.6.3 Limb Inversion Analysis. - The limb experiment has been analyzed for 
the effects of instrumental error on the inversion of the instrument values of 
total CO in the path to yield CO concentration vs altitude. The calculations 
have been performed for the three CO sink models to which the limb experiment is 
relevant. These are the constant mixing ratio model with 0.1 ppm at ground lev- 
el, the upper atmosphere sink, and the low level sink with maximum concentration 
at 9 km. The error sources, which have been introduced, include a random error 
with standard deviation of 2%, 5%, or 10% and systematic errors of ± 2%, ± 5%, 
or ± 10%. 


For these calculations the atmosphere below an altitude, hj, was divided 

into n speherical shells (j = 1, 2, n) , a distance. Ah, apart (Figure 2.16). 

The atmospheric properties within a shell are assumed to be constant. The alti- 
tude of the first shell is given by 

hj = nAh, 

and the altitude of the jth shell is 

hj = hi - (j-1) Ah 
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The ith ray from the sun to the instrument passes horizontally through the 
lower edge of the ith shell (j = i) . The height of this ray is 

h . . = h . - Ah 
11 1 


The rays passing through the shells can be divided into lengths, a^^, over 
which the atmospheric properties are assumed to be constant. The lengths a. . 


are given by: 


11 


.. =|(r. h.)2 . (r.h.p2| 

^ij = ^ h.)2 - Cr + h.p2| 


1/2 


1 2 


if i = j 


k = i 

Z 

k = j + 1 


a. 1 

ik 


if i ± j 


A square matrix (n x n) of the elements a^^ can be formed. 

In a shell j, the average concentration of CO is given by: 

c. = cfh. - 1/2 Ah) 

J 1 

The integrated amount of gas in a horizontal path is given by: 

{u^} = 2{c^}[a^^] 

where 


(up 

E column 

(cp 

= column 

ta,.l 

= square 

The calculations were 


2 km distance between shells. The matrix elements u. were calculated for the 

1 

three CO models (1, 4, and 5) in Figure 2.1. The correlation interferometer 
gives an instrument output which is proportional to the total CO in the path, 
or u^. In this case the known parameters will be u as a function of altitude, 

and the a^ . elements will be known. A set of equations can then be set up to 

solve for the concentrations c. . 

1 
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where is the inverted matrix [a..] 

1] '■ ij-* 


To test for the effects of random or systematic instrument errors, the 
values of u were computed for each of the three concentration models from the 
relation {u.} = 2{c.}[a. .]. The computed values were then changed by fixed 
amounts ^ J 


a) random error with 2%, 5%, or 10% standard deviation 

b) systematic error of ±2%, ±5%, or ±10%. 

For each of these types of errors, the concentration profile was calculated and 
compared with the exact profile used as the original input. Some results of 
the calculation are shown in Figure 2.17. 

With measurement errors of the order of 5 and 10%, the computations show 
that it is possible to distinguish the carbon monoxide constant mixing ratio 
concentration profile from the concentration profile which assumes a high alti- 
tude sink. The computations with a 2% random error show the same result with a 
much smaller error. 

The computations indicate that instrumental errors of the order of 10% in 
the measurement of the total CO in the horizontal sight path through the atmos- 
phere do not result in significantly larger errors in the computed CO concentra- 
tion profile at high altitudes. While the measurement errors are more signifi- 
cant in the case where an inversion in the CO concentration profile occurs (i.e., 
the low-level sink) , they do not invalidate the limb transmission experiment at 
the higher altitudes. The limb experiment was primarily intended for searching 
for an upper atmosphere sink, and in this case the effects of reasonable mea- 
surement errors are not of sufficient magnitude to permit the upper atmosphere 
sink profile to be mistaken for the standard constant mixing ratio CO profile. 

2.6.4 Multi-Line Model Calculations. - 


2.6.4. 1 The program: The first step in the theoretical determination of 

the feasibility of using the correlation interferometric technique for the mea- 
surement of CO levels in the atmosphere is the computation of theoretical trans- 
mission spectra of the atmosphere in the wavelength regions of interest. Theo- 
retical spectra were calculated by a fairly elaborate computer program, devel- 
oped by GE. Tliis program has been described in greater detail elsewhere (ref. 
5). 

Additional subroutines, patterned after Cooley (ref. 31) and Gentleman 
(ref. 32) then compute the Fourier transforms of the resultant spectra. These 
Fourier transforms which are closely related to interferograms form the data 
base for the theoretical feasibility study. 

The computer program was used to generate a large number of theoretical 
spectra and transforms of CO in the presence of H 2 O. Spectra were generated 
in the 2.3 nm spectral region (using geometry case 2) and in the 4.6 urn region 
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Figure 2.17 Computed Atmospheric CO Concentrations with Simulated Instrument Error 



(using geometry case 1). The regions covered were 4230 - 4330 cm“^ and 2100 - 
2200 cm”^, respectively. Spectra included combinations of CO models 1, 2, 3, 
and 4, as described in Section 2. Ground emissivity was varied in the 4.6 yn 
spectra, and albedo was varied in the 2.3 ym spectra. In all, 53 theoretical 
spectra and transforms were computed at 2.3 ym and 19 at 4.6 ym. Although there 
are many additional combinations which would have been desirable, extremely long 
running time of the program (about 10 - 15 minutes on a GE 635 or IBM 360) makes 
a full gamut of combinations prohibitively costly. 

Following generation of these theoretical Fourier transforms, analyses were 
performed by use of a weighting function computer program. This program, (ref. 
5) patterned after a program furnished by Barringer Research, is an implementa- 
tion of the weighting function concept described in detail in Section 6. 

The weighting function can be considered a correlation function which, when 
multiplied by the Fourier transform of a spectrum gives a measurement of the CO 
level . Inputs to the weighting function program consist of a basic set of sev- 
eral Fourier transforms of spectra generated from a constant amount of CO and 
various combinations of other variables, e.g., temperature profile, H 2 O level, 
etc. One of these transforms is designated as the nominal case. One additional 
transform is input, designated as the target case, which is identical to the 
nominal case with the exception of a variation in the CO level. The weighting 
function generated from these cases is then applied to other transforms which 
were not included in the basic set. It is this operation of generating weight- 
ing functions from one set of data and applying them to another set which de- 
termines the feasibility of the correlation operation. 

Tne theoretical spectra Fourier transforms described above were used in 
various combinations as basic sets for the weighting function program. Two 
groups of studies were made, at 2.3 ym and 4.6 ym. The results of these studies 
are described below. 

Calculations were made using the principles outlined in Sections 5 and 6. 
The intent of these calculations is to determine the sensitivity of the method 
in the determination of CO densities, the accuracy to be expected in the mea- 
surements, the effects of various atmospheric parameters, the optimum wavelength 
region to be used, and in general, to establish the feasibility of the technique 
in obtaining the data needed to determine the CO sink. 

While the theoretical feasibility studies will be necessary, the establish- 
ment of the feasibility will only be accomplished after the instrument is shown 
to be capable of measuring CO in the atmosphere. 

2. 6. 4. 2 Calculations: It must be remembered in examining the results that 

in the measurement of CO burdens as such or for their measurement to find the CO 
sink, it is important to be able to see an effect produced by a rather small 
variation from normal CO column density, since, although at ground level the 
concentration may drop off by a large factor, the total column density may only 
have small drops, as little as about 10%. Further non- sink areas should not ap- 
pear as sinks and the accuracy for very low CO column density may be poor, as 
long as it is good enough to show it to be less than the normal atmospheric 
amount. It must be emphasized that the most accurate measurements are needed 
at densities near those of a normal CO profile. 
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Further, the results should not be affected significantly by variations in 
atmospheric or ground properties, e.g., atmospheric temperature profile, temper- 
ature inversion layer, ground temperature, ground emissivity, reflectivity, and 
the shape of the CO profile. 

The results which are presented in the following table are given as an ex- 
ample of the results of the many similar calculations made. This table (2.14) 
gives the results for all interferograms , using a specified set of interfero- 
grams to determine the weighting functions and then using these weighting func- 
tions to determine the CO density in all the interferograms. The table gives 
the percentage difference of the calculated value from the actual amount of CO 
in the model. The four left-hand columns give data on this model, the first 
column giving the actual amount of CO and the second column giving the CO pro- 
file shape, the CO model numbers being those of Figure 2.1. The third column 
gives the total water burden and the fourth column gives the water profile 
shape. Water models correspond to water models 1, 2, and 3 of Figure 2.4, mul- 
tiplied at all altitudes by a constant factor. The middle four columns are for 
the 2.3 ym (overtone) band and the right four columns are for the 4.6 ym (fun- 
damental) band. These last eight column headings note the four corresponding 
temperature models of Figure 2.5. 

The table presents the results using a specific set of runs to obtain the 
weighting functions and gives the accuracy of the calculated CO for each inter- 
ferogram of either the 2.3 or 4.6 ym band as noted. The sets used in obtaining 
the weighting functions are indicated by the x's. 

Some of the results presented in Table 2.16 show accuracy which will fit 
the accuracy requirements detailed above while some show inaccuracies which are 
too large to be acceptable. This is as expected. The calculated accuracy of 
the calculations are dependent on the conditions which were used for determining 
the weighting functions. It is evident that, for the overtone band, the results 
are quite acceptable. The data of this table are typical of the many calcula- 
tions carried out on this band. The calculated errors near the standard pro- 
file, that is for burdens near that of the normal clean atmosphere, are less 
than 5% for all temperature models. Although the errors are greater for the 
very small and very large CO burdens they are certainly not so great that their 
size is not readily distinguishable from that of the standard. This is not the 
case for the fundamental band. Here small and large burdens give excessive er- 
rors which make them, in many cases, indistinguishable from the standard. And 
for the case of a temperature inversion the errors are extremely large. This 
example for the fundamental band shows smaller errors than most other calcula- 
tions made on that band. Such calculations were so unfavorable that many fewer 
calculations were made for this band than for the overtone band. 

Many more conditions were calculated and are reported in detail elsewhere 
(ref. 5). All show similar results - quite useable for 2.3 ym and very poor 
for 4.6 ym. 

Some general conclusions cari- be drawn from the results for the 2.3 ym band. 
If the conditions used in obtaining weighting functions for CO are within the 
range of CO where absorption is about linear, that is within the range from very 
low CO densities to slightly above our standard model, the results will be good 
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TABLE 2.16 CALCULATED ERRORS IN CO MEASUREMENT 


% Error (4.6 ym meas.) 
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within the range of water density and temperature covered in the weighting func- 
tion runs and probably not good outside such range. That is, if water density 
models 1 and 3 are used, those for in-between water densities have the desirable 
accuracy, but if models 2 and 3 are used, the results for water model 1 are in- 
accurate. That is, interpolation gives reasonable results but extrapolation 
does not. 

In these models the total water content was varied by a factor of fifteen. 
Such a variation in total column density is probably as much as the variation 
over the year for any one region. Seasonal variations for a region range up to 
a factor of four (ref. 33) while shorter term variations may be somewhat larger. 
Hence, in treating the data for any one area, weighting functions based on the 
range of conditions for that area can be used rather than those which bracket 
the range of conditions for all areas. 

Conslusions about sinks and sources can also be made. These results show 
that sinks (even those at low altitude having only a 9% lower CO optical thick- 
ness than the standard model) are detected readily and regions of high CO con- 
centration are also readily seen. It is also to be noted that the presence of 
a temperature inversion layer does not interfere with the results for this (2.3 
urn ) band. The effects of surface reflectivity, temperature, and emissivity 
and of atmospheric temperature are negligible. 

In all 4.6 pro cases, the source runs show up with large negative errors. 

In almost all urban model runs, negative CO densities are calculated. These are 
meaningless and undoubtedly smaller sources would show up as sinks. This is an 
intolerable situation. 

Certain other tests were made to help show the effects of atmospheric and 
earth-surface parameter variations. Results show that surface temperature and 
emissivity and atmospheric pressure all appreciably affect the results and would 
have to be included in the interpretation of any measurements with the 4.6 pm 
band . 


2.6.5 Conclusions . - The accurate measurement of atmospheric trace species 
by the correlation interferometric technique has been shown by analysis to be 
feasible. Specifically, it has been shown to measure atmospheric amovints of 
carbon monoxide with an accuracy of about 10%. Thus the technique should be ca- 
pable of picking out CO sinks which are only about 10% lower in CO column den- 
sity using measurements which are made by determining the absorption in a part 
of the first overtone band of CO in the 2.3 pm region of the spectrum. It has 
been determined that the use of this band is preferable to the use of the funda- 
mental band in the 4.6 pm region. The latter is unsuitable because of the ef- 
fect of atmospheric parameters including the atmospheric temperature profile 
(and the associated atmospheric emission), the ground temperature, the ground 
emissivity at wavelengths being used, atmospheric pressure and atmospheric path 
length. For that band, atmospheric emission causes an enhancement of the radi- 
ation which is a function of temperature and must be taken into account. If at 
some altitude the atmospheric temperature is greater than the ground temperature 
the effect of absorption below is minimized. Effectively, the instrument cannot 
see well below an atmospheric temperature peak. Then the absorption does not 
follow the CO profile at low altitude but rather there is less net absorption 
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per CO molecule near the ground even for CO profiles when there is more CO in 
that region (see Figures 2.14 and 2.15). These difficulties arise for any spe- 
cies with all spectral techniques which are primarily involved with radiation 
of wavelengths greater than about 3.5 ym. This is because of the predominance 
of earthshine over reflected solar radiation at these wavelengths (see Figure 
2.6). The overtone band is not affected significantly by variations in the pa- 
rameters noted above and absorption in this band follows closely the atmospheric 
CO profile for any reasonable atmospheric model. Atmospheric scintillations 
have little effect at either wavelength region. 

The conclusions for the two wavelengths studied may be summarized on the 
basis of the following considerations. 

Optical Thickness: 

4.6: For much of the range, the absorption curve is not linear. 

The higher optical thicknesses may be so high that there is 
little sensitivity in that range. 

2.3: The optical thickness is suitable for the range of concen- 

trations of interest. The only limit is the upper limit of 
the limb experiment but it will be suitable for a sink at 
any reasonable expected altitude. 

Interpretation: 

4.6: Very difficult because of the need for accurate atmospheric 

temperature data as a function of altitude, ground tempera- 
ture, and ground emissivity at 4.6 ym. 

2.3: Data directly presented in terms of CO density. 

Sensitivity to Low Altitude Sink: 

4.6: Sensitivity to low altitude CO is very low and errors intro- 

duced are likely to be larger than the CO decrease effect 
being sought. 

2.3: Sensitivity at low altitude same as the similar effect at 

high altitude. Sinks of less than 10% CO density decrease 
appear to be detectable. 

Temperature Inversion Layer: 

4.6: Drastic effects on calculated CO densities prohibit detec- 

tion of significant sinks. 

2.3: No appreciable effect. 

Line Strengths : 

4.6: Strong. 

2.3: About 1% of 4.6 lines. This is strong enough. 
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Atmospheric Emission: 

4.6: Considerable emission in this band causes much variation 

from straight absorption model and leads to difficulty in 
interpretation as mentioned. 

2.3: Emission negligible. 

Night Use: 

4.6: Possible. 

2. 3: Not possible. 

Mapping Experiment: 

4.6: Emission effects disastrous. 

2.3: Interferents cause problems but these are overcome by tech- 

nique used. 

Limb Experiment: 

4.6; Atmosphere optically thick at lower altitudes of interest. 

2.3: Sensitivity limits altitude but altitude range reasonable. 

Spectral interferents can be overcome by the correlation interferometry 
technique. The chief spectral interferents in the 2.3 ym spectral region are 
water and methane. By using the interferogram directly, carbon monoxide can be 
accurately measured in the presence of atmospheric amounts of these gases which 
would prevent accurate measurement of CO by ordinary spectral methods . In order 
to accomplish the interferometric measurement it is necessary to calibrate the 
instrument over the entire range of density of interferents for which it will be 
used and to do so with variation of important atmospheric parameters over which 
it will be used. Thus the calibration must cover the range of methane of about 
2 to 5 atmosphere cm, since this is the range that is expected to be encoun- 
tered, and to cover the range of about 0.2 to 3 percipitable cm of water, since 
variations from dry to wet atmospheres include this range. Further, since the 
population of the rotational water levels are appreciably affected by tempera- 
ture, variations corresponding to changes in atmospheric temperature profiles 
must be included in the calibration, thus necessitating calibration for condi- 
tions of a cold dry, a hot dry, a cold wet, and a hot wet atmosphere as well as 
conditions in between these. 

The calibration determines the weighting function which multiplies the sec- 
tion of the interferogram in such a way as to minimize the effect of spectral 
interferents and maximize the effect of the gas to be measured. The choice of 
conditions used to determine the weighting functions is critical. The use of a 
wide range of conditions and interpolation between these gives much better re- 
sults than use of a narrow range with extrapolation. In practice under flight 
conditions, it may be best to use a weighting function derived from a wide range 
of conditions to obtain an approximate measure of CO density and conditions and 
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then use a weighting function derived from a narrower range of density and con- 
ditions to obtain a more accurate CO density measurement. 

The feasibility analysis thus shows the correlation interferometry tech- 
nique to be capable of the measurement from a remote platform of carbon monoxide 
in the atmosphere over the desired range of density using the first overtone 
band. The measurement can be made in the mapping mode (observing sunlight re- 
flected by the earth) and the limb mode (observing direct sunlight through the 
earth's limb). 
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3. CORRELATION INTERFEROMETRY 


3.1 Principles of Interferometry 

The use of spectral techniques for the remote measurement of concentrations 
of trace atmospheric species is dependent on separating out the effects of the 
species being measured from those of all other species present in the optical 
path. In techniques where a part of the spectrum is measured, the separation 
must be obtained by spectral resolution. Thus, some separable part of the spec- 
trum must show an appreciable effect of the species being measured and any sig- 
nificant effects of other species must be such that they can be eliminated. If 
a different technique is used effects must still be separable but the separation 
criteria is no longer spectral resolution but a type of resolution peculiar to 
that technique. Such a technique is interferometry. In this technique the sep- 
aration is accomplished by a resolution of path differences, provided the data 
are obtained directly from the interferogram rather than from a spectrum ob- 
tained by transforming the interferogram. The instrument developed for this 
program is a correlation interferometer (ref. 34) . The following discussion 
will describe the basic theory of its operation. 

An interferometer is a fairly simple device (Figure 3.1) . The essential 
elements are a beam splitter and two mirrors, plus a detector to measure the ra- 
diation output. Light from the source is incident on the beam splitter, B. At . 
the beam splitter it is divided into two paths; one portion of the light goes 
to one mirror, M^; the other portion of the light goes to the other mirror, M 2 . 
The two portions recombine at the beam splitter and the intensity of the light 
once they recombine is registered by the detector, D. The intensity of the ra- 
diation received by the detector will depend on the difference between the 
lengths of the paths traveled by the beams in the two arms. The length of the 
Path F'B-Mi-B-D can be different from that of the path F-B-M 2 -B-D. If the two 
optical paths are exactly the same, the path difference (delay*) is zero, and 
there is a peak in intensity. If monochromatic radiation enters the instrument, 
and if the path difference is increased by one-half the wavelength, the inten- 
sity reaching the detector goes down essentially to zero. Then as it increases 
again towards one wavelength path difference, another peak occurs. This sinus- 
oidal oscillation about a mean level repeats at intervals of one wavelength if 
the light is monochromatic. The instrument actually does a Fourier transforma- 
tion of the spectrum of the radiation entering. The way in which the delay is 
scanned, that is the way in which the path difference is changed, is in most in- 
struments a matter of shifting one of the mirrors. One of the features of the 
COPE. instrument is that the problem of having to scan and maintain the position 
of the mirror accurately is avoided by not scanning the end mirror, but, in- 
stead, scanning a plate of refractive material in one arm of the interferometer 
which is generally left in a fixed position. If this plate is rotated, the path 
length in that arm of the interferometer will vary. This will accomplish the 
same effect as moving one of the mirrors back and forth without the alignment 
problems. This is a specific advantage of the technique of this instrument. 

*The term "delay” refers to the temporal variation of the interferometer sig- 
nal. This is related to the variation of the path difference by the scan 
velocity. 
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Other advantages, some of which are advantages of interferometry in general, in- 
clude a large light throughput, multiplexing of spectral elements, a compact yet 
flexible instrument, and a handy output. 


3.2 Relationship Between Spectra and Inter ferograms 

In Figure 3.2 a few particular cases of the relationship between spectra 
and Interferograms are shown. As shown on the previous figure, monochromatic 
spectrum produces an interferogram which is essentially a sinusoidal variation 
as the delay is scanned. If instead of a single monochromatic line there is a 
pair of monochromatic lines, then the two sinusoidals due to lines will beat to- 
gether. Again a characteristic maximum is obtained at zero delay. As the os- 
cillations beat together, they get out of phase and go through a minimxjm, then 
maximize again at a point characteristic of the line separation. If there are 
a number of regularly spaced lines, the beat patterns will all reinforce at zero 
delay again, but they very rapidly get out of step with each other and decay to 
very small values. Then there appears a point where they all reinforce again at 
a delay which is inversely proportional to the spacing between the lines. This 
is almost exactly the effect seen in the case of carbon monoxide as shown in 
Figure 3.3, which gives the actual spectrum* and interferogram. The spectrum 
displays fairly regularly spaced lines. There is certainly a noticeable change 
in the spacing from one end to the other, but it is certainly not a random spac- 
ing. As seen in the previous figure, there is reinforcement for zero delay 
which dies out rather rapidly, and then, at a delay which is characteristic of 
the spacing between the lines (about 3 cm~^), the interferogram amplitude peaks 
up again. In general the relationship between the spectrum of the radiation 
and the interferogram of the radiation is given by a Fourier transformation, the 
cosine Fourier transformation. That is, the interferogram, apart from the con- 
stant term, is the cosine Fourier transform of the spectrum (Figure 3.4) . The 
interferogram signal, as a function of path difference is given by: 

I(X) = 2/” (S^(o)) cos (2ttox) da 

where 

a = frequency (cm”^) ^ 

S^(a) = spectral input at a (ergs-cm^-sr"^ (cm- -sec"^) 

X = path difference (cm) 

Several relationships exist between the interferogram and the spectrum 
from which it is produced. The interferogram has a carrier frequency directly 
related to the mean frequency of the spectrum producing it. It will be ampli- 
tude and phase modulated, the amplitude and phase modulation being described by 
the band envelope shown. From an interferogram up to a given delay, the spec- 
trum could be reconstructed with a resolution given inversely by the delay out 
to which the interferogram is obtained. So, from an interferogram from zero to 
1 centimeter, inverse Fourier transformation reproduces a spectrum to resolu- 

*This spectrum is that given by Plyler (ref. 17) with the mercury emission 
lines omitted. 
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Figure 3.3 CO Spectrum and Interferogram 

(of part of P branch of spectral band) 




SPECTRUM 



Figure 3.4 "Heterodyning" o£ Interferogram 



tion of one wave number (cm"^). Conversely, if the information characteristic 
of one wave number resolution is not needed, if the information need only be 
characterized by a 3 cm~^ resolution, then it is only necessary to scan the in- 
terferogram out to 1/3 cm (3 millimeters) to obtain information equivalent to a 
resolution of three wave numbers. It is very important to do this. Significant 
improvements in signal to noise ratio are obtained by limiting ourselves to 
scanning only the portion of the interferogram which gives the best signal to 
noise ratio. If essentially all the information on any given species (all the 
effect of that species on the interferogram) occurs over a small part of the in- 
terferogram, only that range of path difference need be scanned. The operation 
of the correlation interferometer involves the treatment of the interferogram 
data directly to obtain data on species densities rather than the use of the 
spectrum obtained by the Fourier transform of the interferogram. With such an 
operation the concept of spectral resolution loses meaning. Although the spec- 
trum which can be obtained from an interferogram has a resolution which is the 
reciprocal of the path difference scanned, the ability to obtain data on one 
species in the presence of others by direct examination of the interferogram, is 
not simply dependent on the length of path-difference scan. 

Probably the main objection to interferometers is that most people think in 
terms of using a spectrum and in order to obtain a spectrum from the interferom- 
eter, the interferogram must be transformed. This problem is avoided completely 
by not looking at the spectrum at all. It is not really necessary to use the 
spectrum. The measurement can be made quite adequately on the interferogram 
itself. This is a major advantage of the correlation interferometer. 

Another relationship between the spectrum and the interferogram is that the 
spectriom of the incoming radiation can be severely band- limited by optical fil- 
ters. If it is band- limited, and if an interferogram that would reproduce the 
spectrum only to a limited resolution is taken, then the spectrum could be char- 
acterized by 'a farily small number of points. On the other hand, to character- 
ize the interferogram in this form with its very rapid oscillations (the oscil- 
lations are essentially proportional to the mean frequency and not to the width) 
with any degree of accuracy, a relatively large number of points, on the order 
of at least one per cycle of oscillation, would be needed. This would be many 
more points than are necessary to characterize the spectrum to a corresponding 
resolution. There is a lot of redundancy in the information because of the 
band limitation. The method used in the correlation interferometer to elimi- 
nate the redundance is simply to take the amplitude and phase modulated sinus- 
oid and heterodyne it down with a local oscillator, where the local oscillator 
is the interferogram of radiation somewhere around the mean frequency of the 
spectrum. By heterodyning the interferogram down with a cosine or a sinusoidal 
variation, the interferogram is reduced to its essential variations. This is 
illustrated in the bottom part of Figure 3.4. All the information necessary to 
characterize the interferogram can be retrieved by sampling a relatively small 
number of points, a few points for each of the much longer cycles shown in the 
bottom figure. There is a slight difference whether the interferogram is beat 
with a sinusoid at a given wavelength or a cosine; slightly different beat-down 
interferograms are obtained. Actually both of them are used but, through most 
of the remaining discussion this fact will be ignored. In generating the local 
oscillator for carrying out the heterodyning, radiation which passes through the 
same interferometer as does the signal radiation is used. This insures that the 
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local oscillator always has the correct phase relationship with the signal in- 
terferogram. It also relaxes some of the accuracy needed in the knowledge of 
the plate drive (the scan drive) . If this were not done, the scan drive would 
have to be known accurately to within a very small fraction of the wavelength. 

By having the local oscillator going through the same interferometer, this prob- 
lem is alleviated. 


3.3 The Measurement in the Presence of Interferents 

Figure 3.5 outlines very briefly the basic principles by which the measure- 
ment on the interferogram is made. There are essentially two problems involved 
in a measurement of this type. We are trying to measure the CO burden on the 
basis of the radiation received from a satellite. The problems in making the 
measurement are: (1) can the measurement be made with sufficient accuracy in 

view of the noise limitations, and (2) can the measurement be made under condi- 
tions where the radiation received is affected not only by the gas that we are 
trying to measure, in this case carbon monoxide, but is also severely affected 
by other gases. In fact, the radiation is much more affected by such things as 
water vapor and methane. Ignoring for the moment this problem of the interfer- 
ence species, consider how a measurement would be made. Consider an interfero- 
gram due to a target signal, say CO, such as shown. There are different ways 
by which one could measure the CO burden which produced this signal level. One 
might sit at a constant delay and measure the signal level at that point. If 
the CO burden were doubled, the effect of CO on the intensity we would measure 
would be doubled (for the moment assuming variations are linear with the gas 
burden); the intensity would essentially be a measure of the carbon monoxide 
burden. If there were a constant noise level associated with the measurement of 
any point then the best point to make a measurement would be at the point where 
the signal level is maximum. If a finite range of delay were scanned, then the 
measurements made at all points would be combined to get some sort of average 
measurement of the carbon monoxide burden. In general, the optimum measurement 
that can be made in such circumstances is given by combining all the measure- 
ments, all the intensities of the various points, in a manner related directly 
to the intensity of the signal shape. That is a weighting function or correla- 
tion function is generated and this is multiplied together with the signal. The 
measurement made is the integral of the correlation function times the signal, 
the integral over the delay range which is scanned. It is integrated and the 
measurement obtained is directly proportional to, in this case, the carbon mon- 
oxide burden. The signal-to-noise ratio in such a measurement can be shown to 
be optimum when the correlation function looks exactly in shape like the target 
signal. That result holds when there are no interferents. Now consider another 
gas species which is affecting the signal received, for example, water vapor. 

If a point measurement of the intensity were used, the signal received would not 
be simply the signal due to carbon monoxide but the sum of the signals due to 
the carbon monoxide and to the water vapor. Even if the carbon monoxide level 
were to remain the same, drastically varying measurements due to variations 
in water vapor might be obtained so that a very poor measurement of CO would be 
obtained by sitting at that point. A region where a signal is unaffected by all 
interferents may not be available or if it is, the signal-to-noise ratio achiev- 
able at the point might not be adequate for the measurement. This problem is 
handled as follows. The correlation function used is not matched exactly to the 
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Figure 3.5 Application of Correlation Function 



target gas signal shape. The correlation function (W) is adjusted so that when 
it is cross-correlated with the interferents and the result is integrated over 
the range scanned, all the positive correlation regions are balanced exactly by 
the negative correlation regions, so that the total area under the curve, the 
total result of the measurement, comes out to zero as illustrated, still main- 
taining the correlation function as close as possible to that constraint to the 
signal interferogram (target interferogram) in order to still get as large a 
positive correlation between the correlation function and the target gas as pos- 
sible. The result of the measurement is still proportional to the carbon monox- 
ide burden. In principle, if there are a number of interferents of this sort, 
rejection of these interferents can still be achieved so long as we have at 
least as many points to describe our correlation function as we have gases which 
are affecting the radiation. The measurement of the target gas can be made in 
real time. The first gas, in this case CO, and any other gases which signifi- 
cantly affect the interferogram can be measured. 

In the correlation interferometer, the path difference scan is limited 
principally to that portion of the interferogram which is most affected by the 
gas whose density is being measured (excluding that portion near zero path dif- 
ference) . In the case of CO that portion of the interferogram obtained with 
path differences between 2.70 and 3.95 mm has been selected. Since the lines 
in the first overtone band of CO are separated by about 3 cm-^ on the average 
this shows up in that part of the interferogram centered at about 1/3 cm"^ = .33 
cm = 3.3 mm. Thus, the CO spectrum is as shown in Figure 3.3 and the interfer- 
ogram (using a spectral filter centered at 4287 cm-^) is also shown in Figure 
3.3. In this the major effect of the CO is seen to be in the region of 0.3 cm. 
Taking this portion of the interferogram between 0.27 and 0.395 cm gives essen- 
tially all of the information on CO that can be obtained from the interferogram. 
The determination of CO is not, however, the simple measurement of this peak in 
the interferogram since other gases, particularly CH 4 and H 2 O also have some ef- 
fect in this region. Thus, this peak is divided up into a number of sections, 

32 in the case of our instrument. The effect of each gas on each of these 32 
sections is determined in effect, by calibration.* 

Thus, in the simplest case, a portion of the interferogram affected only 
by H 2 O is chosen to determine the simount of H 2 O** and the remainder of the in- 
terferogram is, in effect, corrected using these data and calibration informa- 
tion. CH 4 is dealt with similarly. Thirty-two different effects could be dealt 
with this way. Since in a real case portions of the interferogram affected by a 
single gas are not always available, the set of 32 simultaneous equations is 
solved to determine the amounts of individual gases causing the combined effects. 

*Actually, the effects of specified gases, other than that being measured (such 
as CO), need not be specifically known by calibration. Tlie calibration can 
' consist of 31 different atmospheric conditions with one amount of CO and a dif- 

ferent amount of CO with one of these conditions. The measurement of CO is 
then made by determining the effect on each section of the interferogram of 
other gases of the atmosphere as determined by a best fit of the 32 points. 

**Since the H 2 O spectrum and hence its interferogram is affected by temperature, 
two or more portions of the interferogram must be used to take this effect 
into account. 
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The output of the interferometer will be voltages, each related to the con- 
centration (column density) of one species. The calibrations can be thought in 
terms of these voltages. Thus, to be able to separate out the effects of indi- 
vidual species it is necessary to resolve parts of the interferogram which show 
the effects of these species. This, resolution is not attainable with the spec- 
trum which would be produced from the interferogram obtained with this instru- 
ment. Thus, by directly using the interferogram the effects of specific species 
and hence its density is measurable whereas it would not be if the interferogram 
were converted to a spectrum and that used for the measurement. This compara- 
tively simple measurement can be made very much faster than can a measurement of 
species by a high-resolution spectrometric technique, whether by a spectrometer 
technique or by an interferometric technique which employs the Fourier transform 
of a interferogram obtained over a large range of path difference. 


69 



Page Intentionally Left Blank 



4. 


CORRliLATION TRIiATMIiN'I' 


4.1 Basic Principles 

4.1.1 Instrument Output . - The interferometer receives radiation, the 
spectral characteristics of which arc affected by a number of parameters. We 
wisli to measure one of these parameters (CO burden) , with as high accuracy as 
possible and as little interference as possible from variation in the other pa- 
rameters. Within the instrument, optical filters limit the range of the radia- 
tion to a fairly narrow spectral region where CO absorption has the greatest 
relative effect on the radiation. 

The correlation interferometer operates by varying the delay between the 
two beams of a Michelson interferometer over a defined range. As the delay is 
swept, a sinusoidal type signal is generated at the detector output. This sig- 
nal is both amplitude and phase modulated. The band center and width are di- 
rectly related to the center and width of the optical spectral filter, while 
the modulations are characteristic of the more detailed spectral information. 

4.1.2 Preprocessing . - An important function of the interferometer elec- 
tronics hardware is to reduce the high sampling rate implied by the relatively 
high center frequency of the interferogram signal, by synchronously detecting 
the signal with a reference signal, (similar to the local oscillator in a homo- 
dyne system) . This reference is generated from the interferogram of radiation 
having stable spectral characteristics. The reference radiation passes through 
the same interferometer system as the signal radiation, so that to a very large 
extent the reference derived from its interferogram is fixed in phase relative 
to the interferogram of any given signal radiation, regardless of any irregular 
ities in scan, or variations in the interferometer arm lengths. 

When this synchronous detection is carried out, the result is a signal con 
taining only the information of the sidebands of the original signal, beat down 
to zero frequency means. (Actually two such signals are produced, from refer- 
ences both in phase and in quadrature phase with the reference interferogram.) 
The signal's information is still band limited, by virtue of the original spec- 
tral band limitations. The signal is integrated and digitized at an interval 
sufficient to retrieve this information, producing a preprocessed, digitized, 
interferogram. 

4.1.3 Final Processing .- Final processing consists of applying a linear 
digital filter ('weighting function' or 'weights') to the preprocessed interfer 
ogram. The filtering operation consists of taking the weights, one for each 
point of the interferogram, multiplying them together with the corresponding in 
terferogram points, and summing the results. The weights are chosen, in a man- 
ner to be described below, so as to give a final result which is insensitive to 
variations in all parameters except the desired one. (The weights are also cho 
sen to give a result which is optimized against a mixture of noises.) The re- 
sult may be directly converted using a defined zero point and scale factor, to 
units of target gas burden (atm-cm, ppm-m) , etc. 
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4.2 Determination of Weights 

4.2.1 Basic Philosophy and Theory .- The determination of the weights is 
based on the possibility of representing any observed interferogram as a linear 
combination of component interferograms . The number of such component interfer- 
ograms required to represent a range of actual observed interferograms to some 
accuracy, will depend on a number of factors depending on the range of condi- 
tions under which the interferograms are observed. The number will depend main- 
ly on the number of parameters which vary significantly over the range of condi- 
tions, and to a lesser extent on the degree of non-linearity in the variations 
of the interferograms with these parameters. (Parameters will characteristical- 
ly represent factors such as water vapor burden, methane burden, and temperature 
profile, as well as the target CO burden). 

Consider then any interferogram, I, within the defined range. It may be 
written as a linear combination of N constituent interferograms. Using the in- 
dex k to define the point on the interferogram, and j to define the particular 
constituent interferogram we may write 

N 

I(k) = Z qj Ij(k) 

j = 1 

The qj's are the strengths of the constituent interferograms in the particular 

interferogram observed. Let us assume one of these, q.., is the target gas bur- 
den. 


Now the final processing consists of the application of a weighting func- 
tion, H(k), to the interferogram. The resulting measurement is 

M = z H(k) I(k) 
k 

= Z q. L H(k) I (k) 
j J k ^ 


From this, it can be seen that if we can find weights H(k) that give 


Nk 

Z H(k) I (k) =1 j = jl 

^ ^ 0 j jl 


Then the measurement 


will be M = 



the target burden. 


If 


varies from 1 to 


N, then this is a set of N equations in as many unknowns, H(k), as there are in- 
terferogram points (Nk) . In order for these equations to be satis fiable, there 
must be at least as many interferogram points as there are constituents (Nk ^ N). 
In the case that Nk = N, there are N linear equations in N unknowns, and the 
solution for the H(k) is straightforward. 


In the case that there are more interferogram points than there are con- 
stituent interferograms, there are many possible weighting functions which will 
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satisfy the N equations. It is, however, possible to choose a unique set of 
weights by considering the noises in making the measurement, and minimizing 
their effect on the measurement. 

4.2.2 Noises in Measurement .- The noises operate in the following manner. 
The total integration time for measurement of a given point may be written D(k) 
(its duration, assumed unity up till now) . The contribution to the final mea- 
surement for that point is actually AI4(k) . 

AM(k) = [H(k) D(k)] I(k) 

The noise contribution for that point is significant in how it can vary this 
from its "correct" value. The noises may or may not decrease with integration 
time; that is, they may be either random or synchronous with the scan. They may 
or may not depend on the expected level of the interferogram itself (an 

RMS value for the constituent interferogram at a given point); that is, they may 
be either additive or multiplicative. 

A discussion of various physical origins for these noises will be presented 
in Section 4.2.4. However, it is sufficient at this point to be able to charac- 
terize them by four numbers, ai....a 4 , which represent the severities of the 
four possible combinations of noise type: 


ai - 

Random Additive 

(RA) 

82 - 

Random Multiplicative 

(RM) 

as - 

Synchronous Addivive 

(SA) 

84 - 

Synchronous Multiplicative 

(SM) 


The noises will combine in an RMS manner to give the mean square noise contribu- 
tion AN^(k) to the measurement at a particular interferogram point, and we may 
write 


AN2(k) = [H(k) D(k)]2 X [Cai + 82 <I>2)/D(k) + (ag + 84 <I>2)] 

= [H(k) D(k)]2 G(k). 

G(k) is an effective mean square interferogram error for point k (although it 
may originate through error either H, D, or I) . 

4.2.3 Derivation of Weights . - 

The total measurement is thus 

M = Z AM^(k) = E [HD] I 
k k 

Q 

and the total squared noise is 

o n 2 = E AN2(k) = E [HD] 2 G. 
k k 
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It is possible to choose, given the various 1^ (k) and the duration D(k), those 
values for [HD]^ which minimize the above noise term; subject to the constraint 
that M = q.^, the target quantity, regardless of what values the other quanti- 
ties may have (to the extent that the assumed linear representation is valid) . 
The solution, derived elsewhere (ref. 35) is: 

[HD], = E at! I (k)/G(k), 
j J » J ^ J 

where A"^ is the inverse of the matrix 

A = E I (k) I (k)/G(k). 

J » k J 


(Knowing D(k), H(k) = HDj^/D(k).) 

A couple of points should be noted. First, the weights produced in this 
manner may be multiplied by an arbitrary scale factor, to give a new set of 
weights; however, the signal to noise ratio will remain unchanged. Secondly, 
if G(k), (the squared noise at a given point) is multiplied by a scale factor; 
(i.e., if all noise terms a^ are increased by the same factor), then there is 

no change to the weights produced. In fact, G(k) may be written as: 

^ G(k) = const X [1 + R 

where R = [a 2 /D(k) + a 4 ]/[ai/D(k) + as] 

represents the relative importance of the multiplicative terms with respect to 
the additive terms. If we assume D(k) is constant, then this one parameter will 
determine the weights produced; whether they optimize against multiplicative 
noises (R>>1), additive noises (R<<1), or some intermediate mixture. Analysis 
of the way R varies with the a's for reasonable values of the noises should in- 
dicate which terms will have most effect on the shape of the weights produced 
when varied. 


Similarly, it is possible to examine the relative importance of the two 
multiplicative terms (in the numerator). In fact, the ratio of their contribu- 
tion to the total final noise is given by a 4 x D/a 2 . The same is true for the 
two additive terms. On the other hand, to examine the total contribution of 
both multiplicative terms in comparison to the contribution of the additive 
terms, it is necessary to first calculate the weights. Once this has been done 
the absolute importance of the four noises may then be calculated, in the same 
units as those in which the target burden is measured: 


N? = 

ai 

E 

(HD)2/D 

Nf = 

a2 

E 

(HD)2 

n2 = 

as 

E 

(HD <I>)^/D 

< = 

a4 

E 

(HD <I>)2 
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4.2.4 Physical Origins for the Noise Terms .- The various noise terms may 
each arise due to one or more sources. Some possible sources are listed belov/: 


RA: - Detector noise. 

- Photon noise. 

- Other electrical noises. 

Digitization of integrated outputs. 

RM: - Scintillation or rapid variations in target, albedo, 

or illumination. 

Random errors in derivation or reference signal. 
Random variations in scan waveform. 


SA; 


SM: 


Presence of spectral signatures which were not repre- 
sented in sample used to derive weights. 

Crosstalk from reference into signal channel. 

Inaccuracies in carrying out multiplications 
(truncation of weights) . 

Change in scan waveform from the shape when the weights 
were derived. 
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5. BREADBOARD MODEL CORRELATION INTERFEROMETER 


Within the objective of the work on the Carbon Monoxide Pollution Experi- 
ment (COPE) program, a major goal was the development of the correlation inter- 
ferometer instrument for quantitative measurement of carbon monoxide, and other 
trace constituents in the atmosphere. The details of the design, construction, 
and testing of the breadboard correlation interferometer are given in a previ- 
ous report (ref. 2) . 

The specific objectives of this phase of the program were, by building and 
testing a breadboard interferometer, to establish the feasibility of the method, 
to gain operational experience and experimental data on the interferograms of 
the gases of interest and to obtain inputs for the design of an engineering mod- 
el. As these objectives imply, the breadboard has provision for independently 
varying several of its parameters, so that the effects of these parameters could 
be investigated. This means that the breadboard instrument is more complex than 
later versions, in which these parameters will have optimum, fixed, values se- 
lected on the basis of experience with the breadboard. One such parameter was 
the wavelength. The breadboard was built for operation at both 2.35 ym (first 
overtone of CO) and at 4.6 ym (fundamental of CO). Based on the results ob- 
tained and on analytical results previously presented, the 2.35 ym band was cho- 
sen for the engineering model . The various parameters of the breadboard model 
are given in Table 5.1. 

The laboratory tests should cover the range from 1 to something less than 
0.1 atm-cm, a reasonable lower limit being 0.02 atm-cm, the limits being imposed 
by the mapping experiment. The ground level sinks considered cause changes of 
optical thickness which require that a decrease of the order of 10% be detect- 
able. That is, the absolute limit of sensitivity should be such that a column 
density of 0.02 atm-cm (1/10 that of the nominal unpolluted atmosphere) is mea- 
surable. 

The experiments were planned to test the ability of the breadboard model 
correlation interferometer to measure carbon monoxide in its optical line of 
sight to these levels. Since the capability of making these CO measurements 
with sufficient accuracy in the presence of atmospheric constituent gases which 
interfere spectrally with CO is the primary question to be answered, the tests 
set up realistic combinations of carbon monoxide and interferent gases over as 
wide a range as conditions as is practical . 

The first series of tests were designed to verify the instrument operation 
in the laboratory using the General Electric multiple chamber system. In this 
facility the carbon monoxide and selected interferent gases can be introduced 
into the optical line of sight of the instrument in known quantities, and the 
effects of gas temperature and pressure on the measurements can be studied. 

The limitation of these tests is that not all of the atmospheric constituent 
gases can be introduced into the chambers in sufficient quantities to approxi- 
mate atmospheric amounts. Water vapor is a major interferent whose effect can 
not be simulated in these laboratory tests. However, the laboratory tests pro- 
vided a first check on the instrument operation at the minimum CO level which 
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TABLE 5.1 BREADBOARD SYSTEM PARAMETERS 


Aperture 

Interferometer 

Telescope 

Field of View 
Interferometer 
Telescope 

Spectral Bands 

. Delay Scan Range 

Number of Sample Points 
Sample Length 
Scan Rate 

Number of Scans Accumulated 

Noise Equivalent Power 

Noise Equivalent CO Amount 
(2% Albedo, x = 1 sec) 

Detector 2.3 pm 

4.6 pm 

Weight 

Interferometer 

Telescope 

Electronics 

Power 


6.6 cm diameter 
22.0 cm diameter 


0.12 radian diameter 
0.034 radian diameter 

4240 - 4340 cm'^ (2000 - 2200 cm- 

2.5 - 4.0 mm 

0- 64 

1- 63 fringes 
1 Hz 

1 - 500 

1.6 X 10-^^ watt/Hz^/^ 

0.004 atm- cm 

LN 2 cooled PbS immersed on SrTi 03 
LN 2 cooled InSb 


15 lb. 

10 lb. 

60 lb. 

100 watts 



one wishes to detect. Tlie laboratory tests were carried out using the 2.3 um 
overtone band and using the 4.6 urn fundamental band. 

The second phase of the tests was a series of outdoor measurements with 
the instrument looking through an actual atmospheric path at sunlight. These 
outdoor tests indicated the sensitivity of the correlation instrument measure- 
ment of CO over a range of atmospheric parameters which were covered in this 
test period. However, these tests did not permit the instrument response to be 
evaluated over the full range of atmospheric parameters to be encountered on a 
long term seasonal and latitudinal basis. The outdoor tests were carried out 
using the 2.3 ym overtone band. 

The third phase of the tests was a series of measurements with the instru- 
ment looking through a 430-foot tunnel at NASA/LaRC. The composition of the gas 
in the tunnel, including an appreciable amount of water, was controlled. 

The testing of the correlation interferometer showed that this technique 
is capable of making the desired measurements over the important range of CO in 
the presence of atmospheric gases. 

It was found the atmospheric amounts of CO can be measured with an error 
of 10% or better. The instrument test confirmed the conclusion reached ana- 
lytically that the use of the overtone band at 2.3 ym gives much better results 
than does that of the fundamental band at 4.6 ym. The range of sensitivity at 
2.3 ym is much larger, the absorption being nearly linear from 0 to above 1 atm- 
cm and usable to at least 10 atm-cm; whereas at 4.6 ym near linearity is very 
limited, to a factor of about five in the range of densities of interest. The 
effects of gas temperatures and pressures and source temperatures are negligible 
using the 2.3 ym band but quite significant using the 4.6 ym band. 

Carbon monoxide can be measured accurately in the presence of the major 
interferents for this spectral region - methane and water. The methane effect 
has been found to be completely separable from that of CO. These same data can 
be used to measure methane as well as CO with errors of less than 10% over the 
range of 1 to 5 atm-cm, which covers that to be expected in atmospheric measure- 
ments. The tunnel tests showed that, for water densities up to 230 atm-cm, the 
effect on CO accuracy is separable and minimal. The outdoor tests showed that 
the effect of atmospheric amounts of water can be overcome by a proper choice 
of correlation functions. The choice of conditions used to determine correla- 
tion functions is critical. The conditions should cover the range of variation 
of interferents and atmospheric parameters which affect the measurement. The 
use of a wide range of conditions in obtaining correlation functions and inter- 
polation between these gives much better results than use of a narrow range 
with extrapolation. The outdoor tests showed that using a broad range of ap- 
propriate conditions for these functions greatly improves results. 

The breadboard experimental results gave good and necessary data for ob- 
taining design parameters for the engineering model . 
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6. 


THE DESIGN, FABRICATION, AND GROUND-BASED TESTING OF 
THE ENGINEERING MODEL CORRELATION INTERFEROMETER 


6.1 System Design 

The design of the engineering model was based "on that of the breadboard 
instrument, modified as indicated by the results of breadboard testing or as 
required by the more extended objectives of the engineering model program. 

The functions to be performed by the system, subject to the constraints 
of aircraft and balloon operations, are as follows: 

1. To collect radiation, in sufficient quantity, from a field of 
view appropriate to the intended flight test program. 

2. To limit this radiation to the narrow spectral band where absorp- 
tion by CO occurs. 

3. To form the interferogram of the radiation as a function of time 
and convert it to a time-varying voltage for signal processing. 

Since the radiation is limited to an optical band of ± Ao the 
voltage is band- limited to frequencies So^ ± SAo where S is rate 
of scanning path difference. 

4. To normalize the signal amplitude in terms of the total received 
power in the spectral band. This is required because the observ- 
able which is related to the CO amount is the transmission of the 
atmosphere. This normalization should be capable of dealing with 
intensity changes over the entire information frequency bandwidth 
of the instrument (i.e., 0 to SAo) . 

5. To move signal frequency band (from So^ ± SAa to 0 ± Aa) by syn- 
chronous detection and smoothing. The use of synchronous detec- 
tion (or hetrodyning) gives better S/N ratio for small signals 
than linear or square law detection (refs. 36,37) and does not 
require that the signal frequency be maintained exactly at the 
peak of a narrow band electrical filter. Transforming to zero 
frequency minimizes the position accuracy required on subsequent 
sampling. 

6. To sample the smoothed signal, at known points in the interfero- 
gram, over a range of delay centered around the CO signature. 

The spacing of the points to be that dictated by the bandwidth 
and the sampling theorem, i.e., at least one point every 1/2 AAo 
units of delay. 

7. To digitize the value of each sample and output this data in a 
form suited to telemetry and/or magnetic recording. 

8. Measure and output, in suitable format, the housekeeping data 
required during flight testing. 

The means by which these functions are implemented are shown in Figure 6.1 and 
described below in general terms and in detail elsewhere (ref. 3) . 
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Figure 6.1 System Block Diagram 





























The incident radiation is collected by 4.8 inch diameter, £/4, Maxwellian 
fore-optics with a 7.5° field stop. The radiation then traverses the field 
widened, refractively scanned interferometer as a 7° beam and is then imaged 
onto the immersed, dry-ice-cooled, lead sulphide detector by a f/1.5 silicon 
lens. The spectral band is defined by a temperature controlled interference 
filter located at the detector housing and a rotary- solenoid-actuated shutter, 
directly in front of the filter, serves to interrupt the radiation, when re- 
quired, for automatic gain control (AAGC) . 

The path difference in the interferometer is scanned by tilting the scan 
plate by means of a torque motor, using a drive waveform shaped so that the 
path difference decreases linearly over 3/4 of the scan and the plate returns 
rapidly to zero during the remaining 1/4. During the, linear portion of the scan 
the main detector sees a portion of the radiation as an amplitude modulated 1 
kHz sinusoid. This is the interferogram of the absorption spectrum due to the 
atmospheric path traversed by the radiation. 

A second optical path, co-axial with the main beam, takes radiation from 
a quartz/iodine lamp through a collimator, a spectral filter and the interfer- 
ometer, exiting, just in front of the main silicon lens, to a gas cell, a tele- 
scope and a room temperature PbS detector. The output of this detector is an- 
other 1 kHz sinusoid, the reference interferogram, whose shape is determined by 
the absorption in the gas cell. The reference interferogram establishes the 
path difference scale for the scan. 

The signal from the main detector preamplifier goes to the dividend input 
of an analog divider, for which the divisor input is a voltage proportional to 
the irradiance level at the instrument. This signal, which monitors the total 
irradiance and its fluctuations up to 20 Hz, is generated by a bore-sighted, 
radiometric (fast AGC) channel, operating in the same spectral band as the in- 
terferometer. 

The quotient output of the divider is then proportional to the fraction 
of the received radiation which is modulated by the interferometer. To main- 
tain the proportionality constant, in spite of long-term changes in the rela- 
tive responsivity of the main and FAGC channels, the radiation to the main de- 
tector is, periodically, interrupted and the amplitude of the resulting pulse 
is measured, at the quotient output of the divider, digitized (as the recipro- 
cal), displayed and recorded. If the relative responsivity of the main and 
FAGC channels remains constant, so does the pulse, if not, the pulse amplitude 
is proportional to the relative responsivity and is used, during data reduc- 
tion, to correct for any relative responsivity changes. 

From the FAGC divider output, the main signal is filtered by 10% bandwidth 
filter centered at 1 kHz and then input to two synchronous rectifiers, which are 
just electronic switches, opened and closed by square waves derived from the 
peaks and zero crossings, respectively, of the reference interferogram. Two 
channels are used, separated by 90° in phase, because the main interferogram is, 
in general, out of phase with the reference interferogram. 

The outputs of the synchronous rectifiers are supplied to two dumping inte- 
grators whose outputs are sampled and reset to zero at fixed increments of path 


83 



difference as determined by counting the zeros of the reference interferogram. 
The integrator outputs are sampled, sequentially, by a 12-bit A/D converter. 

The digitized values (2 per sample) are buffered and clocked out to telemetry 
while the next sample is being taken. 

The portion of interferogram sampled is delimited by the "RUN” pulse, 
which goes high when the scan plate passes through a given position, as deter- 
mined by a LED/photodiode sensor, and goes low after n(m + 1) - 1 zero crossings 
of the reference interferogram, where m is the length of a sample and n is the 
number of samples/scan. The integrators are allowed to run for m cycles of the 
reference, the remaining cycle being used to dump the accumulated charge. As a 
check on scan operation, the length of the run pulse is measured in terms of an 
internal crystal clock, displayed on the LED display and is also output. 

The temperature of the main channel spectral filter and of the enclosure 
containing the interferometric components and the reference channel filter, are 
both controlled above ambient by proportional and integral controllers and 
thermistor sensor. Even though the path difference in the interferometer is 
stabilized by the use of an ultra low expansion structure and by temperature 
control, so that the delay is essentially controlled by the position of the scan 
plate, it is also advisable to have a means of monitoring the absolute delay, so 
as to check for any long term drifts or displacements due to mechanical shock. 
The number of reference interferogram peaks occurring between the mechanically 
generated start of scan pulse and a pulse denoting the position of a local maxi- 
mum in the reference interferogram is counted and displayed. 

For test and checkout purposes, the system can be interfaced with a Nova 
mini-computer, so that the resulting system very closely resembles the bread- 
board model of the correlation interferometer. 

As well as storing and accumulating the sample values, the Nova can be 
used to control the engineering model system, as was done in the breadboard in- 
strument. The control program is conversational, via a teletype, and permits 
the operator to set the number of samples/scan and the number of scans to be 
co-added and to have the accumulated sample values either printed out and/or 
punched on paper tape or to have them multiplied by any one of a set of 8 
stored correlation functions and output the sum of the resulting products. The 
stored correlation functions are input, as required, from paper tape and there 
is also the facility to re-enter a previous set of sample values from tape and 
treat it with correlation functions. All operator input points in the program 
are trapped to prevent accidental overwriting of stored program or data and the 
A/D converter data is trapped for overflow. 

Comparing the engineering model with the breadboard model of the correla- 
tion interferometer, the major changes in functional blocks are the addition of 
control, telemetry and "housekeeping" blocks to replace the Nova mini -computer, 
the provision of "fast" AGC, the substitution of unit-magnification Maxwellian 
fore-optics for the Newtonian telescope used in the breadboard and the separa- 
tion of the main and reference channels so that they no longer share a single 
field stop and spectral filter. 
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6 . 2 System Parameters 

The data and operating experience gained during the construction and test- 
ing of the breadboard model correlation interferometer (ref. 2) showed that, 
while there were several areas in which the design could be improved, the major 
radiometric and interferometric parameters were well chosen. Accordingly, these 
parameters (etendue, detector NEP, operating wavelength and scan range) were 
carried over to the engineering model design. The external field of view was 
changed from 2 x 2° to 7 x 7° (equal to the internal field of view of the inter- 
ferometer) to give a ground resolution element compatible with the V/h charac- 
teristics of the intended flight vehicles. This change in field permitted a 
reduction in instrument size by the elimination of the fore-optics telescope 
used with the breadboard. 

The Engineering Model Correlation Interferometer is shown in Figure 6.2. 

The light-colored object in the left center is the frame made of ultra-low-ex- 
pansion material. The gray object seen through the frame opening is the oscil- 
lating compensating plate. The two fixed mirrors are shown on the extreme left 
and in the center. The detector is on the upper right. 

The major system parameters for the engineering model are listed in Table 
6.1. The bulk of the changes from the breadboard design occurred at the sub- 
system level. 


6,3 Engineering Model Correlation Interferometer Testing 

In order to properly evaluate the performance of the engineering model 
correlation interferometer a series of tests of the instrument was carried out. 
These consisted of tests in the laboratory, outdoors at Valley Forge, Pa. and 
Hampton, Va., on board a Falcon Fan-Jet and on board a Jet Ranger Helicopter. 

The flight tests are described in sections 7 and 8 of this report. 

6.3.1 Multiple Chamber Test Facility.- The multiple chamber test facility 
was designed to permit known amounts of test gas to be introduced into the opti- 
cal path of the instrument between the source and the fore-optics. Since it was 
necessary for the complete 7° field of view of the instrument to be filled by an 
extended area source of nearly uniform radiance, the cells were constructed with 
an inside diameter of 30.5 cm (12.0 in.). This permitted cells with a total 
length of 150 cm (59.05 in,) to be connected together, and still allowed the 
conical field of view of the instrument to clear the inside surfaces of the cells 
and view the source. The details of the test facility are shown in Figures 6.3 
and 6.4. 

Each of the three cells which comprise the multiple chamber system is 
equipped with a vacuum gauge which is used for testing for leaks and for deter- 
mining when the cell is sufficiently evacuated to be filled with test gas. The 
cells can be isolated from a common manifold which is used for evacuating the 
cells and for filling them with the proper gas mixture. The pressure of the gas 
in the pumping mainfold, and hence the pressure in the individual cells with the 
appropriate valves opened, is determined by means of a mercury manometer. 
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Figure 6.2 Engineering Model Correlation Interferometer Sensor 



TABLE 6.1 ENGINEERING MODEL SYSTEM PARAMETERS 


Aperture 

Interferometer 

Field of View 

Interferometer 

Spectral Band 

Delay Scan Range 

Number of Sample Points 
Sample Length 
Scan Rate 

Noise Equivalent Power 

Noise Equivalent CO Amount 
(2% Albedo, t = 1 sec.) 

Detector 

Weight 

Interferometer 

Electronics 

Power 


6.6 cm diameter 

0.12 radian diameter 
4260 - 4400 cm-1 

2.5 - 3.9 mm 

0- 64 

1- 63 fringes 
1 Hz 

1.6 X 10-^^ watts/Hz^/^ 

0.004 atm- cm 

Cooled C195 °K) PbS 
Immersed on SrTiOs 

55 lb. 

41 lb. 

168 watts at 24 volts 





GAS BOTTLES WITH 



Figure 6.4 Multiple Test Facility 






The individual cells are constructed from aluminum pipe with an inside 
diameter of 30.5 cm (12.0 in.). Aluminum flanges are welded to each end of the 
pipe to make each cell a self-contained unit which can be used individually or 
bolted together with other cells to form a longer optical path. The windows, 
which are 2.54 cm (1.0 in.) thick polished pyrex plate glass, are clamped be- 
tween the cells and sealed with 0-rings to provide a leak-free system. Each 
cell forms a unit with an optical path length of 50 cm (19.68 in.). 

A tungsten halogen lamp provides the infrared energy source for the tests. 
The lamp illuminates a set of three sheets of pyrex plate glass with a rough 
ground surface. These scatter and diffuse the energy to provide a uniform dif- 
fuse source for the interferometer to look at. A variable voltage transformer 
is used to control the lamp power and hence the light output. For outdoor test- 
ing the sun's rays are directed parallel to the axis of the cell system onto 
the diffuser plates by means of plane or concave mirrors. 

The cells, the source, the pumping manifold, and the correlation interfer- 
ometer were all mounted on a laboratory table provided with wheels. This makes 
the complete system mobile and quite easily moved outdoors for tests. 

6.3.2 Laboratory Tests of the Correlation Interferometer.- The laboratory 
tests of the instrument were designed to determine the sensitivity of the corre- 
lation interferometer to carbon monoxide in the presence of interferents when 
the measurement light levels were equivalent to an earth albedo of 30 to 40%. 

The laboratory tests of the breadboard (ref. 2) demonstrated the ability of the 
interferometer to make the type of carbon monoxide measurements required. The 
lead sulphide detector initially installed in the engineering model correlation 
interferometer is the same detector which had been used in the breadboard in- 
strument. It became necessary to use the same detector when the detector manu- 
facturer was unable to supply a new detector with higher sensitivity in time for 
the instrument tests. 

All of the interferent gases cannot be introduced into the test cells in 
the quantities present in the atmosphere. Water, which is a significant inter- 
ferent at 2.35 ym cannot be introduced in significant quantities in the 150 cm 
(59.05 in.) cell length available. However, the effects of methane and carbon 
monoxide can be studied with realistic atmospheric amounts. A study of the cor- 
relation interferometer response to combinations of these two gases is suffi- 
cient to determine its operating parameters. 

The gases used for these tests were various premixed gases (30% methane in 
nitrogen, 5.21% methane in nitrogen, 2.1% carbon monoxide in nitrogen, 1% carbon 
monoxide in nitrogen, and 0.1% carbon monoxide in nitrogen) together with pure 
nitrogen for use in diluting the available mixtures when necessary. The optical 
thickness of gas desired is introduced into the system by (1) evacuating the 
system, (2) filling one or more cells with the required mixture, (3) re-evacu- 
ating the remaining system, and (4) repeating the procedure until all needed 
cells are filled. 

The nominal background amounts of carbon monoxide and methane in a double 
vertical path through the atmosphere are given as: 
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1. 0.18 atm-cm of carbon monoxide (corresponding to a constant 
mixing ratio of carbon monoxide of 0.1 ppm) 

2. 2.80 atm-cm of methane (corresponding to a constant mixing 
ratio of methane of 1.6 ppm) 

These gas quantities have been computed from atmospheric concentration mod 
els for carbon monoxide and methane which have been discussed previously (ref. 
5). 


The laboratory tests were conducted with a small matrix of methane and car 
bon monoxide conditions. The methane consisted of 3 amounts from 1.75 to 3.50 
atm-cm and the carbon monoxide consisted of up to 5 amounts ranging from 0 to 
1.0 atm-cm. The test matrix with a range of carbon monoxide and methane mix- 
tures is necessary to permit computation of a correlation function applicable 
over the range of test gases selected. The basic procedure in generating a set 
of data to be used in determining a suitable correlation function for the meth- 
ane and carbon monoxide mixture, is the following: 

1. A test condition is set up in the gas cells which corre- 
sponds to a nominal background amount of the gases used. 

In this case 2.6 atm-cm of methane and 0.2 atm-cm of car- 
bon monoxide was used since this corresponded to approxi- 
mately the amount of gas in a vertical atmospheric column 
for a double pass through the atmosphere. 

2. A series of interferogram scans is then obtained and co- 
added at some predetermined sampling interval. In this 
case interferograms each consisting of 31 sampling points 
equally spaced over the selected path difference interval 
were collected. A minicomputer data handling system col- 
lects and formats the data, and outputs a punched paper 
tape record. 

3. The nominal amount of carbon monoxide in the cell is then 
increased by a predetermined amount, known as the target 
amount, which was 0.30 atm-cm for the purposes of the tests. 

The interferogram for this methane and carbon monoxide con- 
dition is then sampled as above. 

4. This is repeated for a series of test conditions in which 
the carbon monoxide is maintained at 0.2 atm-cm and the 
methane is varied. For these experiments the three meth- 
ane amounts varied from 1.75 to 3.5 atm-cm. 

The measurements were run twice to provide duplicate measurements to be 
used as test cases. In addition, some other combinations of methane and carbon 
monoxide were run as test cases . 

The data is all recorded on punched paper tape and then is input through 
the GE time-share computer system for storage on disc. Using these stored in- 
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terferograms , various correlation functions can be computed and applied to in- 
terferograms of test cases to obtain a measure of the carbon monoxide in the 
path. 


The procedure for computing the correlation functions is described in Sec- 
tion 4.2. The application of the correlation function to a measured interfero- 
gram to yield a corresponding value of CO or CH 4 is as follows: 


1. The interferogram consists of a set of 31 in-phase points, 
[I(k)], and 31 quad-phase points, [Q(k)], due to the syn- 
chronous detection process. 


2. The correlation function similarly has values H[I(k)] and 
H[Qfk)] corresponding to each of the interferogram points. 

3. The measured value of CO, for example, is given by the ex- 
pression 

N 


q(CO) = V |H[I(k)]I(k) + H[Q(k)]Q(k)l x^ 

{ target ' ^^^^^nominal} ^ nominal 


k=l 


where : 
q(CO) 
q(CO) 


target 


is the measured value of CO 

is the value of CO in the interferogram used 
as a target increment in the correlation 
function matrix 


qfCOl , is the value of CO in the interferogram used 

•^nominal , , , 4 . ...u ? 4 .- 

as a background amount in the correlation 

matrix. 


N 

c 


is the number of scans co-added to form the 
matrix of correlation interferograms 


N 


m 


is the number of scans co-added in the in- 
terferogram to be measured. 


In a similar manner it is possible to set up a matrix of test cases to 
measure methane (or other gases) and compute a correlation function. 


Following the September 1973 flights in Ottawa, Canada, the interferometer 
was found to have had a mechanical displacement in the torque motor which drives 
the oscillating refractor plate. This was corrected and some modifications were 
also made in the detector pre-amplifier circuit and in the temperature control 
circuits which stabilize the interferometer temperatures. When these modifica- 
tions were completed the new lead sulphide detectors were available and one of 
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these was installed before further laboratory testing was initiated. The new 
detectors (designated #001 and #992) were found to have a nitrogen hold time of 
about 6 hours at room temperature, and when installed in the interferometer the 
hold time with dry ice and methanol was better than 12 hours. Methanol or ace- 
tone was used in the detector in place of the glycol when it was found that the 
glycol became quite viscous at dry ice temperatures raising the possibility that 
the detector temperature might not be constant. 

With the new detector installed, matrices of test cases were run to again 
verify instrumental noise and offsets from sources other than detector noise. 

The new detector and its spare were replaced and interchanged to determine the 
magnitude of inaccuracies which might result if replacement of a detector became 
necessary. In addition, a series of measurements was conducted with the inci- 
dent light level chopped at frequencies between 5 and 60 Hz to attempt to simu- 
late the effects of varying ground albedo on the measurements. These data are 
given later in the section describing the results of laboratory tests. 

6.3.3 Outdoor Correlation Interferometer Tests.- After the operation of 
the engineering model correlation interferometer had been verified in the labo- 
ratory, it was taken outdoors to make measurements using radiation which has 
traversed a real atmospheric path. The outdoor tests were designed to demon- 
strate the ability of the instrument in making measurements of carbon monoxide 
and methane in the presence of large amounts of water vapor, which is an inter- 
ferent at 2.35 pim. 

The primary purpose of the outdoor tests was to obtain a set of interfer- 
ograras from the correlation interferometer covering a wide range of atmospheric 
conditions. Some of these interferograms were then used to generate correlation 
functions which were optimized for detecting carbon monoxide and methane against 
their background atmospheric inter ferent gases over the range of test conditions 
for which the data were generated. 

The multiple chamber system used during the laboratory tests was taken out- 
doors for use in adding increments of known amounts of carbon monoxide and meth- 
ane to the optical path which included a single atmospheric path. For these 
tests the light source is the sun itself. A 35.6 cm (14.0 in.) square mirror 
was used to direct the sunlight into the cells approximately parallel to the ax- 
is of the cells and the instrument. The sunlight was diffused at the entrance 
window of the cells by directing the sunlight onto a ground pyrex diffuser plate 
which covered the entire cell aperture. A total of three diffuser plates was 
employed to attenuate the sunlight at the interferometer to a value which corre- 
sponds to that seen by the instrument looking at a sunlit body with an albedo 
in the range of 30 to 40%. 

A number of outdoor measurements were made with the apparatus shown in Fig- 
ure 6.5 using sunlight in the vicinity of the General Electric Valley Forge 
Space Center. The solar evaluation angles during the test were such that the 
atmospheric air mass covered varied by a factor of about 3, with a corresponding 
variation in the range of carbon monoxide and interferent gases. 

The test procedure was as follows: 
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Figure 6.5 Outdoor Experimental Setup 


1. An inteferogram was obtained with the cells evacuated 
and only the background atmospheric amounts of carbon 
monoxide and interferents in the optical path. 

2. The carbon monoxide in the path was then successfully 
incremented by 0.5, 1.0, and 1.5 atm-cm by introducing 
the bottled carbon monoxide mixture into the cells, 
and the corresponding interferograms obtained. 

3. The carbon monoxide was pumped out of the cells and the 
methane was then incremented by 1.5, 3.0, or 4.5 atm-cm 
and the interferograms obtained. 

4. Other combinations of carbon monoxide and methane in- 
crements were also run, and interferograms for the at- 
mospheric background amounts of gas were also obtained 
at regular intervals. 

The purpose of collecting a large number of interferograms during the 
course of a day was to obtain a matrix of test cases which covered the range of 
carbon monoxide, methane, and water anticipated during flight testing of the 
instrument. The outdoor tests provided a set of measurements which was used in 
computing correlation functions for carbon monoxide and methane, which were then 
applied to the outdoor data (not used in generating the correlation functions) 
as a test of the instrument's measurement capabilities. 

6.3.4 Corroborative Measurements.- At the same time that outdoor measure- 
ments were being made with the correlation interferometer, auxiliary solar spec- 
tral absorption measurements were being made in order to define the atmospheric 
amounts of carbon monoxide and methane at the time of the measurements. This 
measurements program was carried out by NASA Langley Research Center personnel 
operating from their instrumented van located at Valley Forge, Pa. The instru- 
mentation within the van consisted primarily of a Fourier transform spectrometer 
with a maximum spectral resolution of 0.1 cm”^. A sun tracker mounted on the 
roof of the van permitted the sun's rays to be continuously directed into the 
spectrometer to obtain atmospheric absorption spectra. The solar absorption 
spectrum covered the carbon monoxide fundamental band at 4.6 ym and the tail of 
the methane band in the region of 3.7 ym. These spectra were obtained at a res- 
olution of 0.25 cm-^ using a liquid nitrogen cooled HgCdTe detector. Although 
the interferometer spectrometer has a maximum resolution of 0.1 cm“^, the spec- 
tra were obtained at a resolution of 0.25 cm-^ in order to provide spectra with 
acceptable signal to noise in a reasonable measurement time. Since the sun 
tracker did not provide perfect tracking it was necessary to compromise on the 
spectral resolution. The measurement time for the 0.25 cm"^ spectra was up to 
2 minutes. This resolution was found to be acceptable for the purpose of deter- 
mining the atmospheric levels of carbon monoxide and methane. The atmospheric 
spectra were obtained at periodic intervals during the day while the measure- 
ments with the correlation interferometer were being made. 

In order to relate the spectra with specific amounts of carbon monoxide 
and methane, calibration absorption spectra of known amounts of carbon monoxide 
and methane in a cell were obtained. These were used to generate calibration 
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curves of line absorption vs gas amount for the carbon monoxide and methane 
lines which were least overlapped by water vapor lines in the atmospheric spec- 
tra. These calibrations were then used to determine the amounts of carbon mon- 
oxide and methane in the atmospheric absorption spectra taken during the time 
that atmospheric measurements were being made with the correlation interferom- 
eter. 


Following the Canadian flights in April 1974, the correlation interferom- 
eter was taken to NASA Langley Research Center for calibration and determination 
of weighting functions to be applied to the flight data. The calibration was 
carried out using the multiple chamber test facility in a manner similar to that 
described in Section 6.3.3. The instrumented NASA van was also employed at this 
time to provide the spectral data needed to determine the CO and CH 4 in the at- 
mospheric path at the same time that the calibration measurements were being 
made with the correlation interferometer 


6.4 Engineering Model Test Results 

6.4.1 Laboratory Tests.- The laboratory tests conducted with the engi- 
neering model were mainly designed to determine the instrumental measurement 
noise and stability. Since the lead sulphide detector used in the instrument 
during these initial tests was that from the breadboard interferometer, the 
noise and sensitivity characteristics were not as good as calculated, based on 
the specifications of the detectors on order. 

The laboratory data in Table 6.2 were computed from the test matrix with a 
nominal carbon monoxide value of 0.197 atm-cm and a target carbon monoxide val- 
ue of 0.55 atm-cm. In these tests, conducted in May 1973, the nominal case was 
chosen to have carbon monoxide and methane values close to their expected val- 
ues during outdoor measurements. The target case then increments only the car- 
bon monoxide by a pre-selected amount. In this way the instrument response to 
a known amount of carbon monoxide was determined. In this case the test data 
were computed on the basis of 128 co-added scans. Table 6.3 indicates that the 

measurement accuracy is generally much better than 10%. Extrapolation of the 

carbon monoxide values to 1.00 atm-cm yields a lower measured value of the car- 
bon monoxide by 5% or less. However, the breadboard tests indicated that the 
carbon monoxide absorption is becoming non-linear for these larger amounts of 
carbon monoxide and hence the lower measured value of carbon monoxide is ex- 
pected. Extrapolation of the carbon monoxide to low or zero values results in 
a slight offset to the data which amounts to about -0.021 atm-cm. The rms mea- 
surement noise is of the order of 0.005 to 0.008 atm-cm. 

Table 6.3 shows the same data matrix computed on the basis of 8 co-added 
scans. This data shows the same trends as the data discussed in Table 6.2 ex- 
cept that the measurement noise level appears to be higher in general . The rms 
measurement noise is 0.015 atm-cm of carbon monoxide. 

In March and April of 1974 a series of laboratory tests of the interferom- 

eter were conducted with the new lead sulphide detectors installed in the in- 
striiment. A matrix of test cases similar to those obtained for the breadboard 
detector gave the data in Tables 6.4 and 6.5. Figure 6.6 is a plot of the 31 
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TABLE 6.2 


CH 4 

atm-cm 


1.75 

1.75 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

2.60 

3.50 

3.50 


LABORATORY CARBON MONOXIDE MEASUREMENTS 
MAY 1973 


NOMINAL CO =0.197 atm-cm 
TARGET CO = 0.500 atm-cm 
NOMINAL CH 4 = 2.60 atm-cm 
OTHER CH 4 = 1.75, 3.50 atm-cm 
NEASUREMENT TIME =128 Scans 
BREADBOARD DETECTOR (193 K) 


CO 

atm-cm 

MEASURED 

CO 

atm-cm 

ACCURACY 

% 

.197 

.211 

7.1 

.197 

.211 

7.1 

0 

- .030 

— 

0 

- .010 

— 

0 

- .025 

— 

.025 

.005 

- 80.0 

.025 

- .006 

124.0 

.025 

- .005 

120.0 

.197 

.178 

- 9.6 

.197 

.186 

- 5.6 

.500 

.491 

- 1.8 

.500 

.493 

- 1.4 

1.000 

.950 

- 5.0 

1.000 

.959 

- 4.1 

1.000 

.971 

- 2.9 

.197 

.192 

- 2.5 

.197 

.209 

6.1 
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TABLE 6.3 LABORATORY CARBON ^ft)NOXIDE MEASUREMENTS 

MAY 1973 


NOMINAL CO = 0.197 atm-cm 
TARGET CO = 0.500 atm-cm 
NOMINAL CH 4 = 2.60 atm-cm 
OTHER CH 4 = 1.75, 3.50 atm-cm 
MEASURENENT TIME =128 Scans 
BREADBOARD DETECTOR (193 K) 


CH 4 

atm-cm 

CO 

atm-cm 

MEASURED 

CO 

atm-cm 

ACCURACY 

% 

1.75 

.197 

.209 

6.1 

1.75 

.197 

.196 

- 0.5 

1.75 

.197 

.211 

7.1 

2.60 

0 

- .042 

— 

2.60 

0 

- .006 

, 

2.60 

0 

- .021 

— 

2.60 

.025 

.027 

8.0 

2.60 

.025 

- .015 

160.0 

2.60 

.025 

.008 

- 68.0 

2.60 

.197 

.168 

- 14.7 

2.60 

.197 

.192 

- 2.5 

2.60 

.197 

.211 

7.1 

2.60 

.500 

.553 

10.6 

2.60 

.500 

.490 

- 2.0 

2.60 

.500 

.494 

- 1.2 

2.60 

1.000 

.915 

- 8.5 

2.60 

1.000 

.974 

- 2.6 

3.50 

.197 

.185 

- 6.1 

3.50 

.197 

.158 

- 19.8 

3.50 

.197 

.193 

- 2.0 
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TABLE 6.4 LABORATORY CARBON MONOXIDE MEASUREMENTS 

APRIL 1974 


NOMINAL CO = 0.208 atm- 

cm NOMINAL CH^ 

= 2.60 atm-cm 

TARGET CO 

= 0.525 atm- 

cm OTHER CH 4 = 

1.74, 3.48 atm 

NUMBER OF 

SCANS = 128 

DETECTOR #002 (193 K) 

ACTUAL 

ACTUAL 

MEASURED 

ACCURACY 

CH 

CO 

CO 

ACCURACY 

atm- cm 

atm- cm 

atm-cm 

% 

1.74 

0.208 

.214 

2.9 

1.74 

0.208 

.214 

2.9 

2.60 

0 

.031 

— 

2.60 

0 

.019 

— 

2.60 

0 

.028 

— 

2.60 

0.208 

.209 

0.4 

2.60 

0.208 

.208 

0.0 

2.60 

0.525 

.480 

- 8.6 

2.60 

0.525 

.473 

- 9.9 

3.48 

0.208 

.185 

- 11.1 

3.48 

0.208 

.164 

- 21.2 

TABLE 6 . 

5 LABORATORY 

CARBON MONOXIDE 
APRIL 1974 

MEASUREMENTS 


NOMINAL CO 

= 0.208 atm-cm 

NOMINAL CH 4 = 

2.60 atm-cm 

TARGET CO 

= 0.525 atm-cm 

OTHER CH 4 = L 

.74, 3.48 atm 

NUMBER OF 

SCANS = 8 

DETECTOR #002 

(193 K) 

ACTUAL 

ACTUAL 

NEASURED 


CH 

CO 

CO 

ACCURACY 

atm-cm 

atm-cm 

atm-cm 

% 

1.74 

0.208 

.242 

16.4 

1.74 

0.208 

.219 

5.3 

1.74 

0.208 

.234 

12.5 

2.60 

0 

.015 

— 

2.60 

0 

.019 

— 

2.60 

0.208 

.206 

- 1.0 

2.60 

0.208 

.193 

- 7.2 

2.60 

0.208 

.178 

-14.4 

2.60 

0.525 

.438 

-16.6 

2.60 

0.525 

.480 

- 8.6 

2.60 

0.525 

.468 

-10.9 

3.48 

0.208 

.181 

-13.0 

3.48 

0.208 

.170 

-18.3 

3.48 

0.208 

.182 

-12.5 
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INTERFEROGRAM QUAD-PHASE AMPLITUDE INTERFEROGRAM IN-PHASE AMPLITUDE 


Figure 6.6 Interferogram In- and Quad-Phase Sample Values for 
128 Co-added Scans for Laboratory Measurements 
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interferogram sample points for the in-phase and quad-phase components of the 
interferogram co-added for 128 scans for two gas conditions in the cells. This 
actually represents 112 data scans due to the operation of the AAGC shutter. 

The interferogram output consists of 31 separate values for each phase. In Fig- 
ure 6.6 these points arc connected to show the general shape of the interfero- 
gram. However the values between the 31 sample points do not represent inter- 
ferogram values for a different sampling condition. The continuous interfero- 
gram curve shown would change for each sampling condition chosen. This is also 
true for the other interferograms plotted in Figure 6.8 and Figure 6.10. One 
case had 2.605 atm-cm of CH 4 and the other case had the same amount of CH 4 with 
an added increment of 0.525 atm-cm of CO. The second case is designated by x's 
in Figure 6.6 and indicates the change in the interferogram due to CO. Table 
6.6 indicates the weighting function values derived from the data matrix shown 
in Table 6.7 The weighting function consists of one value for each of the in- 
terferogram sample points. 

6 . 4. 1.1 Effect of removing or changing detectors: In the testing of the 

interferometer it was important to determine the effect on the measurement of 
removing a detector and then re-installing it again, or of replacing one detec- 
tor with another if one of them should fail. That is, determine if a calibra- 
tion were performed v\/ith a particular detector, would the measurements with the 
other detector, or the re-installed detector be acceptable. The measurement 
sequence followed was: 

1. To obtain a matrix of measurements with detector #002. 

2. Remove the detector and replace it with detector #001. 

3. Obtain another measurement matrix. 

4. Remove the detector and replace it with detector #002 again. 

5. Obtain the same measurement matrix again. 

6 . Generate a set of weights based on the last test matrix 

and apply them to the other measurements. 

These data are tabulated in Tables 6.4, 6 . 8 , and 6.9. 

Comparing the measurements obtained when the detectors are changed it is 
evident that no significant changes have occurred in the measurement interfer- 
ograms. The measurement accuracies are generally better than 10%. The few 
cases which are more than 10 % do not indicate a failure of the weights to apply 
but rather an infrequent transient noise in the system which causes a greater 
than normal error in the measurement. In Table 6.9 the zero carbon monoxide 
measurements show an offset of 0.050 atm-cm in the measured carbon monoxide 
value. However, the measurement of zero carbon monoxide in Table 6.4 (i.e., 
the matrix used in generating the weights) shows an offset of 0.026 atm-cm. 

Thus the real difference in the two sets of measurements is only 0.024 atm-cm 
of carbon monoxide. In actual measurements a correction could be found and 
applied. 
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TABLE 6.6 WEIGHTING FUNCTIONS APPLIED TO LABORATORY, 
OUTDOOR AND FLIGHT DATA 


CO LAB 

WEIGHT 

CO LaRC/LAB WEIGHT 

CH 4 LaRC WEIGHT 

IN-PHASE 

Q- PHASE 

IN-PHASE Q- PHASE 

IN-PHASE Q-PHASE 


0 

0 

0 

0 

0 

0 

-2165 

-13 

-2423 

5046 

-963 

3103 

1447 

-5013 

-3357 

808 

-1353 

1402 

2643 

-1708 

263 

-2797 

-3981 

2667 

6034 

1905 

7172 

-3169 

-3368 

-3360 

998 

7554 

-3118 

7618 

3290 

-4383 

-7792 

3059 

-10438 

1436 

3475 

-165 

5756 

-9043 

4038 

-8492 

-2617 

5495 

10589 

-6428 

4938 

-2436 

-2210 

548 

7659 

4528 

1395 

3866 

882 

-2750 

-2975 

13121 

-3807 

3549 

3722 

1643 

-7627 

8580 

-4692 

3303 

-8067 

8794 

-11382 

6634 

-7506 

7780 

-2338 

-4170 

-4515 

-5240 

-6128 

-8631 

-3811 

3307 

6032 

-7914 

9589 

-7476 

735 

-2367 

7895 

1203 

11056 

-1294 

-2050 

2362 

7006 

4862 

5262 

7026 

864 

237 

61 

10301 

-970 

7264 

-2480 

6821 

-4634 

11529 

-8191 

10228 

-8923 

4832 

-9129 

8134 

-12446 

11504 

-4557 

-3308 

-14185 

5133 

-13089 

2978 

1594 

-4736 

-9889 

-5624 

-5367 

-5119 

3099 

-1960 

-4613 

-6557 

-260 

-6594 

638 

0 

789 

-9319 

4442 

-10431 

1041 

1422 

9605 

-6995 

4721 

-1067 

-5283 

2221 

7186 

-1327 

-1023 

494 

2804 - 

- -7160 

3604 

-9 

842 

-6295 

6864 

-1164 

604 

-417 

9633 

183 

210 

6098 

4043 

-2618 

-4186 

7823 

-7198 

5277 

159 

3084 

6102 

-11294 

0 

-1475 

-398 

835 

6011 

-2681 

-1296 

2481 


N.B. The weighting functions are derived from interferograms with 
II2 scans of actual data. The values should be reduced by a 
factor of 2 X 10^ when applied to interferograms. 
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TABLE 6.7 INTERFEROGRAMS EMPLOYED IN THE COMPUTATION 
OF WEIGHTING FUNCTIONS 



CO LAB 

WEIGHT 

CO LaRC/LAB WEIGHT 

CH 4 LaRC 

WEIGHT 


CO 

CH 4 

CO 

CH 4 

CO 

CH 4 


atm- cm 

atm-cm 

atm-cm 

atm-cm 

atm-cm 

atm-cm 

NOMINAL 

.208 

2.60 

.154 

1.78 

.115 

1.44 

TARGET 

.525 

2.60 

.487 

1.69 

.135 

2.38 

OTHER 

.208 

1.74 

.115 

1.44 

.422 

1.46 


.208 

3.48 

.120 

1.55 

.198 

2.48 




.124 

3.50 






.128 

2.93 






.135 

2.38 






.141 

1.83 






.208 

2.60 






.208 

3.48 




TABLE 6 . 

8 LABORATORY 

CARBON MONOXIDE MEASUREMENTS 



MARCH 1974 

NOMINAL CO = 0. 

207 atm-cm 

NOMINAL CH 4 = 2.60 atm-cm 

TARGET CO = 0. 

525 atm-cm 

OTHER CH 4 = 1.74, 3.48 atm-ci 

NUMBER OF SCANS 

1 = 128 

DETECTOR #001 (193 K) 

ACTUAL 

ACTUAL 

MEASURED 

CH 4 

CO 

CO ACCURACY 

atm-cm 

atm-cm 

atm-cm % 

1.74 

0.207 

0.206 -0.5 

1.74 

0.207 

0.205 - 1.0 

1.74 

0.207 

0.209 - 1.0 

2.60 

0 

- .011 

2.60 

0 

- .010 

2.60 

0 

- .011 

2.60 

0.207 

0.185 -10.6 

2.60 

0.207 

0.184 - 11.1 

2.60 

0.207 

0.178 -14.0 

2.60 

0.525 

0.478 - 9.0 

2.60 

0.525 

0.510 - 2.9 

2.60 

0.525 

0.527 0.4 

3.48 

0.207 

0.165 -20.3 

3.48 

0.207 

0.202 - 2.4 

3.48 

0.207 

0.191 - 7.7 
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6. 4. 1.2 Effects of incident light modulation: When the interferometer is 

mounted in a flight vehicle looking down at the earth, 'the energy incident on 
the instrument will vary depending on the type of terrain being covered. The 
fast automatic gain control system (FAGC) is intended to compensate for light 
intensity changes which occur with frequency components up to 20 Hz. 

A set of tests were performed to exercise the instrument FAGC system over 
a range of signal frequencies and amplitude modulations which might be expected 
in flight. These tests were carried out using the multiple chamber test facil- 
ity with the light source modulated by a variable speed chopper. Figure 6.7 
illustrates the source modification made to vary the light intensity. The lamp 
in front of the chopper was adjusted to provide a mean light level for the mea- 
surements. Another lamp behind the chopper illuminated a diffuser which pro- 
vided a light source which could be chopped at a predetermined frequency. The 
amplitude of the chopped light was controlled by varying the voltage to the 
lamp. A gas test condition amounting to 0.208 atm-cm of carbon monoxide and 
2.60 atm-cm of methane was set up in the cells and the lamps were adjusted to 
give a mean light ground albedo of about 30%. The light level as indicated by 
the FAGC panel display was 8.10. A set of measurements of the carbon monoxide 
in the cells was made by co-adding 8 scans of data. The rms measurement noise 
level with this unchopped light condition was .0179 atm-cm of carbon monoxide. 
The chopped light amplitude was then adjusted to give an FAGC reading of 8.13 ± 
1.13 or a light variation of about 14%. Sets of 8 scan measurements were then 
made as the chopper frequency was varied from 5, 10, 15, 20, 30, 40, and 60 Hz. 
The amplitude of the chopped light was then changed to a FAGC value of 8.0 ± 2.0 
or a light variation of ± 25% and the measurements repeated. The rms noise lev- 
els of these measurements were computed and are tabulated in Table 6.10. 

The FAGC system design was intended to compensate for light level changes 
as fast as 20 Hz. However it is seen that the higher frequencies did not cause 
significantly higher measurement noise. The amplitude of the light variation 
appears to have a greater effect on the measurement noise than does the frequen- 
cy of the light variation. Thus in flying over non-uniform terrain with greatly 
differing values of albedo one might expect the measurements to show more scat- 
ter . 


6. 4. 1.3 Variation of measurement noise with indicent light: Most of the 

measurements were performed at a FAGC value of about 8 which corresponds to a 
ground albedo of about 30%. Since the FAGC value during flights had been ob- 
served at times to fall below this value, it was of interest to determine the 
measurement noise associated with lower FAGC values. The incident light was ad- 
justed to give values of 2,4 and 4.06 and a set of 8 scan carbon monoxide mea- 
surements with O.208 atm-cm of carbon monoxide and 2.60 atm-cm of methane in the 
cells was made for each of the light levels. The rms noise levels at the lower 
light levels were as follows: 


FAGC 

RMS NOISE 


(atm-cm) 

8.0 

.018 

4.0 

.059 

2.4 

,065 
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TABLE 6.9 LABORATORY CARBON MONOXIDE MEASUREMENTS 

MARCH 1974 


NOMINAL CO 

= 0.207 

atm-cm 

NOMINAL CH 4 = 

2.60 atm-cm 

TARGET 

CO 

= 0.525 

atm- cm 

OTHER CH 4 = 1, 

,74, 3.48 atm-cm 

NUMBER 

OF 

SCANS = 

128 

DETECTOR #002 

(193 K) 


ACTUAL 

ACTUAL 

MEASURED 


CH 

CO 

CO 

ACCURACY 

atm-cm 

atm-cm 

atm-cm 

% 

1.74 

0.207 

.220 

6.3 

1.74 

0.207 

.212 

2.4 

1.74 

0.207 

.197 

- 4.8 

2.60 

0 

.058 

— 

2.60 

0 

.058 

— 

2.60 

0 

.062 

— 

2.60 

0.207 

.224 

8.2 

2.60 

0.207 

.218 

5.3 

2.60 

0.207 

.218 

5.3 

2.60 

0.525 

.460 

-12.3 

2.60 

0.525 

.468 

-10.9 

2.60 

0.525 

.458 

- 12.8 

3.48 

0.207 

.208 

0.5 

3.48 

0.207 

,190 

- 8.2 

3.48 

0.207 

.206 

- 0.5 


TABLE 6.10 RMS 
WITH 

FREQUENCY 

NOISE LEVELS IN CARBON MONOXIDE MEASUREMENTS 
THE INPUT LIGHT CHOPPED (8 SCANS) 

RMS NOISE (atm-cm) 

(Hz) 

FAGC = 

8.13 ± 1.13 

FAGC = 8.0 ± 2 

0 


.0179 

.0179 

5 


.0282 

.0456 

10 


.0314 

.0491 

15 


.0109 

.0352 

20 


.0260 

.0184 

30 


.0322 

.0410 

40 


.0140 

. .0228 

60 


.0342 

.0130 



GAS BOTTLES WITH 

MANOMETER PREMIXED TEST GASES 


O 03 
lU CJ 

o. cc 
a. 
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Figure 6.7! Multiple Test Facility with Chopped Source 




The noise values do not scale linearly with the light level, but since the com- 
puted noise values have a tendency to show some scatter even for similar sets 
of data, the data do show a general inverse scaling with light level. 


6.4.2 Outdoor Tests. - Following the flight tests of the correlation in- 
terferometer in May 1973, the instrument was returned to the General Electric 
Valley Forge Space Center for a series of outdoor measurements which were then 
used for a determination of correlation functions for carbon monoxide and meth- 
ane. On June 14th and 15th approximately 60 interferograms were obtained for 
various combinations of carbon monoxide and methane in the cells as the atmo- 
speric background amounts of carbon monoxide, methane, and water varied during 
the course of the measurements. 

The procedure used for processing the outdoor calibration data was the fol 
lowing : 


The solar absorption spectra were examined and the carbon monoxide 
lines in the 4.6 ym region, which were found to be the least overlapped 
with water and nitrous oxide lines, were located. The absorption for 
the individual lines was measured and compared with the measured carbon 
monoxide calibration using known carbon monoxide in a gas cell. This 
provided a carbon monoxide measurement for each of the seven carbon mon- 
oxide lines used. An average carbon monoxide level was then determined 
and a standard deviation computed. The primary carbon monoxide . calibra- 
tion was performed with the gas at 760 torr. Since the real atmosphere 
contains carbon monoxide at total pressures ranging from 760 torr to 
very low pressure, and the line absorption for amounts in the range from 
0.1 to 0.2 atm-cm are quite strong, it was believed that some correction 
to the data might be necessary (ref. 7) . As a first order correction, 
carbon monoxide calibrations were obtained for 0.2 atm-cm of carbon mon- 
oxide at 760 and 380 torr, and the ratio of the line absorptions calcu- 
lated. From this data the computed atmospheric values of carbon monoxide 
were increased by a factor of 1.12. A more accurate correction was not 
considered justified since the scatter in the spectral data were much 
greater. The atmospheric carbon monoxide in the path was then plotted 
as a function of time during the day and a smooth curve drawn through 
the data points. This permitted interpolation between the data points 
to obtain carbon monoxide in the atmosphere at times when measurements 
were made with the correlation interferometer. Table 6.11 shows the 
background carbon monoxide as a function of time on June 14-15, 1973. 

The solar spectral data were treated in a similar manner to extract 
information regarding the methane in the atmospheric path. The methane 
spectrum in the vicinity of 3.7 ym was examined for this purpose. Table 
6.11 also shows the background methane as a function of time on June 14- 
15, 1973. 

These spectral data were then added to any carbon monoxide or meth- 
and increment which may have been in the gas cells at the time that the 
interferograms were being collected with the correlation interferometer. 
The interferometer along with carbon monoxide and methane amounts in the 
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path were input to the computer system for computation of correlation 
functions for carbon monoxide and methane. It should be noted that 
not more than 10 test cases were used in computing the correlation 
function which means that there were approximately 50 interferograms 
available to use as tests of the applicability of the correlation 
function. 

The computed methane in the path (using the correlation function 
which appeared to be the best over the complete range of test data) is 
tabulated in Table 6.12. The corresponding carbon monoxide data (for 
atmospheric measurements within the range of applicability of the best 
carbon monoxide correlation functions) are tabulated in Table 6.13. In 
each case the errors in the spectral measurements of the background gas 
are shown. If the errors involved in determining the amount of gas in 
the path spectrally are taken into account, the correlation interferom- 
eter measurements of carbon monoxide and methane agree quite well with 
the known amounts. 

Following the April 1974 flights a series of outdoor measurements were made 
at Hampton, Virginia for calibration purposes. From the matrix of test condi- 
tions in the multiple chamber cells, interferograms were obtained from which 
weighting functions were derived for application to the flight data recorded on 
the Falcon and helicopter flights. Figure 6.8 illustrates the interferograms 
collected during the outdoor calibration. One interferogram is for the amount 
of CO and CH 4 in the atmospheric path with no gas in the cells, i.e., 0.154 
atm-cm CO + 1.78 atm-cm CH 4 . The second interferogram includes the addition of 
0.333 atm-cm of CO for a total amount of 0.487 atm-cm CO + 1.69 atm-cm CH 4 . 

Table 6.6 indicates a weighting function for CO (CO LaRC/LAB) derived from 
a combination of the outdoor calibration interferograms and some laboratory in- 
terferograms which do not have any effect from water vapor (Table 6.7). Table 
6.6 also lists a methane weight (CH 4 LaRC) derived from the Langley calibration 
interferograms . 

A series of flight measurements were made at NASA Langley Research Center 
with the correlation interferometer mounted on a Bell Jet Ranger 206 helicopter 
(Figure 6.9). These flights are described in Section 8 . 

The interferograms collected during the Falcon and helicopter flights have 
the same general shape as those collected during the outdoor calibration tests. 
Figure 6.10 is an interferogram measured during the Churchill flights described 
in Section 7. 
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TABLE 6.11 CARBON MONOXIDE AND METHANE MEASUREMENTS 
USING THE IF -6 INTERFEROMETER 
VALLEY FORGE, PA. JUNE 14-15, 1973 


DATE 

LOCAL TIME 

CO (atm-cm) 

CH 4 (atm-cm) 


6/14/73* 

1100 

.090 ± 

.016 

1.03 ± 

.16 

*Strong 


1255 

.088 ± 

.020 

1.11 ± 

.19 

winds . 


1445 

.107 ± 

.023 

1.16 ± 

.19 



1555 

.119 ± 

.033 

1.11 ± 

.23 



1632 

.167 ± 

.044 

1.28 ± 

.23 


6/15/73 

0909 

.165 ± 

.044 

1.68 ± 

.18 



TABLE 6.12 CORROBORATIVE MEASUREMENT OF METHANE 
SOLAR LOOKING EXPERIMENT 
VALLEY FORGE, PA. JUNE 14-15, 1973 


BACKGROUND + 
CELL INCREMENT 
atm-cm 

1.02 ± .16 


1.09 ± .18 

1.34 ± .23 
1.40 ± .23 
1.68 + .18 


2.08 ± .20 
2.55 ± .20 


3.07 ± .18 
3.31 ± .20 
3.80 ± .20 
4.03 ± .20 

4.50 ± .22 
5.27 ± .23 
6.73 ± .23 


CORRELATION 

INTERFEROMETER 

atm-cm 

0.93 

1.02 

1.03 

1.07 
1.11 
0.87 
0.91 
1.22 

1.25 
1.38 
1.46 

1.54 

1.68 

2.08 

2.45 

2.46 

2.55 

2.90 
2.82 

3.90 
4.09 
4.22 

4.25 
5.24 
6.70 
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TABLE 6.13 CORROBORATIVE MEASUREMENT OF CARBON MONOXIDE 

SOLAR LOOKING EXPERIMENT 
VALLEY FORGE, PA. JUNE 14-15, 1973 

BACKGROUND + CORRELATION 

CELL INCRE^^ENT INTERFEROMETER 

atra-cm atm- cm 


.082 

± .020 

.082 



.108 

.084 

± .020 

.084 

.085 

± .020 

.049 



.085 

.088 

± .020 

.076 

.095 

± .020 

.124 



.144 

.100 

± .023 

.063 



.100 

.122 

± .033 

.112 

.124 

± .033 

,156 

.155 

± .040 

.167 

.202 

± .045 

.260 

.220 

± .045 

.265 

.267 

± .045 

.235 



.250 

.583 

± .020 

.552 

.584 

± .020 

.591 

.593 

± .020 

.558 



.616 

.605 

± .023 

.620 

.630 

± .035 

.680 

.667 

± .044 

.658 



.681 

.682 

± .045 

,657 

1,08 

± ,016 

.991 

1.09 

± .016 

1.06 



1.06 

1.58 

± .020 

1.56 


1.57 

1.57 
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Figure 6.8 Interferogram for In- and Quad-Phase Sample Values for 128 
Co-added Scans for Outdoor Measurements at NASA Langley 
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Figure 6.9 The Helicopter Platform Over the NASA Van 
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INTERFEROGRAM QUAD-PHASE AMPLITUDE INTERFEROGRAM IN'-PHASE AMPLITUDE 



Figure 6.10 Interferogram In- and Quad-Phase Sample Values for 

128 Co-added Scans for the Flights Over Churchill, Canada 
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7. 


ENGINEERING MODEL AIRCRAFT FLIGHT TESTS 


7.1 Platform and Interface 

The Engineering Model Correlation Interferometer was mounted on board a 
Falcon Fan-Jet Aircraft. This aircraft was owned and operated by the Airborne 
Sensing Unit of the Canada Centre for Remote Sensing. It was equipped with two 
windows in the floor which permitted downward viewing. The interferometer was 
mounted over the more forward of these windows which was in line with the door 
of the aircraft. The Falcon is shown in Figure 7.1. The red box in the open 
doorway is the sensor. A mounting plate was made to fit into mounting positions 
above the window. Four shock mounts (Barry Controls Model 2K2-BA-25) were 
placed on the plate. The sensor is shown in position in Figure 7.2. The window 
can be seen below the instrument. The interferometer was mounted on a second 
plate which was mounted on shock mounts. Over the aft window was positioned a 
Panasonic TV camera which looked directly down with a 75° field of view. The 
7° field of view of the interferometer was always within that of the TV camera. 
The TV camera output was recorded on tape which also contained a voice track on 
which could be put orally any comments and information desired during the 
flight . 

The output of the interferometer was, in general, recorded by a Mincom 14- 
channel magnetic tape recorder. The following information was recorded: 

Interferogram data 
2 timing pulses 
Interferometer scan time 
Fast AGC 

Ratio of main channel to FAGC 
Time from time code generator 
Reference interferogram 

The electronics for the interferometer were installed in standard 19” 
racks which were built into the aircraft. These electronics comprised 3-units 
-- the power supply, the processor, and the NOVA interface. The NOVA 1200 
mini-computer and the teletype units were mounted in the aircraft to permit an 
occasional on-the-spot determination of the approximate amount of CO burden and 
a check by use of the standard cell. 

A Systron-Donner Corporation time-code generator of IRIG-B format was also 
mounted on board in such a way that the time -- using a pre-set reference time 
-- was input onto the data tape and visually displayed. The visually displayed 
time could then be orally put onto the video tape. For the helicopter flights, 
to be covered later, a NASA, 36-bit, 100 pulse per second, 1000 cycle carrier, 
time-code generator was used. 


7.2 Data Processing 

Processing of the data from the FM recorder tapes recorded during the 
flights is a two step procedure. The first step is to play the recorded analog 


115 




117 


Figure 7.2 Correlation Interferometer Installed Above Falcon Window 


tapes into a computer through an interface unit to produce another data tape 
which contains the interferogram points and timing data in binary format. This 
binary tape is then available for further processing of the interferogram on 
another computer. 

The four analog data tracks used in the processing are: 1) the digitized 

interferogram samples in NRZ serial form, and also a single 16-bit word indi- 
cating the duration of scan or absolute AGC value which is output during the 
scan plate flyback; 2) a clock track which provides the timing pulses to tell 
the computer when to sample the interferogram data track, the clock track pro- 
vides 24 sampling pulses for each of the 31 interferogram sample points; 3) a 
clock track which provides the timing pulses to tell the computer when to sample 
the duration of scan or absolute AGC data, this clock track provides 16 sampling 
pulses for the single data word; 4) a time code track which provides absolute 
timing information to permit the interferogram data to be related to positional 
information as provided by the TV camera system, the time code is NASA 36-bit 
with a 100 pps code on a 1000 cps carrier. 

The interferogram data track and the two clock tracks are played into the 
computer through an interface which combines the two clock tracks and adds a 
sync word for each frame of data. The computer then checks for a sync word 
within each frame and counts the clock pulses to ensure that the correct number 
of data bits are generated within the duration of a single frame of data. This 

error check eliminates almost all questionable data which may be generated by 

noise pulses on the clock track. 

The time code track is played into a reader which is also interfaced with 
the computer. This adds the time of the data frame in binary form as part of 

the data record for a single scan of the interferometer. 

The digital tape which is thus obtained can then be processed further on 
an SDS-930 computer. The data tape is first unpacked to yield another binary 
tape which formats the interferogram data into the in-phase and quad-phase com- 
ponents, similar to the format output by the NOVA when used with the instrument, 
and separates out the time, duration of scan information and absolute AGC data. 
This tape can now be processed to output interferogram values on a single scan 
basis or after co-adding any specified number of scans. In addition, correla- 
tion functions can be generated on the GE 605 time share computer system using 
previously input calibration interferograms obtained with the NOVA computer on 
paper tape. The correlation functions can then be applied by the SDS-930 to 
the recorded interferogram either on a single scan basis or after co-adding any 
number of scans to yield a measurement of CO or CH as a function of time. 


7.3 Canada Flights - May 1973 

In May 1973, two flights were made on board the Falcon Fan-Jet aircraft. 
These were undertaken mainly to test the operation of the correlation interfer- 
ometer under flight conditions. It was felt that it would also be helpful to 
determine the response of the instrument over various surfaces with different 
albedos. The routes for these two flights are shown in Figure 7.3. Both 
flights generally followed the Ottawa River (with data recording being begun 
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around Petawawa) and then going over the North Bay and Sudbury areas. Flight 
one then went south over Georgian Bay; Owen Sound, Ontario; Perth, Ontario and 
Cornwall, Ontario, to the Montreal area, then back to the Ottawa area; over the 
city of Ottawa, and finally over the area east of Ottawa. Flight two turned 
northeast from Sudbury to Val d'Or, Quebec and back to Ottawa. During flight 
two there was extensive cloud cover and the data were much less useful than in 
flight one. In flight one, the video tape was improperly recorded so that the 
location could only be approximated by use of time and flight speed information 
together with a rough plot of the flight path. 

The useful data were those obtained by use of the NOVA and teletype. Table 
7.1 gives the column density measurements of CO and CH 4 along with the approxi- 
mate locations where these data were obtained. Each of these measurements (de- 
scribed briefly by Goldstein, ref. 38) were averaged over about two minutes. 

The actual measurements are of the total column densities, expressed in atmos- 
phere-centimeters; the quantities noted as ppm are obtained by assuming a con- 
stant mixing ratio throughout the column, undoubtedly a poor assumption, but 
suitable to give a rough idea of the lower atmosphere concentrations. The first 
measurements, particularly over the Petawawa area, were over a cloud covered re- 
gion. Tlie measurements noted as Georgian Bay were over the water of that bay 
and provided very little signal and hence the measurements are not meaningful. 

The areas around Cornwall, Montreal, and, to some extent, Ottawa, are industrial 
areas. The high burdens of CO reflect that condition. The Ottawa area measure- 
ments were made over areas south and west of the city, quite removed from the 
city. The measurements near Alfred show a high level of CO which may or may not 
be true. For these measurements there were no other measurements with which 
these could be compared. Therefore, the accuracy of the measurements cannot be 
determined. 

From these measurements it was concluded that on the first airborne tests 
of the correlation interferometer the instrument operated properly. The measure- 
ments, both of CO and CH 4 were reasonable and showed logical variations with the 
environment . 


7.4 Test Flights - September 1973 

In September 1973 a second series of flights were undertaken in order to 
further test the operation of the correlation interferometer under aircraft 
flight conditions with a variety of atmospheric and surface conditions. The 
scan time of the oscillating plate of the correlation interferometer changed 
both during the flights and the times of calibration. Thus, none of the cali- 
brations could be used nor could any calibration suitable for the treatment of 
the flight data be obtained. In addition to the major cause of error and noise, 
the varying scan velocity, electronic cross-talk developed contributing further 
to the noise. Although both of these problems were subsequently detected and 
corrected, the flight data could not be said to be meaningful. 
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7.5 Canada Flights - April 1974 


A third series of flights were made in April 1974 aboard the Falcon. This 
series consisted of three flights. The first was made in the Churchill, Mani- 
toba area; the second was made in southern Ontario countryside; the third was 
made over Toronto. 

7.5.1 Churchill Area Flight .- The flight in the Churchill area was made 
for the purpose of obtaining data over a large area having a uniform surface, 
thus providing a test for the repeatability of the instrument operation. Al- 
though the flight originated in Winnipeg, Manitoba, the instrument was not put 
into operation until shortly before reaching the actual area selected for the 
repetitive measurements because, due to the very low temperatures above the win- 
dow and, hence, around the instrument, the scan time was significantly lower 
than normal operating values. Even during about the first third of the first 
pass this condition prevailed. Three passes were made over the same area, as 
shown in Figure 7.4. This area was almost completely covered with fresh snow 
which could be distinguished from older (presumably recrystallized) snow by its 
much greater reflectivity (of the order of 10% and 25%, respectively). The ter- 
rain in this area was essentially free of any vegetation. Most of the rivers 
and creeks were still frozen and snow-covered, although the ice in the center 
of some of the larger rivers had started to break up, giving the only variation 
in albedo. However, such covered only a small part of the field of view which, 
for the altitude of 31,000 feet, was about 3700 feet, and thus the variation 
was not significant. In comparing the data all of the straight portions of the 
three passes was used except for the first one-third of pass one. That was 
eliminated because of the still increasing scan time and because this was over 
a slightly different area, some of which had significant vegetation. 

The weighting functions used in obtaining the CO and CH 4 burdens from the 
interferograms were ones obtained from the use of laboratory data only. These 
had essentially no water in the path. For the cold atmosphere of the Churchill 
region this was thought to be most suitable. For all subsequent measurements 
discussed in this report, the weighting functions employed several points of 
calibration data obtained at NASA Langley in May 1974, and several points ob- 
tained in the laboratory. It is believed that these spanned the burdens of wa- 
ter and other gases encountered in the measurments. 

The results of these measurements is shown in Table 7.2. These are pre- 
sented as single-scan data, and as data averaged over a number of scans — 2 , 

8 , 16, and 32. The average absolute deviation from the average is given in each 
case. For CO this is 44% for 1 scan, 31% for 2 scans, 18% for 8 scans, 12% for 
16 scans and 10% for 32 scans. For CH 4 the average deviation is 10% for 1 scan 
and 2% for 32 scans. The single-scan data for both CO and CH 4 for the three 
passes of the Churchill flights are given in Appendix I as an example of data 
of this type. 

The values of the total burdens (in atm-cm) of CO and CH 4 was converted to 
values of concentration (in ppm) by assuming a constant mixing ratio thoughout 
the atmosphere and by taking into account the total amount of gas in the path, 
here equivalent to 2.1 vertical columns of atmospheric gas. This is a reason- 
able assumption in this remote area although the mixing ratio is likely to de- 
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TABLE 

7.1 FLIGHT 

MEASUREMENTS OF 

CO AND 



CH 4 OVER CANADA 

, MAY 24 

, 1973 



CO 


CH 4 




atm-cm 

ppm 

atm-cm 

ppm 


Petawawa 

.074 


1.50 

0.84 

80% cloud cover 

North Bay 

.176 


1.45 


40% cloud cover 

Sudbury 

.216 

.141 

1.43 

0.84 

10 % cloud cover 

Georgian Bay 

.091 

.060 

1.09 

0.71 

Over water 

Owen Sound 

.275 

.180 

1.63 

1.07 


Perth 

.148 

.094 

1.62 

1.03 


Cornwal 1 

.486 

.307 

2.53 

1.59 


Montreal Area 

.479 

.296 

2.08 

1.28 

Some water 

Ottawa Area 

.149 

.104 

1.38 

0.96 

Descending 

Ottawa 

.389 

.288 

1.72 

1.27 


Alfred 

.431 

.315 

1.63 

1.19 

Bog 

*Clouds at 8000 

ft. 






TABLE 7.2 MEASUREMENTS IN THE CHURCHILL AREA 


CO CH 4 

atm- cm atm- cm 


32 

SCANS 
Pass 1 

.278 

± .048 

3.71 

+ 

0.08 


Pass 2 

.275 

± .011 

3.74 

+ 

0.03 


Pass 3 

.267 

± .022 

3.74 

+ 

0.06 

16 

SCANS 
Pass 1 

.279 

± .053 





Pass 2 

.275 

± .011 





Pass 3 

.267 

± .032 




8 

SCANS 
Pass 1 

.279 

± .060 





Pass 2 

.271 

± .035 





Pass 3 

.267 

± .049 




2 

SCANS 
Pass 1 

.278 

± .096 





Pass 2 

.271 

± .077 





Pass 3 

.267 

± .082 




1 

SCAN 
Pass 1 

.278 

± .135 

3.71 

+ 

0.33 


Pass 2 

.271 

± .103 

3.74 

+ 

0.41 


Pass 3 

.267 

± .125 

3.74 

+ 

0.32 
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crease in the stratosphere. Thus the concentrations in the stratosphere should 
possibly be slightly lower than given by a constant mixing ratio. The values 
obtained with this assumption are compared with other measurements made in the 
Arctic in Table 7.3. The measurements of Wilkness (ref. 39) were made in the 
Arctic Ocean area at approximately latitude 70 - 80 °N, and those of Seiler 
(ref. 40) were made in the polar area. Thus these were at higher latitudes than 
were those of this work (57°40' - 58°25' N) . It can be seen that the agreement 
for both CO and CH 4 is good. The Churchill measurements show that the correla- 
tion interferometer gives consistent values of the gas burden. 

7.5.2 Southern Ontario Flights .- The flights in the Southern Ontario area 
were made for the purpose of determining the effect of rapidly varying surface 
characteristics on the measurements. In this area the surface cover does vary 
rapidly with various types of vegetation, fields plowed in various ways, roads, 
rivers and creeks, etc. It was planned to traverse the same path several times 
each with a different value of V/H (speed t altitude in sec-^) by changing the 
values of one or both of these parameters. The actual six passes are shown on 
the map in Figure 7.5. The paths are actually different because of the diffi- 
culty in exactly repeating a path in a high-speed aircraft and because of a re- 
quest by air-traffic control, after the first three passes, to change the area. 
The changes did not change the surface characteristics appreciably. 

For each of the six passes a region which was judged to encompass the most 
rapidly varying surface characteristics was chosen for comparision. In Table 
7.4 these data for each of the six passes are given in terms of the average of 
single scan data and the average absolute deviation therefrom. The values of 
V/H for each pass is given. It is seen that the average absolute deviation does 
not increase with V/H even though the values of V/H encountered were more than 
twice as high as the maximum value designed for (.05 sec-^) . It is concluded 
that the instrument can make consistent measurements over rapidly varying ter- 
rain with values of V/H in excess of 0.1 sec-^. 

7.5.3 Toronto Flights .- The flights over Toronto were made at the request 
of and in conjunction with the Ontario Ministry of the Environment. Nine passes 
were made in various regions over the city and its environs. These are shown in 
Figure 7.6. The first pass and part of several others were made over the waters 
in Lake Ontario and a very low signal was obtained. The results are presented 
as 32-scan averages in the map in Figure 7.6. It can be seen that the higher 
values were found in the downtown region and along the heavily traveled free- 
ways, the Queen Elizabeth Way (along the South shore) and Rt. 401 (along the 
north side of the most populous area) . The area to the north shows lower values 
as do, to a lesser extent, the area along the pass which more or less followed 
the Don River Valley. 

Point samples were obtained by the Ontario Ministry of the Environment who 
provided the data listed in Table 7.5. 
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TABLE 7.3 ARCTIC MEASUREMENTS (ppm) 



CO 

CH 4 

THIS WORK* 

0.14 

1.89 

WILKNESS 



TROPOSPHERE 

WINTER 

SUMMER 

0.22 

0.09 

1.62 

1.38 

STRATOSPHERE 

0.12 

1.59 


SEILER 

STRATOSPHERE O.ll 

*Based on constant mixing ratio througout 
the atmosphere. Winter conditions pre- 
vailed at time of these measurements. 


TABLE 7.4 SOUTHERN ONTARIO FLIGHT DATA 



CO 

V/H 

PASS 

(atm -cm) 

(sec“^) 

1 

.248 ± .079 

.014 

2 

.270 ± .111 

.045 

3 

.261 ± .106 

.029 

4 

.277 ± .103 

.044 

5 

.237 ± .076 

.059 

6 

.232 ± .106 

.121 


125 



126 

















TABLE 7.5 


ONTARIO MINISTRY OF THE ENVIRONMENT CO (ppm) DATA 


Station 1100 1200 1300 


A 

B 

C 

D 

E 


2 2 1 

3 3 3 

1 2 1 

2 

12 2 


These results show variations which follow the same trends as do the re- 
mote data. Of course, these data are coarse and no direct comparison can be 
made between column density measurements and ground samples. It is planned by 
the Ontario Ministry of the Environment to make computer runs using a dispersion 
model to calculate what values each of these two types of measurements should 
give and thus how they compare with each other and with the models. At present 
it can merely be said that the results over polluted areas are reasonable. 
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8 . ENGINEERING MODEL HELICOPTER FLIGHT TESTS 


8.1 Platform and Interface 

As a further test of the Engineering Model Correlation Interferometer, the 
instrument was mounted on board a Bell Jet Ranger 206 Helicopter owned and op- 
erated by Keystone Helicopter, Inc. of West Chester, Pa. This aircraft is shown 
in Figure 8.1. Keystone provided a downward- looking port in the helicopter for 
these tests. A pyrex window was installed in this port for the CO and CH 4 mea- 
surements in the 2.3 pm region. A mounting on which the interferometer shock 
mounts were placed was designed and fabricated and the instrument was installed 
thereon in a position over the rear left seat and the floor in front of it. A 
Sangamo Model 3568 14-channel tape recorder was used for this work and mounted 
on the same frame as the shock -mounted interferometer. The sensor and the re- 
corder are shown mounted on the helicopter in Figure 8.2. The electronics. Fig- 
ure 8.3, including a time-code generator, were placed in the luggage compart- 
ment. The time-code generator and a clock on the dashboard of the helicopter 
were set so that the readings were coordinated. A Sony Model AVC 3400 video 
tape recorder was placed in the nose of the helicopter with a downward-looking 
68 ° field of view which included the 7° field of view of the instrument. This 
recorder also had a voice channel using a built-in microphone so that the oper- 
ator in the left- front seat could put times, locations, speeds, altitudes, and 
other comments on the tape. With the microphone located several feet from the 
operator much of the helicopter noise was picked up. Thus, on the recorded 
tapes, the voice was somewhat difficult to pick up but was, in general, satis- 
factory. 


8.2 Pennsylvania Flights 

Two short flights were made over a suburban area about thirty miles west 
of Philadelphia in order to test the operation of tiie correlation interferom- 
eter using the helicopter as a platform. The paths for these flights are shown 
in Figure 8.4. The measurements obtained are given in Table 8.1. It can be 
seen that these data, which are averages and average absolute deviation of sin- 
gle scan measurements, show precisions about the same as that of the data ob- 
tained on the jet aircraft -- of the order of 40%. The one exception to this 
is the value obtained at a high value of V/H (about 1 sec-^) and is presumed to 
be due to the magnitude of V/H which is more than twenty times the V/H for which 
the instrument was designed. 

8.3 NASA Langley Area Flights 

A series of flights was made in the vicinity of NASA Langley with the ob- 
jective of obtaining correlation interferometer measurements to compare with 
other measurements obtained with ground-based instnmentation to provide both 
comparative and calibration data. These flights were also desired in order to 
further test the operation of the instrument on the helicopter, to test the re- 
peatability of the instrument response and to test the response and operation 
over various terrains including areas of high and low reflectivity. 
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Figure 8.2 Correlation Interferometer and Recorder Installed in the Helicopter 




32 


Figure 8.3 Electronics Installed in Helicopter Baggage Compartment 




Figure 8.4 Helicopter Flight Paths - June 1974 





TABLE 8.1 HELICOPTER MEASUREMENTS IN PA. 


CO 

(atm- cm) 


CH 4 

(atm -cm) 


.240 


.084 

1.82 

± .22 

.208 

+ 

.088 

1.86 

± .21 

.233 

± 

.071 

1.32 

± .29 

.228 

+ 

.086 

1.45 

± .33 

.229 

+ 

.067 

1.43 

± .28 

.226 

+ 

.096 

1.17 

± .38 

.242 

+ 

.193 

1.51 

± .69 


v/H 1 


Ground-based column density measurements were obtained with NASA's Idealab 
IF -6 interferometer of 0.1 cm-^ resolution operated by P. J. LeBel. Many passes 
were made over the van containing that instrument and a chromatograph with which 
ground-level concentrations were obtained. The IF -6 interferometer integrated 
signals over a period of about two minutes, and transformed the interferograms 
into spectra, and obtained the CO and CH 4 burdens therefrom. For each of the 
many passes of the helicopter over the van. Table 8.2 gives data averaged over 
a uniform area, shown in Figure 8.5, east and west of the van, which, in gen- 
eral, consisted of 10 to 19 scans except for pass number 18 which was obtained 
by hovering for 210 scans, and the data obtained by the instrumentation in the 
NASA van. The GE and NASA data are also shown graphically -- the CO for June 8 
and 9 in Figures 8.6 and 8.7, respectively, and the methane for June 8 and 9 in 
Figures 8.8 and 8.9, respectively. The agreement is, in general, quite good -- 
with the error limits which are of the order of 20% for each instrument. The 
altitude of the helicopter was only a few hundred feet above the ground so that 
this part of the radiation path added very little to the total burden as com- 
pared with the interferometer looking up from the ground. 

The scatter in the correlation interferometer measurements was found to be 
about the same as those obtained when measurements were made using the Falcon 
Jet as the platform. Measurements made during the six-minute period of hovering 
(Figures 8.10 and 8.11) show no increase in the CO or the CH 4 burden indicating 
no effect of the exhaust. The helicopter thus appears to be a suitable as well 
as convenient and useful platform. 

Further measurements were made during this series of flights. However, 
difficulties arose in the form of frequent dropout of interferogram data and 
culminating finally in the finding of a broken wire. Hence data obtained after 
the first appearance of this trouble were not deemed usable. The broken wire 
was more likely due to rough ground transportation rather than the comparatively 
smooth flights. The only problem which may have been introduced by roughness in 
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TABLE 8.2 DATA FROM PASSES OVER VAN AT NASA LANGLEY 


DATE 

6/18/74 


6/19/74 


GE NASA 


RUN 

RUN 

TIME 

GE 

NASA 

GE 

NASA 

NO. 

NO. 

HR: MIN: SEC 

CO 

CO 

CH4 

CH4 


1 

9:25 


.17 


1.97 


2 

10:00 


.15 


1.72 


3 

11:00 


.105 


1.36 


4 

12:50 


.09 


.90 

1 


14:03:57 

.24 


.86 


2 


14:06:30 

.17 


1.01 


3 


14:09:41 

.15 


.84 


4 


14:11:50 

.13 


.86 


5 


14:15:23 

.15 


.82 


6 


14:18:01 

.13 


1.09 


7 


14:21:46 

.13 


1.05 


8 


14:27:32 

.13 


.96 


9 


14:46:52 

.19 


1.05 


10 


14:50:40 

.15 


1.20 


11 


14:53:45 

.17 


.93 


12 


14:56:51 

.10 


.87 


13 


14:59:31 

.20 


.99 


14 


15:02:49 

.18 


1.05 


15 


15:05:10 

.22 


1.34 


16 


15:07:35 

.14 


1.02 


17 


15:15:33 

.14 


1.33 



6 

15:30 


.11 


1.55 


7 

16:00 


.12 


1.67 


8 

16:15 


.13 


1.82 


9 

16:30 


.15 


1.95 

18 


16:31:30 

.14 


.88 


19 


16:36:36 

.20 


.56 


20 


16:40:12 

.14 


.76 



10 

16:45 


.17 


2.03 


11 

16:49 


.18 


2.22 


1 

9:10 


.19 


2.13 


2 

9:55 


.15 


1.73 


3 

10:20 


.14 


1.62 


4 

10:35 


.14 


1.42 


5 

10:45 


.14 


1.46 

21 


10:50:20 

.16 


1.08 


22 


10:52:13 

.15 


.88 


23 


10:59:40 

.14 


.90 



6 

11:00 


.13 


1.34 

24 


11:01:40 

,11 


.85 


25 


11:05:10 

.10 


.72 


26 


11:07:03 

.13 


.95 



7 

11:33 


.16 


1.30 


8 

13:00 


.12 


1.18 


9 

13:25 


.12 


1.19 


10 

15:00 


.13 


1.36 


11 

15:30 


.16 


1.50 


12 

15:40 


.16 


1.64 
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Figure 8.5 Area of Flights at NASA Langley 
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Figure 8.6 CO Burdens - Flight 1-8 June 1974 
(Average Over Van Area - NASA LaRC) 



Figure 8.7 CO Burdens - Flight 2-9 June 1974 
(Average Over Van Area - NASA LaRC) 
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Figure 8.8 CH^ Burdens - Flight 1 - 8 June 1974 
(Average Over Van Area - NASA LaRC) 
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Figure 8.9 CH 4 Burdens - Flight 2-9 June 1974 
(Average Over Van Area - NASA LaRC) 






flight is thought to be a possible problem with the recorder at the higher 
speeds used in flights away from Langley since the tape recorder was not shock- 
mounted at that time. All the data on the later flights were suspect. These 
flights included ones over Hampton, Newport News, the Hog Island area, the Dis- 
mal Swamp, and the beach area south of Virginia Beach. The data for these were 
therefore not used. 

The main objective of these flights was met. Data were obtained from the 
helicopter platform. These data were comparable with those data obtained from 
the ground by the IF-6 interferometer. These data gave acceptable comparison. 


8.4 Washington Area Flights 

A series of flights were made in the Washington, D.C. area in order to 
test the ability of the correlation interferometer on board a helicopter to 
make meaningful measurements in the urban areas and in doing so to compare such 
measurements with those made by various other instruments in the same area at 
about the same time. In addition to the correlation interferometer, column 
density data were also obtained with the NASA interferometer and sample data 
were obtained with the NASA chromatograph at ground level, the NRL chromato- 
graph at ground levels and at various altitudes including above the mixing lay- 
er, and the State of Maryland infrared gas analyzer which obtained ground-level 
measurements . 

A direct comparison of the correlation interferometer data with that from 
the IF-6 can be made. The part of the path of the radiation between the ground 
and the platform adds only a few percent as compared to the path between the sun 
and the ground. The comparative data are shown in Table 8.3 and in Figures 
8.12, 8.13, 8.14 and 8.15. The GE data are averages for periods of 15 seconds 
over a region extending from the van in a direction away from an adjacent high- 
way. The NASA data were integrated over periods of one to two minutes. The 
agreement is again seen to be quite good. The NRL data shown are total burdens 
calculated for a vertical column using the NRL sample data (taken over a period 
of about an hour) and a logical extrapolation into the stratosphere. 


8.5 Other Flights 

These flights have shown that the helicopter is a very suitable platform 
for tests of and measurements by the correlation interferometer just as is the 
Jet aircraft. A series of flights was made in September, 1974 for Columbia Gas 
Systems to test the ability of the correlation interferometer to detect natural 
gas leaks by measurement of ethane, methane, and the ratio of these two gases. 
For this work the instrument was converted to operate in the 3.3 ym range. The 
flights took place in rural areas of Kentucky, Tennessee, Virginia, and West 
Virginia. The results will be reported elsewhere at a later date. 
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TABLE 8.3 COMPARISON OF HELICOPTER DATA WITH 
GROUND -TRUTH MEASUREMENTS 



GE 

RUN 

NASA 

RUN 

START TIME 

STOP TIME 

GE 

NASA 

GE 

NASA 

DATE 

NO. 

NO. 

HR:MIN:SEC 

HR: MIN: SEC 

CO 

CO 

CH 4 

CH 4 


7/11/74 


1 

15:10:00 

15:11:16 


.10 


1.43 


6-1 


15:31:02 

15:31:15 

.17 


1.05 




2 

15:34:00 

15:36:32 


.12 


1.57 



3 

16:05:00 

16:07:32 


.14 


1.78 


6-2 


16:06:33 

16:06:52 

.21 


0.95 




4 

16:45:00 

16:46:16 


.16 


2.48 

7/12/74 


1 

09:23:00 

09:24:16 


.13 


2.31 



2 

10:03:50 

10:05:06 


.14 


1.84 



3 

10:17:35 

10:18:51 


.16 


1.80 



5 

10:49:00 

10:50:16 


.15 


1.61 


7-11 


11:01:15 

11:01:31 

.20 


1.08 



7-12 


11:03:05 

11:03:17 

.17 


1.10 



7-13 


11:06:06 

11:06:22 

.21 


1.11 



7-15 


11 : 10:02 

11:10:14 

.20 


1.18 




6 

11:10:50 

11:12:06 


.11 


1.52 


7-16 


11:11:41 

11:11:53 

.14 


1.05 



7-17 


11:14-23 

11:14:38 

.18 


1.25 



7-19 


11:17:19 

11:17:35 

.21 


0.97 



7-20 


11:19:08 

11:19:18 

.16 


1.26 



8-1 


11:32:30 

11:32:37 

.19 


1.22 




7 

11:44:50 

11:46:06 


.14 


1.43 



8 

12:27:20 

12:28:36 


.13 


1.43 



9 

13:51:05 

13:52:21 


.14 


1.50 


9-1 


14:19:57 

14:20:12 

.26 


1.31 




10 

14:20:45 

14:22:01 


.17 


1.67 
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Figure 8.12 CO Burdens - D. C. Flights - 11 July 1974 
(Average Over Van Area) 
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Figure 8.13 CO Burdens - D. C. Flights - 12 July 1974 

(Average Over Van Area) 
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Figure 8.14 CH 4 Burdens - D. C. Flights - 11 July 1974 
(Average Over Van Area) 




Figure 8.15 CH 4 Burdens - D. C. Flights - 12 July 1974 
(Average Over Van Area) 



9. A STUDY OF THE FEASIBILITY OF THE REMOTE MEASUREMENT OF 
TRACE ATMOSPHERIC SPECIES BY CORRELATION INTERFEROMETRY 


A correlation interferometer was developed for the measurement of carbon 
monoxide and methane in the troposphere and in the stratosphere using absorption 
in the 2.35 ym region. Based on positive results obtained with this instrument, 
a study was undertaken to determine the feasibility of using this instrument for 
a large variety of atmospheric trace gases. This study will be discussed in de- 
tail after a brief discussion of the instrumental technique. 

The 2.35 ym band of carbon monoxide was chosen since in this region, the 
non-thermal infrared, the main contribution to the observed radiation is that of 
the light reflected from the earth's surface. This makes the contribution of an 
absorbing molecule essentially independent of its altitude and the interpreta- 
tion of the data is simpler than in the thermal infrared. For solar-looking 
measurements this consideration is negligible. It is not always possible to 
find suitable absorption bands in the non-thermal infrared and thermal infrared 
bands must be used. This is the case for three of the gases examined in the 
multi-gas analysis given below. 

The basic interferometer is shown schematically in Figure 3.1. It is to 
be noted that the two mirrors are fixed in the instrument and the variation in 
path length difference is attained by the oscillation of the compensating plate. 
For the results given below, a thirteen degree oscillation was used giving a 
path length difference from 2.7 to 3.95 ym. 

Figure 3.3 shows the spectrum of carbon monoxide in the 2.35 ym region and 
the corresponding interferogram. Since the spectral lines have a separation of 
the order of 3 cm-^, the interferogram has a peak at 1/3 cm-^ or about 3 mm. 

The correlation interferometer uses only this part of the interferogram since 
most of the information on CO is contained therein (2.7 to 3.95 mm). The in- 
strument operates directly on this interferogram section and does not transform 
it into a spectrum. Indeed, a spectrum obtained from the transformation of a 
portion of an interferogram would have little meaning. 32-scan measurements in 
the laboratory and outdoors show accuracies with approximately 10% error for CO 
and 2% error for CH4. 

For the measurement of a specific gas the optimum wavelength range and de- 
lay range must be known. The wavelength range determines the needed spectral 
filter, detectors and optics. The optimum wavelength range is determined by 
the absorption spectra of the gas and of other (interferent) gases of the atmos- 
phere, by the identity and strength of the light source, by detector sensitiv- 
ity, and to a slight extent, by the availability of optical materials. The 
breadboard and engineering models of the correlation interferometer have em- 
ployed optical filters with a full width, at half maximum transmittance, of the 
order of 45 cm-^, centered at about 4280 cm”^ (2.34 ym) . Only the radiation 
transmitted by this filter affects the interferogram. Each gas has a specific 
spectral range which is best to use, although these may be very nearly the same 
for some gases with the effects still being separable by the correlation inter- 
ferometer technique, since the interferogram on which the direct measurements 
are made is determined by line spacing which may vary significantly while still 
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being in the same wavelength range. This is the case with methane and carbon 
-monoxide. ^ 

The other major parameter to be set is the delay range. The line spacings 
of the absorption spectra for each species is somev/hat different. For measure- 
ment it is very convenient if these are such that the interferometric effects 
are nearly but not exactly the same, as with CH 4 and CO, as mentioned above. 

It is also quite convenient if the interferometric effects are nearly the same 
— that is, occur over the same delay range -- but the spectral effects are 
found in a somewhat different spectral range. This occurs when the rotational 
energy level differences of two gases are nearly equal but the vibrational level 
differences are not nearly equal. It is not entirely fortuitous that this is 
often the case -for similar molecules. Thus examination of spectra and interfer- 
ograms for ammonia, NH 3 , and nitric oxide, NO, as well as methane and carbon 
monoxide, each show major effects on the interferogram in the 2.7 to 4.0 mm 
range for spectral absorptions which are in different wavelength ranges. Thus 
these four gases could be measured using the same delay range, i.e., the same 
difference in the path lengths of the two arms of the interferometer, just by 
the use of different optical spectral filters. 

For some other gases, optimum results would be obtained using a different 
delay range. However, it can be expected that any other delay range can be used 
for more than one other gas, again just by the use of appropriate spectral fil- 
ters. Later the optimum wavelength ranges and path differences (delay ranges) 
are determined for several gases. However, it will also be seen that there is 
a single delay range which is suitable for at least most of the gases of inter- 
est. 


In the first study the following steps comprised each calculation: 

1. Spectra with 0.5 cm~^ resolution were obtained for indivi- 
dual gases at specific optical thicknesses using the Digi- 
lab FTS-14 interferometer-sepctrometer and ratioing to the 
background spectrum. 

2. By use of computer techniques, spectra of various combina- 
tions of gases were obtained with appropriate optical fil- 
ters . 

3. From each spectrum, an interferogram was generated. 

4. For each gas of interest an interferogram for the combina- 
tion of the interferent gases and that gas was compared 
with an interferogram for the combination of the inter- 
ferent gases without the gas of interest. The difference 
is the effect of the gas of interest in the presence of 


interferents . 


The gases studied were: 

CH 4 

Methane 

CO 

Carbon Monoxide 

CO 2 

Carbon Dioxide 

H 2 O 

Water 
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NHo 

NO 

NO2 

N2O 

SO2 

C2H6 

C2H4 


Ammonia 
Nitric Oxide 
Nitrogen Dioxide 
Nitrous Oxide 
Sulfur Dioxide 
Ethane 
Ethylene 


Table 9.1 lists the spectra which were obtained, arranged in order of increasing 
wavelength. 


In addition to the spectra of these gases, spectra of water are needed. 
Since it has not been practical, in the laboratory, to obtain spectra for atmo- 
spheric amounts of water (for which a path-length of several kilometers is 
needed) these were obtained by use of a computer. All spectral regions where 
water is a significant interferent in the spectra given above were covered. The 
amount of water used was 2150 atm-cm, as given by Gutnick (ref. 9) (1962). 
Multiples of this from 0.1 to 10 were used. 


Various combination spectra were generated, from which the corresponding 
interferograms were obtained. For these calculations a Lorentzian filter was 
assumed and filter width was taken of the order of 1/2 to 1 % of its center fre- 
quency. 


The interferograms with and without the gases of interest (the gases to be 
measured) are shown in Figures 9.1 through 9.11 for atmospheric burdens of CH 4 , 
CO, CO 2 , H 2 O, NH 3 , NO, NO 2 , N 2 O, SO 2 , C 2 He, and C 2 H 4 , respectively. The effect 
of the gas of interest is shown by the difference in the two interferograms 
plotted for that gas. Quantitative estimates of the sensitivities and accura- 
cies for the gases are obtained by comparing the interferogram showing the CO 
effect (Figure 9.4) and the experimental accuracies of CO measurements from 
studies such as those mentioned above. Table 9.2 gives the estimated sensitiv- 
ities relative to CO. Thus, correlation interferometry can be a very useful 
technique for the stratospheric and tropospheric measurement of these important 
trace atmospheric species. 

To further study the feasibility of the measurement of several of the at- 
mospheric trace gases listed above, in situations simulating their optical 
thicknesses in the atmosphere, a second systematic study has been carried out 
to simulate the performance of the correlation interferometer in the qualitative 
measurement of these gases. Because of the extensive computer calculations re- 
quired, the complete study was not carried out for all of the gases. 

Each calculation consisted of the first three steps of the calculation de- 
scribed above and the following three steps: 

4. Various combinations of these interferograms were used to 
obtain weighting functions. 

5. Using one or more appropriate weighting function(s), each 
interferogram was tested to determine the optical thick- 
ness the correlation interferometer technique would indi- 
cate. 
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TABLE 9.1 SPECTRA USED FOR INTERFEROGRAMS 
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No significant absorption. 
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Figure 9.1 The Effect of CH 4 on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 





Figure 9.3 The Effect of CO 2 on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 




Figure 9.4 The Effect of H 2 O on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 





Figure 9.5 The Effect of NH 3 on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 






Figure 9.7 The Effect of NO 2 on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 
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Figure 9.8 The Effect of N 2 O on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 
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Figure 9.11 The Effect of C 2 H 4 on the Interferogram Signal 

(Arbitrary Units) as a Function of Path Difference 


TABLE 9.2 

ESTIMATED SENSITIVITIES - SOLAR 

LOOKING ■ 


RELATIVE TO CO 

(CO = 0.04 atm- cm) 


FILTER 

CENTER 

DELAY RANGE 

RELATIVE 

SPECIE 

FREQUENCY 


BURDEN 


cm" ^ 

(ym) 

mm 


CH 4 

4280 

(2.34) 

5.5 - 6.4 

5 




All Other 

2 

CO 

4280 

(2.34) 

2.3 - 3.5 
5.6 - 6.3 

7.4 - 9.4 

1 

CO 2 

4835 

(2.07) 

5.0 - 6.0 

2 




6.0 - 8.0 

125 




All Other 

1250 

H 2 O 

3465 

(2.89) 

All 

1250 

NH 3 

4500 

( 2 . 22 ) 

2.6 - 4.0 

5 




6.3 - 7.5 

10 

N 2 O 

3465 

(2.89) 

3.0 - 3.6 

5 




8.5 - 9.7 

.5 

NOt 

1900 

(5.26) 

2.8 - 3.7 

.2 




5.8 - 7.0 

.08 




9.0 - 9.5 

.08 

N 02 t 

1630 

(6.13) 

0.5 - 4.0 
All Other 

.06 

S 02 t 

1370 

(7.30) 

'V 2.0 

.03 




2.5 - 4.0 

.2 




4.5 - 5.7 

.2 




8.0 - 9.5 

.3 

C 2 H 4 t 

2988 

(3.35) 

1.0 - 4.0 

.5 




4.5 - 5.5 

.2 




6.0 - 8.5 

.5 

C 2 H 6 + 

2988 

(3.35) 

1.0 - 4.0 

.02 




5.0 - 7.5 

.008 




7.5 - 9.2 

.02 


tDetector with a D* = 3 x 10^° cm Hz^/^ w-^; all other 
species assume detector with a D* = 10^^ cm Hz^'^ w“^. 



6 . The indicated optical thicknesses were compared with the 
actual optical thicknesses to assess the accuracy of the 
technique and its sensitivity. 

From this study conclusions can be drawn concerning the best delay range, 
and best spectral range, and the expected accuracy and sensitivity for the mea- 
surement of each of the gases. These conclusions are shown in Table 9.3. 

A critical examination of the results of these two types of calculations, 
as given in Tables 9.2 and 9.3, shows the correlation interferometer to have 
the needed sensitivity and accuracy to measure atmospheric amounts of most of 
the atmospheric trace species. 


TABLE 9.3 CALCULATED ACCURACIES AND SENSITIVITIES FOR MEASURE- 
MENT OF ATMOSPHERIC TRACE GASES AS DETERMINED BY 
APPLICATION OF CORRELATION INTERFEROMETER ALGORITHM 


GAS 

SPECTRAL 

DELAY 

RANGE 

ACCURACY FOR 

SENSITIVITY 


REGION 



ATMOSPHERIC 



(cm"^) 

(mm) 

\ 

BURDEN (%) 

(atm- cm) 

CH 4 

4280 

2.58 - 

3.83 

10 


CO 

4280 

2.58 - 

3.83 

10 

.01 

H 2 O 

4865 

2.58 - 

3.83 

25 


NO 

1900 

5.78 - 

7.19 

10 

.001 



2.58 - 

3.83 

20 

.002 

N 2 O 

3465 

2.58 - 

3.83 

10 

.1 

SO 2 

1370 

8.28 - 

9.69 

15 

.002 



2.58 - 

3.83 

25 

.005 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 1 


SCAN 

CO 

CH4 

SCAN 

CO 

CH4 

NO. 

(atm-cm) 

(atm-cm) 

NO. 

(atm-cm) 

(atm-cm) 


. 4414 

3, 7672 

544 

.3136 

3. -9275 

'+95 

• 4986 

4.4-806 

545 

• 4003 

3.9303 

496 

.5834 

2.6820 

546 

• 55Q6 

3,b8b8 

497 

• 546 1 

3.i>357 

548 

•3609 

3 . 5575 

498 

. 1959 

4.0735 

549 

• 3490 

3.7733 

5Q0 

.2991 

3.3077 

55l 

• 3918 

2.9152 

bCl 

.2508 

3.4481 

55? 

• 2251 

3 . 6 8. 4 1 

5(j? 

.5335 

3.5612 

553 

• 3788 

3,1893 

b03 

.0699 

4,u531 

554 

• 3182 

3, 7687 

bQit 

.3953 

3.3726 

556 

.3594 

4,2140 

bo 5 

.5877 

3. 7s 86 

557 

• 3791 

3 . 3 .3 9 0 

b06 

.1340 

3.7339 

55,8 

• 3199 

3 . 33*+ 1 

bo 8 - 

.4385 

3.8813 

559 

• 5116 

3,8152 

b09 

.5378 

3.7410 

560 

• 4895 

3.3166 

bl 1 

.1576 

3.3294 

561 

• 4313 

3,3375 

bU? 

.226 7 

3.4253 

562 

• 4729 

3,9021 

bl 3 

. 4746 

3.1921 

564 

• 4199 

4,3725 

blu 

.5798 

3. 7908 

565 

• 2257 

3.8922 

bl6. 

.3726 

3 . b7 76 

566 

• 3114 

3.8468 

bl7 

. 4Q56 

3.6653 

567 

• 1593 

3,6139 

bl8 

.3220 

3.6270 

568 

• 3274 

3, 6792 

bl9 

. 4874 

3,1220 

569 

• 3465 

3.9123 

bHu 

.4473 

3 .2826 

570 

•2569 

3.9997 

b2l 

.4193 

3.7929 

572 

.4523 

3.5903 

bHP 

.3614 

3.4855 

573 

• 3978 

3.4034 

bE4 

.2707 

3.4323 

57 4 

• 47 86 

3 • 56 05 

bE5 

.4848 

3.7762 

575 

• 3902 

3.6289 

b26 

.5151 

3 • 1855 

576 

• 2325 

3 . 0639 

b27 

.3613 

3.7536 

577 

• 3253 

3,4836 

b2.K 

. 1847 

3,6'+81 

578 

• 2243 

3,9457 

b29 

•2983 

3.7990 

580 

.4323 

3,5026 

b30 

.5388 

2.7415 

581 

• 3865 

3.6278 

b32 

.0715 

3 , 7033 

583 

.3436 

3.6203 

533 

.1201 

3.1972 

584 

• 3223 

3.4526 

53'+ 

.2481 

3.1434 

585 

• 1240 

3,6160 

53b 

.5632 

3,6784 

586 

• 36 8 8 

3.637 4 

536 

.5896 

3.6374 

588 

• 2345 

3.9257 

537 

.2399 

3.7466 

589 

.2110 

3.6411 

538 

• 3963 

3 . 8850 

590 

.3756 

3,2890 

54C 

.3539 

4,1497 

591 

• 3687 

3,3935 

54 \ 

.2059 

4,0133 

592 

• 2288 

4.2'965 

54 8 

.3146 

3,4315 

593 

• 3858 

3.5786 

54 3 

.2150 

3.1389 

594 

• 2403 

3,3953 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 1 - Continued 


SCAN 

CO 

CH 4 

SCAN 

CO 

CH 4 

NO. 

(atm- cm) 

(atm -cm) 

NO. 

(atm-cm) 

(atm-cm) 


• 3420 

3 . 0 4 9 0 

645 

• 1698 

3.4495 

by? 

• 4998 

4.017? 

646 

♦ 2159 

3.9573 

b98 

• 157b 

3.D603 

647 

• 2975 

3 . 9574 

byq 

• 424 7 

3.6813 

648 

•389 7 

4 . 1862 

600 

• 1871 

4.0481 

649 

• 3530 

3, 7326 

601 

•4109 

2.8514 

65c 

• 2701 

3,4851 

60? 

• 07 38 

3 • 7 0 0 7 

65? 

• 4306 

4,3353 

604 

• 2770 

4,3084 

653 

• 2577 

3.9929 

60b 

• 1902 

3.338? 

654 

• 36 38 

3,3630 

606 

• 3259 

3 • o?06 

655 

• 1426 

3.6467 

607 

• 6 ? 6 S 

3.4115 

656 

* 1945 

3 , 04 7 1 

608 

• 2454 

3,9362 

657 

• 1821 

4 . 4800 

609 

• 6068 

3 » b 6 8 7 

658 

.0765 

3.9314 

610 

• 2014 

3,4794 

66 c 

• 0952 

3.4958 

61? 

• 0636 

3 , 9554 

66 1 

•2014 

3.1061 

619 

• 1637 

3.5953 

662 

•1318 

4.0710 

614 

• 1598 

3,3301 

663 

• 1806 

3.4139 

615 

• 4952 

?,8<+75 

664 

-•0237 

4,0238 

616 

•2567 

3,1798 

665 

•3453 

4,1098 

617 

• 1964 

4,U099 

666 

• 376 7 

3,7916 

618 

• 2597 

3,3153 

668 

•1223 

4, 1844 

6?0 

• 3280 

4,0318 

669 

• 3660 

3,0951 

621 

• 1745 

3.6790 

67o 

• 3890 

3,7040 

6 ?? 

•2370 

3 , 08 48 

671 

• 4051 

3,6560 

6 ? 3 

• 2564 

3,b8o7 

67? 

• 3410 

3*2676 

6?4 

• 4429 

3.2533 

673 

• 2633 

3,3534 

625 

• 2260 

3,4534 

674 

• 2179 

4, 1?37 

626 

• 57 85 

3.5117 

6 76 

• 2404 

3,688.8 

628 

• 0345 

^•120-3 

677 

•1029 

3 , 0505 

629 

•0911 

4,2961 

678 

•3956 

3*6307 

63o ■ 

• 3291 

3.8373 

679 

• 3354 

4,6841 

631 

• 5265 

3,7q49 

680 

• 2938 

3,1424 

63? 

• 2200 

3, 681 5 

681 

. •1257 

3,3827 

633 

• 3670 

4 , 3?42 

682 

• 1061 

3,7920 

634 

• 3427 

3 , 3690 

634 

• 2738 

3.7112 

636 

• 0978 

3,7733 

685 

*2990 

3,6155 

637 

• 4875 

3,2260 

686 

• 1042 

4.3178 

638 

•3020 

3,5784 

687 

• 4345 

3,0503 

639 

•3368 

4,1089 

688 

• 5124 

3 , 67 62 

640 

.5475 

4,024? 

689 

• 0850 

4.1894 

64 1 

• 1255 

3,7961 

69o 

*3401 

2.9083 

64? 

• 4205 

3,5003 

692 

• 2436 

3,8055 

644 

•2406 

3.7924 

693 

• 3398 

3.3031 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS I - Continued 


SCAN 

CO 

CH4 

SCAN 

CO 

CH4 

NO. 

(atm- cm) 

(atm -cm) 

NO. 

(atm -cm) 

(atm- cm) 

bSh 

* 30 76 

4,3462 

744 

• 4508 

3.7579 


• 3713 

3.2095 

745 

. 1594 

3.954 7 

636 

• 4556 

4,0562 

746 

-,122b 

3,4149 

637 

• 2S5» 

3,8356 

74 8 

.0680 

3.9513 

63 H 

• 1800 

3.2073 

749 

.4829 

3,9301 

700 

• 3621 

4,0033 

750 

. 1912 

4,1199 

701 

• 2143 

4,0279 

751 

• 3706 

3,5468 

70? 

• 2713 

4,0246 

752 

.4301 

4,0?45 

703 

• 1784 

4,4967 

753 

.4228 

3. 74?0 

704 

• 3720 

3.6151 

754 

.2892 

3,2673 

7 05 

• 1834 

3,6815 

756 

.2149 

4.304 4 

706 

- .0569 

4,1049 

757 

• 0567 

4,5085 

70S 

• 6185 

3 , 3 0 3 8 

758 

• 0615 

3,3264 

709 

• 1904 

4.3701 

759 

.2124 

3.2552 

710 

• 1444 

3,9782 

76o 

• 02 73 

3,5825 

711 

.3881 

3,3063 

761 

- . 1602 

2.7442 

71? 

.3480 

4, 0239 

762 

.4504 

2.6617 

713 

.0640 

. 3.3514 

764 

.2331 

3,6278 

714 

.3616 

3,4333 

765 

.3066 

3,6287 

716 

.0595 

4,2171 

766 

.2902 

3,9451 

717 

• 1433 

3,4260 

767 

*.1197 

3.5703 

718 

.2130 

3,4244 

768 

- .0771 

3,9907 

719 

.2451 

3,9875 

769 

• 1431 

3.5653 

720 

.5476 

3,7732 

770 

.2710 

3.0198 

721 

. 1647 

3.9809 

772 

.7107 

3.5736 

7?2 

.2888 

4,3687 

773 

.0825 

3,7305 

724 

.3865 

3,8327 

774 

.1034 

4,6591 

725 

.3539 

4,1933 

775 

.3308 

3,4203 

726 

.0724 

4,7648 

776 

*3864 

3.7169 

727 

.1117 

3.6538 

777 

.0950 

4,3656 

728 

.3261 

3.2113 

778 

.2872 

4.1862 

729 

.3331 

3,7831 

78o 

• 194b 

3,7804 

730 

.2872 

3.8079 

781 

.0562 

3,5782 

732 

• I860 

4.1111 

782 

. 1747 

3,1294 

733 

.2303 

4,3677 

783 

. 1416 

4,3846 

734 

.2011 

3.5780 

784 

-•0149 

4,2199 

735 

.7204 

3.3603 

785 

• 2157 

3,6093 

736 

• 0884 

3.2708 

786 

.0437 

3,9714 

738 

.?493 

4,4892 

788 

• 1421 

3,7192 

74 0 

• 062 1 

3,6423 

789 

. 1959 

3.6568 

74 1 

.3798 

3 . 824 7 

790 

.3478 

3.9927 

74? 

• 0702 

3.7184 

791 

. 1979 

3,9617 

74 3 

.3426 

3.1718 

792 

.3212 

3,9463 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 1 - Concluded 


SCAN 

CO 

CH4 

SCAN 

CO 

CH4 

NO. 

(atm -cm) 

(atm- cm) 

NO. 

(atm -cm) 

(atm- cm) 


793 

”•1623 

3,8937 

842 

« 3565 

3. o860 

794 

.4649 

3,4827 

844 

,2006 

3 • 0 7 0 6 

796 

• 1711 

4,1553 

845 

, 1495 

4,4679 

797 

.5779 

3,4853 

846 

.5773 

2.9689 

798 

• 2556 

4, 7868 

847 

.1363 

3,o891 

799 

• 1145 

3,4703 

848 

,6451 

2,9528 

aoo 

• 2770 

3,9799 

849 

,1011 

3.0967 

801 

- .1993 

4.1181 

850 

.2870 

4,6828 

80? 

.0291 

3,8873 

852 

.1108 

?.9?07 

804 

*1867 

3.9052 

853 

.4763 

4,6941 

805 

“ *2203 

3.8315 

854 

.3277 

3,5806 

806 

• 2408 

3,8482 

855 

.3452 

3.7967 

807 

» 2968 

3.2358 

856 

.3588 

2,6560 

80S 

• 1217 

4.3579 

857 

• 1715 

4,7654 

809 

.2357 

4 • 6650 

858 

• 3163 

3 , 6580 

810 

• 2469 

3.7579 

860 

.1646 

3.6162 

812 

.4798 

4, 1364 

861 

,3266 

3,3507 

813 

.3800 , 

2,8540 

862 

.0192 

4,2168 

814 

. 1409 

3.8365 

863 

,3085 

3.7707 

815 

-.0038 

4.5899 

864 

-.0857 

4,3306 

816 

.0569 

4.1330 

865 

.6671 

3,7557 

817 

.1821 

3. 7697 

866 

.3210 

3,6373 

818 

.1022 

3,5541 

868 

• 0924 

3,5251 

820 

.4318 

3.9949 

869 

,1156 

4,4565 

82l 

.5142 

3 • o458 

87o 

•0735 

4,0488 

822 

.2843 

4,4393 

871 

• 1672 

3,9825 

823 

.2944 

3,7261 

»72 

.4748 

3.2381 

824 

.3971 

3,6872 

873 

.2498 

3,4621 

825 

. 0944 

4,2770 

874 

, 1974 

4,2047 

826 

.6754 

3.0305 

876 

,4012 

3,7691 

828 

.4526 

3,9432 

878 

,3086 

2,6420 

829 

.0879 

4 . 2668 

879 

. 3653 

3.1910 

830 

* 0456 

3.9784 

880 

*4595 

4,1029 

831 

.4843^ 

3,9239 

882 

,5578 

4,6950 

832 

.0996 

3,4931 

883 

,3745 

3.6817 

833 

.3206 

3,0610 

884 

,2809 

3,1809 

834 

.2953 

3.7912 

885 

.1022 

4,0699 

836 

• 0359 

3,8679 

886 

• 1786 

3,7063 

837 

-.0350 

3,3926 

887 

.3628 

3,4542 

838 

.2541 

3, 7200 

888 

.2206 

3,6038 

839 

.3229 

2,4163 




840 

.2155 

3,7555 




841 

.0374 

3,4919 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 2 


SCAN 

CO 

CH 4 

NO. 

(atm-cm) 

(atm-cm) 


100 b 

#2868 

4.1674 

1006 

.4163 

3,0810 

1008 

#2698 

3,5252 

1009 

#4106 

4.U497 

1010 

#4125 

3.2554 

1011 

#2376 

2.8463 

1012 

#1571 

3,6803 

1013 

#3713 

3,8110 

lOlA 

#4311 

4 ,0535 

1016 

#1959 

3,7383 

1017 

#3348 

3,1565 

1018 

#2712 

3.1903 

1019 

#3131 

3,4222 

10^0 

# 556 7 

4.2075 

1021 

#2495 

3,5897 

1022 

#1119 

4,0217 

1024 

#3633 

3,8737 

102 b 

#5056 

4,0399 

1026 

.2043 

3.7509 

1027 

#0512 

4,2556 

1028 

#2602 

3,8406 

1029 

# 1773 

4.3095 

1030 

# 4623 

3,0569 

1032 

#0762 

3,9604 

1033 

#3713 

3,5240 

1034 

#4023 

3,6930 

1035 

#2207 

3,1706 

1036 

#3655 

3,8404 

1037 

#2895 

4.2679 

1038 

#0719 

3,3185 

1040 

#1879 

3,7050 

1041 

#4165 

3,9070 

1042 

#3008 

3,4111 

1043 

#3597 

3,8848 

1044 

#2708 

4.0655 

1045 

-#0734 

3,6664 

1046 

#1542 

3,4808 

1048 

#6799 

3,0064 

1049 

#3843 

3 , 6562 

1050 

#2082 

4.4672 

1051 

#505 7 

3,7472 

1052 

#2054 

3,4425 

1053 

#4103 

3.7695 


SCAN 

CO 

CH 4 

NO. 

(atm-cm) 

(atm-cm) 


1054 

• 2282 

4,4830 

1056 

#1722 

3,2463 

1057 

#0227 

3.7654 

1058 

#4300 

3,6119 

1059 

.3394 

3,5375 

1060 

#1983 

4,9188 

1061 

#3380 

4,2097 

1062 

#4366 

4.1020 

1064 

# 1444 

4.2077 

1065 

#2345 

3.2104 

1066 

#4154 

3,4784 

1067 

-#0278 

3.5251 

1068 

#2004 

4,0226 

1069 

#2713 

3,4799 

1070 

#4688 

2,7993 

107? 

#2308 

3,5638 

1073 

#2481 

3,2826 

1074 

#3383 

3.6115 

1075 

#3795 

3.6156 

1076 

#0610 

5.0027 

1077 

#2162 

4 • 3663 

1078 

#6051 

3.6921 

1080 

#4799 

2,9864 

1081 

#2043 

3.4243 

1082 

#2707 

3,4847 

1083 

#2524 

3. 7960 

1084 

#3950 

3,2422 

1085 

.2688 

3,6131 

1086 

#2864 

3.7740 

1088 

.3211 

4,3087 

1089 

.1156 

3,6293 

1090 

.0618 

3,5096 

1091 

#0952 

3.4909 

1092 

#1608 

4.3698 

1093 

.2489 

3,9753 

1094 

.3202 

3,4073 

1096 

#2521 

3.6283 

1097 

#2940 

3,7519 

1098 

. 0685 

3,1820 

1099 

#4139 

3.7106 

1100 

#2434 

3,6845 

1101 

-#0613 

3,4410 

1102 

.3622 

3.9003 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 2 - Concluded 


SCAN 

CO 

CH 4 

NO. 

(atm-cm) 

(atm-cm) 


IlOA 

-^0384 

4.3229 

1105 

.2016 

3,7565 

1106 

.2451 

4.0009 

1107 

.3988 

3,3109 

1108 

.3941 

4,6328 

1 109 

.1494 

4,0542 

1110 

. 1759 

3.9237 

111? 

.2208 

4.2733 

1113 

. 1911 

2,9248 

1114 

.5337 

3,1117 

1115 

.0853 

3,7637 

1116 

.2206 

4.1741 

1117 

.2121 

3,5599 

1118 

.3904 

3,9989 

1120 

.0752 

3,4651 

1121 

.2110 

3, 0374 

1122 

.3609 

3.7407 

1123 

.5880 

3.7233 

1124 

.2753 

3.0305 

1125 

.0788 

2,9862 

1126 

.6016 

3,8256 

1128 

.4166 

3.7921 

1129 

.1936 

3,8989 

1130 

.1040 

3,7903 

1131 

.1879 

3.8463 

1132 

.0515 

4,0103 

1133 

.5241 

3.5129 

1134 

.2733 

3,6744 

1136 

.3361 

3,8689 

1137 

.2747 

3,4437 

1138 

.2213 

3.2364 

1139 

.0939 

3,3710 

1140 

.3368 

3,3618 

1141 

.4734 

4,1248 

1142 

.4334 

4,6006 

1144 

. 1556 

4.0145 

1145 

.2519 

4,0727 

1146 

.2339 

3,4130 

1147 

.2511 

3.7755 

1148 

.3221 

3.5119 

1149 

.1358 

3,3880 

1150 

.1780 

3,4276 

1152 

.2006 

3.9702 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 3 


SCAN 

CO 

CH 4 

SCAN 

CO 

CH 4 

NO. 

Catm-cm) 

(atm-cm) 

NO. 

(atm-cm) 

(atm -cm) 


H77 

.4768 

3.3745 

1522 

.0873 

3.4279 

1^78 

.2175 

3.9028 

1523 

.3473 

3.3312 

1479 

.5268 

3 • 8084 

1525 

.1845 

4,2121 

1480 

.0510 

3.5577 

1526 

.2708 

3,5464 

1481 

.3010 

4.4073 

1527 

.1676 

4,4597 

1483 

.3408 

4.U858 

1528 

.0015 

3,5720 

1483 

.1356 

4,3624 

1529 

.4998 

3.8104 

1485 

.1366 

3,9919 

1530 

.3240 

3.4525 

I486 

.0420 

3.9280 

1531 

.1640 

3,5448 

1487 

.4430 

3,1901 

1533 

.2804 

4,6137 

1488 

.3065 

3,5991 

1534 

.2389 

3.9817 

1489 

.5330 

4.1067 

1535 

- .0078 

3,0343 

1490 

. 1082 

4,3858 

1536 

. 1 355 

4,1415 

1491 

.3474 

3,9523 

1537 

.1432 

3,5880 

1493 

. 1574 

3.6289 

1538 

.3692 

4.0141 

1494 

• 1642 

3.6896 

1539 

*2251 

3.6787 

1495 

.45HO 

3.7801 

1541 

.2245 

3.1454 

1496 

.4135 

3,7834 

1542 

i 1874 

3 . 8629 

1497 

. 0583 

3,7813 

1543 

.0626 

4,4396 

1498 

.1686 

3.8327 

1544 

.2862 

4,1641 

1499 

• 1752 

3.6347 

1545 

. 2066 

3,6483 

IbOl 

.2093 

3,1881 

1546 

.2179 

3,6191 

lb02 

.3324 

3,8098 

1547 

.1164 

3,3416 

1503 

.2178 

3.5156 

1549 

»162b 

3.6864 

1504 

.2450 

4.0570 

1550 

.1110 

4,2651 

1505 

.0790 

3,1762 

1551 

• 4476 

4 .1902 

1506 

.2745 

4,7297 

1552 

.5530 

3.5514 

1507 

.6113 

4.0112 

1553 

.1712 

3.5763 

1509 

.4603 

3,5867 

1554 

.3242 

3,7022 

1510 

.5168 

3.6849 

1555 

.2087 

3,5765 

1511 

.2081 

3.5694 

1557 

.2237 

3,6872 

1512 

.2314 

3 . 5827 

1658 

.4336 

3,8788 

1513 

.2073 

3.9886 

1559 

.3563 

4.2715 

1514 

• 0252 

3.9828 

1560 

.4939 

3.4884 

1515 

.2130 

3,9023 

1561 

.4564 

3,6975 

1517 

.3302 

2.7942 

1562 

.3355 

3.8744 

1518 

.2266 

3.fcj323 

1563 

.1661 

4 . 0055 

1519 

.3755 

3,7126 

1565 

• 3041 

3,3584 

1520 

• 0644 

4,4009 

1567 

.1825 

3,8688 

1521 

.4109 

3.8127 

1568 

.2839 

4,1946 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 3 - Continued 


SCAN 

CO 

CH 4 

SCAN 

CO 

CH 4 

NO. 

(atm-cm) 

(atm-cm) 

NO. 

(atm-cm) 

(atm-cm) 

1^69 

.4193 

3.2297 

1616 

.3507 

3.9522 

1570 

.4676 

4.0736 

1617 

.2929 

3.385 7 

1571 

.5647 

3.3441 

1618 

.3855 

3.9930 

1573 

.4298 

3.8105 

1619 

.1323 

3.9168 

1 57 4 

.3105 

4,1?81 

1621 

. 3549 

3,b629 

157b 

.2955 

3,8639 

1622 

. 3696 

3 . 3663 

1576 

-.1501 

4,0385 

1623 

.2620 

3, 7535 

1577 

. 1515 

4.0345 

1624 

• 2991 

3,6925 

1578 

.2033 

3,5?96 

1625 

.4078 

4,0020 

15 79 

.2469 

3.6731 

1626 

.2930 

4,0952 

1581 

•3189 

3,5239 

1627 

.4781 

3,4514 

158? 

.3703 

3.1980 

1629 

.2214 

3.8439 

1583 

. 1622 

3,8813 

1630 

.3075 

3 , bb44 

1584 

.2139 

4,5345 

1631 

.4059 

3,9219 

1585 

• 1871 

4,6585 

1632 

.3370 

3.5830 

1586 

.3040 

3.4062 

1633 

.3876 

3.5376 

1587 

• 2533 

4,4232 

1634 

.2906 

4.4530 

1589 

• 2944 

3.5255 

1635 

.3088 

4,0479 

1590 

.2776 

3.4718 

1637 

• 1345 

3.7847 

1592 

.2295 

3.5109 

1638 

.3436 

3.2425 

1593 

.3083 

3,5508 

1639 

• 2335 

3.O080 

1594 

.2588 

3.8130 

1640 

.3038 

3.5825 

1595 

. 1553 

4.2368 

1641 

• 2734 

4.2928 

1597 

.2526 

3.55?9 

164? 

.1000 

3,3930 

1598 

.2053 

3,8751 

1643 

.2067 

4.1527 

1599 

.8513 

3.5342 

1645 

• 3409 

3,5819 

160,0 

.0728 

3.5339 

1646 

• 1387 

4,0837 

1601 

• 4o43 

3,8289 

1647 

.4962 

4.0801 

160? 

.4416 

3,5767 

1648 

• 4216 

3,7687 

1603 

.3673 

3.3965 

1649 

.1031 

3 , 0630 

1605 

• 1656 

3.3766 

1650 

.5030 

4,0007 

1606 

.3975 

3.2749 

1651 

.4021 

4,1301 

1607 

• 1308 

3.2493 

1653 

• 1253 

3,0392 

1608 

• 2440 

3.967? 

1654 

• 5550 

3.9279 

1609 

• 3484 

3.4985 

1655 

.3558 

4.3770 

1610 

. 1883 

4 . 3553 

1656 

• 1539 

3,5717 

1611 

.2395 

3,8755 

1657 

.3974 

3.1919 

1613 

.5238 

3,6800 

1658 

.2791 

3.0581 

1614 

.2722 

4.2438 

1659 

.3003 

3.2392 

1615 

.2565 

3.5113 

1661 

.3284 

3,9508 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 3 - Continued 


SCAN 

CO 

CH 4 

SCAN 

CO 

CH 4 

NO. 

(atm-cm) 

(atm-cm) 

NO. 

(atm-cm) 

(atm-cm) 

1662 

• 20b 7 

3,2993 

1 707 

,2418 

3.9159 

1663 

.4011 

4.4723 

1709 

• 0265 

4,0695 

1 664 

• 2549 

3,9832 

1710 

• 1649 

3. 3324 

166b 

.3413 

3,8486 

1711 

.6964 

4.1459 

1666 

.2601 

3,8580 

1712 

, 1956 

3 , 6005 

166 7 

.2563 

4,0266 

1713 

• 1561 

3,3155 

1669 

.3478 

3.1449 

1714 

' *2875 

4.2316 

1670 

• 1428 

3.3879 

1715 

,3306 

4,0059 

1671 

.4092 

3,8860 

1717 

,3304 

3,6309 

167? 

.2615 

4.3452 

1 718 

,2696 

3,4067 

16 73 

.0493 

3 , 8809 

1719 

,2311 

3,5805 

1674 

. 1587 

3,6117 

1720 

, 1524 

4.1449 

167b 

.1909 

3.1674 

1721 

.14 79 

3.7665 

1677 

• 2223 

3.6277 

1722 

,2902 

3,2383 

1678 

.3434 

3.6318 

1723 

,3526 

4.3002 

1679 

.124? 

4.1511 

1725 

• 1315 

3 , 6202 

1680 

.4242 

4,1712 

1 726 

,4002 

3,3151 

1681 

.1966 

4,0987 

1727 

,0820 

3 , 0 b 4 6 

1682 

.2446 

4,0355 

1728 

,3533 

3,4064 

1683 

.3025 

3.9236 

1729 

,5358 

3,7462 

1685 

.3292 

3, J373 

1 730 

. 3063 

3 . 3754 

1686 

.2191 

4,1675 

1731 

• 1196 

4,1177 

1687 

.3927 

3,8789 

1733 

,4239 

3.3292 

1688 

.2069 

3,2807 

1734 

• 0977 

3.8958 

1689 

.1119 

3,4802 

1735 

.3787 

3,9266 

1690 

.3089 

3,9911 

1736 

.5194 

4,1238 

1691 

.2975 

3.2313 

1737 

.0385 

3.8903 

1693 

.0465 

3,2169 

1738 

.2405 

3*8801 

1694 

.1650 

4.1514 

1739 

.0694 

3.9135 

1695 

.3675 

3,5476 

1741 

.1990 

3,2134 

1696 

.1479 

3,7486 

1742 

.3508 

3.9128 

1697 

.3809 

4.1470 

1743 

.0186 

3,8950 

1698 

.239 7 

3,6515 

1744 

.3966 

3.5590 

1699 

.1737 

3.9723 

1745 

.0070 

4.8462 

1701 

.2938 

3,8524 

1746 

. 1026 

3,6716 

170? 

• 2279 

3,8987 

1747 

. 3544 

3,7890 

1 703 

.2377 

3,5570 

1749 

.3167 

3,7640 

1704 

.520? 

3,356? 

1750 

.1184 

3.5218 

170b 

- .0024 

4,2962 

1/51 

.2873 

3.9571 

1 706 

.0543 

3,8803 

175? 

. 1813 

4,1515 
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CHURCHILL AREA SINGLE SCAN DATA 
PASS 3 - Concluded 


SCAN 

CO 

CH^ 

NO. 

(atm -cm) 

(atm- cm) 


1 /b3 

.5271 

3,6298 

1 754 

.2216 

3,7671 

1 /55 

.1377 

3.3227 

1 /57 

.0433 

3,2476 

1 /5S 

. 1901 

3,8724 

1 759 

.2451 

3,8785 

1/60 

. 4655 

2.9977 

1761 

.2384 

3,1290 

176? 

• 0368 

3,9787 

176 3 

- . 0076 

3,7299 

1765 

.2749 

3,8975 

1 /66 

• 4816 

4,0721 

1767 

.2008 

3.6595 

1768 

. 1766 

3.0828 

1769 

•2061 

^.1099 

177c 

.3258 

'+.3265 

1771 

.4176 

3,7880 

1773 

.4276 

4.0791 

1774 

.4462 

3.7716 

1775 

• 1772 

3.5431 

1776 

•0023 

3.9505 

1777 

• 3958 

4.3587 

1 778 

• 3135 

3,7452 

1779 

• 5307 

3.5998 

1781 

.4194 

3.6877 

1782 

.3288 

3,6862 

1783 

.3251 

4.0213 

1 784 

.5216 

3,7277 

1785 

. 1629 

3.2293 

1786 

. 1498 

3.4528 

1787 

.4086 

3.4054 

1789 

• 1951 

4,1959 

1790 

.5041 

3,3472 

i/yi 

.3898 

3.‘^420 

1792 

.4213 

3,9977 

1 793 

• 2694 

4,3403 

1794 

.3919 

3.6098 

1795 

.4340 

3.8411 

1797 

• 1741 

3,3817 

1798 

.3412 

3,1748 

1799 

• 0135 

3.9825 

1800 

.7353 

3.9083 
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